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An Analytical Approach to the Fuel Sloshing 
and Buffeting Problems of Aircraft’ 


HAROLD LUSKIN¢ | 
Douglas Aircraft 


ABSTRACT 


Two dynamics problems are discussed: the motion of an air- 
craft with partially full fuel tanks and airplane buffeting. 

In discussing the first problem, it is stated that for small mo- 
tions the fuel in a tank can be replaced by a fixed mass plus a 
simple pendulum. A spring-mass system is an alternative that is 
used in writing the linearized equations of longitudinal motion 
for an aircraft with a partially full fuel tank. The motion is found 
to be described by three essentially independent oscillations: 
the phugoid, the short period, and a fuel oscillation. The fuel 
oscillation has a frequency essentially the same as that of the fuel 
in a stationary tank, and significant damping is obtained only 
by the use of tank baffles. 

In connection with airplane buffeting, it is indicated that it 
appears to be an aeroelastic response problem that should be sus- 
ceptible to analytic approaches. Some ideas on the nature of the 
forcing functions are presented. 
from turbulence and the mechanics of stochastic processes may 


be useful in developing a rational buffeting theory. 


It is suggested that concepts 


INTRODUCTION 


ay AERONAUTICAL ENGINEERS must solve 
many problems without waiting for suitable 
theories to be developed, it is nevertheless true that 
theoretical methods may be of considerable assistance. 
Theory may indicate the correct parameters and some 
of the relations between them. Analysis may suggest 
useful experiments and the proper quantities 
measured. Methods for interpreting the experimental 
data may be found. A successful theory may even re- 
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sult. In such cases certain steps have often been fol- 
lowed. 

(1) The physics of the problem has been understood 
and the problem has been formulated. 

(2) The equations (simplified to the point where 
they can be handled easily but not to the point where 
the essential character of the problem has been lost) 
relating the variables have been written. 

(3) The solution to the equations have been found 
and applied. 

(4) Experimental results have been obtained to 
verify the accuracy of the assumptions. 

As the engineering art has developed, many problems 
have passed through these steps. The problems not 
solved by analysis might have resisted attack for many 
Some of them might never have been im- 
portant. For others, the mathematics may not have 
been available, or the labor of obtaining the numerical 
Sometimes the 


reasons. 


solutions may have been too great. 
problem was not well defined. 

Two aeronautical problems for which complete the- 
ories do not exist are considered here. The problem 
of the motions of an aircraft with partially full fuel tanks 
has not been sufficiently pressing in the past to attract 
The buffeting problem has 
It is the intent 


much analytical attention. 
not yet been completely formulated. 
of this paper to carry these problems through some of 
the steps toward solution. 

The problem of the motion of an aircraft with par- 
tially full fuel tanks (also called ‘‘fuel sloshing’’) has 
been known for many years. Test pilots have been 
aware of it because it was associated with both lateral 
and longitudinal oscillations of disturbing magnitude 
and with unusual fore-and-aft trim changes. It has 
been eliminated in the few cases where it has occurred 
by the installation of empirically designed tank baffles. 


Turbojet and rocket-powered vehicles, however, have 
generally tended to have higher than usual values of 
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the ratio of fuel weight to gross weight. Also, the trend 
to higher wing loadings and thinner wings has removed 
the tanks from the wings, where structural members 
acted as baffles, to the body, where large uninterrupted 
volumes for tanks exist. Recently, the existence of 
unexplained flight oscillations has served to focus at- 
tention on sources of the difficulties. All of these fac- 
tors have indicated the need for analysis of the prob- 
lem. Work on this problem at the Santa Monica Divi- 
sion of the Douglas Aircraft Company, Inc., began in 
1949. The first published analytical contribution!’ 
to the literature of sloshing was made in January, 1951. 

The beginning stage in the Douglas study of fuel 
sloshing developed the concept of replacing the oscil- 
lating fuel in each tank by a simple pendulum plus a 
fixed mass. It was shown by Graham? that this con- 
cept was consistent with the linearized equations for the 
small motions of a fluid in a stationary rectangular 
tank. He later extended the analysis to cover a variety 
of other tank shapes. This work was of primary im- 
portance, since the motion of the fuel and the forces it 
transmitted to the aircraft were the only elements of 
the problem not previously well known from conven- 
tional stability analyses. The aircraft with sloshing 
fuel was then replaced by one with a simple pendulum 
and a fixed mass for each tank. Some of the results 
of this work having to do with the longitudinal motion 
are taken up here. The case of the lateral motion has 
been discussed in reference 1. 

Aircraft buffeting is a problem that has been little 
understood. In fact, a recent review of aeroelas- 
ticity’ had only one statement to make concerning 
buffeting: ‘“‘The aerodynamics underlying the buffet- 
ing phenomena . is still fairly obscure.’’ Buffeting 
first became important when the low-wing monoplane 
was introduced. Tests‘ of that configuration showed 
that a generally satisfactory method for reducing the 
buffeting to tolerable levels was the use of the wing- 
fuselage fillet. Buffeting research was largely directed, 
therefore, at development of fillets for each particular 
aircraft. Such additional work® as was done con- 
sidered the vibrations to occur in the following way: 
‘A stalled wing produces a wake of considerable ex- 
tent and intensity. The wake is essentially a 
Karman vortex trail and thus has a unique frequency 
associated with it. The tail structure responds to the 
fluctuations and is supposed to have its own natural 
frequencies removed from that of the wake in order to 
avoid resonance.’ Other research into buffeting has 
been exploratory (or has merely defined the flight 
buffeting régimes for particular aircraft) and has not 
resulted in a clear formulation of the problem. This 
being the case, it has not been possible to proceed with 
a solution. Furthermore, new types of buffeting have 
appeared which are not corrected by fillets or described 
by older theories. At the present time, developments 
in the fields of dynamic analysis and fluid turbulence 
may be possibly be of aid in providing a better under- 
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standing of the buffeting problem. 


paper this possibility is briefly explored. 


In the present 


Study of the Literature led to some calculations py 
Graham® on the disturbances created by a vortex trai 
Later, Liepmann’ suggested that some of the ideas of 
turbulence could be used to provide insight into the 
aerodynamics of buffeting. He used this to calculate 
the mean square lift response of an airfoil to the random 
angle of attack fluctuations of isotropic turbulence, 
since a turbulent wake is considered to be more realistic 
than the vortex trail. 


NOMENCLATURE IN ORDER OF APPEARANCE 
W, = weight of equivalent simple pendulum 
which replaces the fuel in the equiva. 

lent dynamical system, lbs. 


W, = actual fuel weight, lbs. 

h/l = ratio of fuel height to fuel-tank length 

L,/t = ratio of simple pendulum length to 
fuel-tank length 

Y,/i = ratio of pendulum height above un 
disturbed fluid surface to fuel-tank 
length 

= period, sec. 

g = gravitational acceleration constant, ft. 
sec.” 

F = horizontal force on fuel-tank sides, lbs 

0 = pendulum displacement angle; also 
fuel free surface displacement angle, 
rad. 

W, = total weight, airplane plus fuel, lbs. 

my = total mass, airplane plus fuel, slugs 

My = mass of equivalent pendulum, slugs 

u, Ww = perturbation velocities of the aircraft 
along the x, z axes, ft. per sec. 

q = angular velocity of the aircraft about 
the y axis, rad. per sec. 

Xaero., Zaero., Maero. = aerodynamic forces parallel to the x 
and z axes, and aerodynamic mo- 
ment about the y axis, lbs. and lb. ft 

6 ‘ = angle between the x axis and the hori- 
zon, rad. 

Jo = initial value of 6, rad. 

i] = perturbation in 6, rad. 

Uo = air speed in undisturbed flight, ft. per 
sec 

V = air speed at any instant, ft. per sec. 

Ler % = coordinates of the airplane plus fixed 
mass of fuel center of gravity rela- 
tive to the x, z axes, ft. 

eS = coordinates of equivalent moving fuel 
mass center of gravity relative to 
the x, z axes, ft. 

L = lift, lbs. 

s = x coordinate of the equilibrium posi- 
tion of the moving fuel mass center 
of gravity in undisturbed flight, ft 

I, = moment of inertia of the system, ait- 
plane plus all fuel, about the y axis 
slug-ft.? 

Fa = moment of inertia of the airplane plus 


the fixed fuel about an axis through 
the point (Xa, 2a), slug-ft.? 

rs = moment of inertia of the moving fuel 
about an axis through the point 
(Ep, Zp), slug-ft.? 
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FUEL SLOSHING AND 


= radius of gyration, ft. 


: = mean aerodynamic chord, ft. 

k = spring constant of the equivalent fuel 
mass-spring system, lbs. per ft. 

D: = airplane rotary damping constant, lb. 
ft. sec. 

Dy = fuel viscous damping constant, lb. sec. 

per ft. 

N = root of the characteristic equation 

7. = (OXnero./OU), 0 

a = (OXuaero./OW)w 0 

a = (OZnero. /OU), =O = aerodynamic 

a = (OZnero./OW)w =o resistance 

Ma = (OMacero./Ow/ Uo)w =0 derivatives 

M, = [OMsero./O(gc/Uo)| = aerodynamic 
rotary derivative 

D = drag, lbs. 

Cp = drag coefficient 

Cbp = parasite drag coefficient 

Ch = lift coefficient 

“ = lift curve slope 

a = angle of attack, rad. 

Vy = vertical velocity amplitude, ft. per sec. 

I = circulation, ft.?/sec. 

b = spacing between vortices, ft. 

h = spacing of vortex rows, ft. 

Ci? = mean square lift coefficient 

w? = mean square vertical velocity, ft.?/sec.? 

( = airfoil chord, ft. 

L = scale of turbulence, ft. 


AIRCRAFT WITH PARTIALLY FULL 


FuEL TANKS 


DYNAMICS OF AN 


The Dynamics of the Fuel Motion 


The one new component occurring in the problem 
being considered is the fuel, which is free to move. 
As many new degrees of freedom are thereby added to 
the ordinary rigid body aircraft dynamic analysis as 
there are coordinates describing the fuel configuration. 
For a completely general fuel motion, an infinite num- 
ber of coordinates are needed. Such a description of the 
dynamics is complex, and only the configuration for 
the fuel moving under small disturbances (which is con- 
sistent with the assumptions usually made in aircraft 
dynamic stability analysis) in its fundamental mode 
will be studied here. With this assumption each tank 
adds two new degrees of freedom to the dynamical 
system. Considering longitudinal motions only, a 
tank adds only one new degree of freedom. 

Insight into the dynamics of the fuel motion can be 
had by observing the motion of water in a container. 
An arbitrary large disturbance to the container will 
create a confused motion of the free surface, but soon 
the smaller waves will die out, leaving only a moder- 
ately damped more or less periodic motion. Repeating 
the experiment with a smaller disturbance along only 
one axis will show the periodic motion more clearly. 
Finally, the use of a small container will show 
heavy damping and thus indicate light damping in 
large tanks. A similar example of fluid in motion is 
that of the manometer, which is easily shown to have a 


simple oscillatory mode.*® 
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[i 
~ 
8 = 
4 _— 
oe — 
~ 
i 
Fic. 1. Definition of fluid surface displacement angle and 


pendulum displacement angle. 


It is indicated that the fuel in a large aircraft tank 
without baffles will, after a small disturbance, have a 
lightly damped periodic oscillation, and that it can be 
approximated by a simple pendulum. This possibility 
was studied by Graham.” His results 
tionary? rectangular tank are based on the linearized 
expression for the velocity potential given by Lamb. 
He shows that the fuel in its fundamental mode can 
be replaced by a simple pendulum plus a fixed mass for 
pendulum angular displacements, 6, like the fluid sur- 
Sinu- 


for the sta- 


face angular displacement as shown in Fig. 1. 
soidal waves are found to exist on the free fluid sur- 
face. The expressions for the ratio of pendulum weight 
to fuel weight, pendulum length to tank length, and 
height of the pendulum axis above the undisturbed 
surface to the tank length are found to be functions 
only of the ratio of undisturbed fluid height to tank 


length. 
W,/W, = [tanh (rh//)|/(rh/1) (1) 
Lg/¥ = (1/m) coth (rh//) (2) 
Y,/l = (1/m) [2 — cosh (rh//)|/sinh (rh/l) (3) 


These relations are plotted in Fig. 2. 


Some intuitive feeling for the equivalent pendulum 
can be gained by studying the extreme cases of an al- 
most empty tank and a deep tank. For the nearly 
empty tank, the pendulum has nearly the same weight 
For the deep tank, on the other hand, the 
The geometry of the two 


as the fluid. 
effective mass is small. 
pendulumis is shown in Fig. 3. 


+A study” of the nonstationary tank has shown that the 
stationary tank is an adequate representation for this problem. 
t A tank that is completely full or nearly so and which has a 


fixed top may not be described by this analogy. 
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3 The period, P, of the fundamental mode is given by sai 
PENDULUM WEIGHT P/(2l V gh) = [(rh//) coth (xh/1) (4 
FUEL WEIGHT ; , ; 
aL while the amplitude of the horizontal force on the tank Iq 
sides, F, is given by 
nares oe ae . 
F/W, = [tanh (rh//) tan 6] /(rh/1) (x 
Fig. 4 is the plot of these relations. These formulas 
— simplify for cases where the fuel is shallow, h/] < |/9 Int 
0 n n J and where the fuel is deep, h// > 1/2. For example, att 
al 4 HEIGHT 190 a for the shallow tank the period is squ 
TANK LENGTH P = 2//gh (6 eat 
“lb POSITION OF PEND. AXIS a 
ABOVE FUEL SURFACE and for a deep tank 
TANK LENGTH , m,| 
Fic. 2. Characteristics of equivalent pendulum. P = 2V al/g (7 
Eq. (6) gives the time for a shallow water wave to make id 
a trip across the tank and back. The velocity of this My 
wave is V gh. 

The small motion approximation gives a vertical force In | 
on the pendulum which is equal to its weight and a side tem 
force which is proportional to its displacement. Those 
considerations show that a mass and a spring are com- I 
pletely equivalent to the pendulum in simulating the ” 
small motions of the fuel. The spring stiffness, the force the 
for unit displacement, is equal to the pendulum weight part 
divided by the pendulum length. The mass has the 

NEARLY sesctie anil aids ieee lt Do oe a 
same magnitude as that of the pendulum. ~ 

EMPTY DEEP of 1 
TANK TANK The Longitudinal Motion of an Aircraft with Sloshing Fuel ae 

The motion of an aircraft in its plane of symmetry is not 
to be studied for the case where fuel is in motion in one E 
tank. It is supposed that the fuel is replaced by its viel 

=, | —F>| equivalent mass and spring system. To eliminate a 

Fic. 3. Pendulum analogy to fuel in two-dimensional tanks. coupling with the lateral motion, the fuel will be as- nary 
Extreme cases. sumed to remain symmetrically disposed with respect and 

to the plane of symmetry. The 

A set of rectangular coordinates fixed relative to the orcs 
body is selected. The origin is assumed fixed at the the « 
position of the center of gravity of the entire system duce 
under the initial equilibrium conditions. The longi- or a 

2.0 tudinal, x, axis is assumed aligned with the initial rela- and 
tive wind, and, for simplicity, the motion of the fuel | fatio 
equivalent mass is taken as parallel to the aircraft It 

1.5 longitudinal axis. Fig. 5 illustrates the a indi 
system and its geometrical relation to the aircraft an ie 

RATIOS equivalent fuel ie centers of gravity. ao 

1.0 The system contains four degrees of freedom: trans- nea 
lation along the two coordinate axes lying in the plane the « 
of symmetry, a rotation about the axis perpendicular ine 

0.5 to the plane of symmetry, and a translation of the fuel Anti 
parallel to the aircraft longitudinal axis. By applica- and 
tion of Newton’s laws, the four independent differential ies 

oy 25 50 75 100 125 relations between these degrees of freedom may be de- quen 

FUEL HEIGHT rived in the following form: mode 

Fic. 4. Period and eine oy <a mode of m [te + w4q] a — . 1G 
a liquid in a two-dimensional rectangular tank. X aero. — MpXp — W, sin 0 + myxpq° (98) rectio 
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f 
M,| 0 


w— q(Uo + u) | = 
/ (Sb) 


Dios. W, cos 0 + myXpG + 2m,Xpq 
IG + MpFrXp — Mpl2(s — Xa — Xy)Gky + 
x,w + wx,] = Macro. — Wpx, cos @ (8c) 
— mt, = W, sin 6 + kx, + m,(s — Xa — Xp)q? + 
My,ZpG + m,(% + wq) (Sd) 
Introducing @ = 3% + 8, where the perturbation in 


attitude, J, is small, and neglecting products and 
squares of the perturbation variables, the linearized 
form of the equations of motion result. 

mili = Xeero. — MpX,p — W, sin do — Ww cos B (9a) 


= Lucro. + W, cos 3d) — WW sin 3) (9b) 


(9c) 


m(w — Uog) 


W Xp COS Jo 


Iq ~ My2pXp = |) 


m,i, = W, (sin do + 9 cos do) + Rxp + mp2 4- 


m,u (9d) 


In Eqs. (Sc) and (9c) the moment of inertia of the sys- 
tem has been taken as 


I, = Ig + Ip + ta(Xa? + 20”) + my(E,* + 2,7) (10) 


where J, is the moment of inertia for a rotation about 
the point (x,, 2.) of the aircraft including the “‘fixed”’ 
part of the fuel and J, is the moment of inertia for a 
rotation about the point (£,, 2,) of the equivalent mass 
of moving fuel, both parts of the fuel considered as 
rigid bodies. The correction to the moment of inertia 
that accounts for the fact that the fuel is not rigid has 
not been considered. ¢ 

Expansion of the matrix of coefficients of Eqs. (9) 
yields a sextic for the characteristic equation as con- 
trasted to the quartic that results from the ordi- 
nary aircraft stability equations [Eqs. (9a), (9b), 
and (9c) with equivalent fuel mass terms omitted]. 
The roots of the quartic equation generally describe the 
aircraft motion as composed of two oscillations. In 
the sextic the appearance of two additional roots intro- 
duces an additional mode, which may be an oscillation 
or a combination of subsidences and/or divergences 
and which also modifies the two original aircraft oscil- 
lations. 

It is desired that the nature of the new motion and the 
extent of its influence on the original motions of the air- 
craft be investigated. Ordinarily, the phugoid and 
short-period motions of an aircraft may be studied 
separately because, having widely different frequencies, 
the coupling between these motions is slight and the 
dynamic equations may be modified in accordance. 
Anticipating that the fuel motion natural frequency 
and the aircraft short-period frequency may be of the 
same order and thus far removed from the phugoid fre- 
quency, investigation of the nature of the additional 
mode is facilitated by consideration of the special case 

f Graham, in an extension of reference 2, has shown the cor- 
tection is small. 
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Geometry. 


Fic. 5. 


in which the long-period motion is eliminated. This 
may be accomplished by assuming a constraint such 
that « = 0 and w = Uy, sin @. Taking, for simplicity, 
Zp = J = 0, Eqs. (9c) and (9d) become 
(m,r?)g = (w/Uo)M, + (qce/2Uo)M, — W,x, (11a) 
—m,X, = W,0 + kx, (11b) 
In these equations —M, and — WM, are the aerody- 
namic spring and damping constants. Using the fuel 
and airplane short-period motion natural frequencies 


wp? = k/m,, w, = —M,/m,r’ 
and taking 
D, = —cM,/2Uom,r’ 
Eqs. (11) may be written 
8+ Dw + w28 + (W,/W, (g/r2)x,) = 0 (12a) 
Xp + wy*X, + gd = 0 (12b) 


To account for the possibility of viscous damping of the 
fuel motion by baffles, say, Eq. (12b) may be modi- 
fied by the introduction of a viscous damping term, 
D,Xp, 


&, + D,t) + w,2x, + gd = 0 (12c) 


The characteristic equation for the set of simul- 
taneous differential equations, Eqs. (12a) and (12c), is 


(2 + Dir + w,”) (A? + DprA + @,?) 


(W,/W,) (g/r)? = 0 (13) 


When the ratio of the effective fuel weight to the total 
weight approaches zero, the meaning of this equation is 
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simply that the motion that exists is the aircraft short- 
period oscillation and an uncoupled oscillation deter- 
mined by the natural period and damping of the fuel. 
For finite values of the parameter (W,/W,) (g/r)? = », 
the frequency and damping of each of these modes are 
modified. To determine the extent of this modifica- 
tion it is necessary to extract the roots of the quartic. 

For the special case of zero aerodynamic and fuel 
damping, extraction of the roots is simple, 


Mt + (w? + w,”)r\? + (w?w,? — vy) = O- (14) 
P (w,? + wy”) cat — en? - 
= + wn t ) ! 
r 1 3 = = ( ‘ u ) +y (15) 
and, for w, = w; = a, 
i = ae — ot & 4/5 (16) 


Only the frequencies of the oscillations are affected. 

While for specific cases Eq. (13) may be factored by 
numerical methods, it is desirable to have at least ap- 
proximate relations from which the change in damping 
and frequency from the uncoupled modes may be ob- 
tained. This may be done by assuming a power series 
expansion in orders of v for the coefficients of the bi- 
quadratic equal to the original quartic and then evalu- 
ating the coefficients of the power series expansions. 
Another way is to assume a perturbation in the auxiliary 
function, \ = A* + yu, where X* satisfies the original 
quartic with vy = 0, and to evaluate wu in terms of the 
coefficients of the original quartic and \*. The first 
method results in more convenient relations and will be 
used below. 

By power series development of the coefficients of 


the biquadratic, 


(A? + ad + b) (A? + ad + B) = 
(\ + od + bo) (A? + aX Tr Bo) —v=0 (17) 


where 
1 
a=atayvtayv+...6= 2. a;v 
1 0 
5b = Sob, 
0 


t 
a= Ya;v' 
0 
1 
dBi’ 
0 


Carrying the power series development only to terms of 


B = 


first order in p, 


aq = 
(do = Qo) 

Bo(do a ao)? — Gite — ay) (bo — Bo) + (bo = Bo)? 
(18a) 

a, = —-ay (1Sb) 

Q,| ao(do = @) — (bo — Bo) ] 
A= -— (1Sc) 
(do — a) 
by = A;(ayo — Qo) — Bo (18d) 
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This approximate factorization will fail at the point 
where two real repeated roots occur. Here, one of the 
real repeated roots may be obtained from the rea] root 
of the cubic resulting from differentiation of the original 
quartic. Comparison of Eqs. (13) and (17) and ap- 
plication of Eq. (15) lead to 


AD, = 
(D, — D,)v 
w,?(D, — D,)? — D,(D1 — Dy) (a2? — wy”) + (w2 — w,2)? 
(19a) 
AD, = —AD, (19b) 


P , D,\? AD, : 
2,7 = wr" — > —_— D, oe D, (w," — W,*) (19¢ 
a oe 2 y Pes ‘ 
tp = Wp 2 D, — D, Wy Wy”) 


where 2,” are the squares of the frequencies of the new 
motions to first order in v. [w,;? — (D;/2)?] 
squares of the frequencies of the uncoupled modes, and 
w,” are the squares of the natural frequencies of the short 


(19d) 


are the 


period and the fuel motion. For small v, the approxi- 
mate roots found from Eq. (19) agree well with those 
obtained by numerical methods. 

The new motion is an oscillation in which the damp- 
ing and the square of the frequency vary linearly with 
the coupling quantity v. For typical values of », the 
coupling between the new mode and the aircraft short- 
period oscillation is small. This is illustrated by Fig. 6 
for a specific case of an aircraft with a short-period nat 
ural frequency of 1 cycle per sec. and damping such 
that the amplitude is diminished to one-half in two 
cycles. Zero viscous fuel damping was assumed, and 
a typical value of (W,/W,) (g/r)? = v = 1.0 was se- 
lected. It is seen that even at resonance the damping 
of the short period mode is only slightly diminished 
(6 per cent). The frequency was found to be altered 
by only a fraction of a percent. 
result for the fuel motion.| At a certain small value 


Fig. 7 shows a similar 


+ The actual motion of the airplane and of the fuel are each 


the sum of two motions. The words “short period” and “fuel 


motion” are merely used to identify the modes associated with 


certain roots of the characteristic equation. 





2, FUEL NATURAL FREQUENCY 
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Fic. 6. Per cent change in airplane short-period damping 
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Ol W/L 


a p 
be positive. 


AND 


it is possible to obtain a pair of real roots and 
air of complex roots, and one of the real roots may 
Physically, this situation corresponds to 


the introduction of a weak spring in the equivalent 


fuel mass-spring system. 


The restoring force of this 


spring 1s insufficient to overcome the weight force com- 
ponent acting on the fuel upon rotation of the airplane 
and the excursion of the fuel mass increases without 


limit. 
into a subsidence and divergence. 


Here, the fuel motion oscillation degenerates 
This case arises 


when the fuel tank is nearly empty (the weak spring 


implies a long natural period or a long pendulum). 


Even though no viscous damping of the fuel motion 
was assumed, the coupling introduced a small damping, 


as shown by Fig. 7. 


For practical purposes this damp- 


ing is inadequate, however, and diminution of the fuel 
motion must be secured by some physical apparatus 


such as baffles. 


The influence of a baffle that will pro- 


vide a damping such that the fuel motion amplitude is 
reduced to one-quarter in one cycle has been calculated. 
It was found that with fuel frequencies close to that of 
the airplane short period, the presence of the baffles 
improves the airplane damping by about 10 per cent. 


Returning to Eq. (9), it is recalled that with motion- 
less fuel the phugoid motion of the aircraft may be ob- 
tained, approximately, by considering flight at con- 
stant angle of attack and with small inertia of the air- 


craft. For this case, from Eq. (9), 


mu = uX, — Wad 
—m,Uid = uZ,, 


and the characteristic equation is 


. ae W Zu 
it~ Sy Sg 
mM, m,?Uo 


By substituting the usual relations 


: 2 #1 : 2D 
X ,= ; (; pl 5) Cp = n 


om.  N oW, 
Uo (5 elas Cr ™ U5 


into Eq. (21), there results 


— g | D ss ay 
2. a ee & 


(20a) 
(20b) 


(22?) 


which gives the customary approximate relation for the 
damping and frequency of the long-period oscillation. 


With sloshing fuel, the angle of attack in the phu- 


goid is no longer constant. 
and g and g small, 


(w/Us)Ma — Wpx, = 0 


Qa 


while from Eq. (9d), with #, and # small, 


W,0 + kx, = 0 


From Eq. (9c), for z, = 0, 


(23a) 


(23b) 


With small m,/m,, Eqs. (9a) and (9b) may be written 
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Fic. 7. Damping of the fuel mode. Airplane natural fre- 
quency, w% = 1 cycle per sec. Airplane damping, 2 cycles to '/2 
amplitude. (W,/W:)(g/r)? = 1.0. 

Xx Xe 
“=u w — gv 
mM, mM, 
or 
; ' # Aw w,’M, 
“=u w a (24a) 
mM, mM, UW, 
since by Eqs. (23a) and (23b) 
—k —w,’M, 
v= — x, = =e OY 
W, gUW, 
and 
w,?M, : : 
(« + — w) = uZ, + wZy (24b) 
gw p 


The characteristic equation for the set, Eqs. (24a) and 


(24b), is 


v-Pee(2 as 
mM, m,\1 + (w,?M./gW,) 
| Lw ( ) Zu 
; : = Xx 
m,m,\1 + (w,?M,/gW,) mM, 


( l \(= ns) | ~ 
; . = Q (25) 
L + (w,?M,/gW,)/\m, UW, 


and, with 


wy?M ,/gWy = wp?w,?/v 

Ad = —2D/Uo 

Ze = —2W,/U 

K = 1 + (w,’w,"/r) 
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: = (W, U) = [(OD Oa) Uo] 
Le = —(D, U5) [(OL Oa) Uo] 


gf [DQK+ : bus a] 
i= se 
U, Res +o |* 


D (2K + ) 
i 4/2 — 
Ei aK )* 


For typical values of airplane short-period natural fre- 
quency, w,, fuel natural frequency, w,, and the ratio, v, 
the difference between the dampings and frequencies 
given by Eqs. (26) and (22) is negligible. 

In the rather unlikely case that the fuel natural fre- 
quencies and long-period frequencies are nearly in 
resonance, the fuel and phugoid motions must be con- 
sidered simultaneously. A quartic characteristic equa- 
tion will be obtained which must then be analyzed by 
methods similar to those already described. 

Thus it is seen that the longitudinal motion of an air- 
craft with sloshing fuel is composed of three essentially 
the phugoid, the airplane 
The coupling 


independent oscillations: 
short period, and the fuel oscillation. 
between these motions is, for practical purposes, neg- 
ligible. Satisfactory damping of the fuel oscillation 
must be obtained by the use of a damping device, such 
as baffles, within the fuel tank. These results are in 
agreement with less approximate analyses in which fac- 
torization of the characteristic sextic was performed by 
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Cp = Cp, + (C,? ‘wA.R.) 


it is found that 


C, a C 1 


+5 | 7) ( 4 “1s 4 (1 _ 2Cra )] ie 
2K K iS ce. rA.R./ |f 


numerical methods. All the variables in this problem, 
u, w, 3, and x,, will have solutions of the same mathe. 
matical form. Flight experience has shown that the 
new oscillatory mode is apparent to an observer within 
the airplane as a pulsing in the acceleration 7. This 
may, perhaps, be more readily comprehended by visu- 
alizing the motion of a cart containing a mass attached 
to it by a spring, as shown in Fig. 8. Both the cart 
and the mass oscillate longitudinally about the system 
center of gravity while that point is transported (with 
phugoid and short-period motions) in space. 


BUFFETING 


A General Description of the Aircraft Buffeting Problem 


Buffeting is a forced vibration in which the external 
forces are generated by eddying flows.¢ Fluctuating 
forces may act on the component producing the eddy- 
ing flow, or they may be induced on some other part 
of the aircraft which is in or near the eddies. The 
eddying flow can be produced by wing or tail shock- 
stall at high subsonic speeds; by separation on the 
wing at the low speed stall; by separation from dive 
flaps, bomb-bay doors, landing-gear doors, wing flaps, 
canopies, radomes and guns; by the mixing region 
created by a jet or by a propeller slipstream; and by 
numerous other things. The most important cases of 
low-speed buffeting resulted from separation at the 
wing-fuselage juncture of some low-wing monoplanes. 
At higher speeds the principal buffet problem 1s con- 
cerned with transonic wing and/or tail shock-stall. 

At the present time, the problem of buffeting has not 
been solved. No design methods are available whereby 
stall warning or buffet loads can be predicted. It 1s 
not possible to estimate the effect of wing or tail geom- 
etry, tail position, Mach Number, and other variables 
on buffeting. No completely satisfactory methods for 
investigating buffeting during routine wind-tunnel tests 
have been developed. Aircraft designers have only 
experience on generally similar designs (where such 
are available) for guidance, and the buffeting char- 

{ Flutter is a self-excited vibration not generated by external 


forces. Buffeting is a response problem; flutter, a stability 


problem. 
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acteristics therefore remain unknown until the flight- 
test period. This situation has made it necessary for 
some aircraft to undergo long development programs 
in flight and even then to be limited in performance 
by placards. What is needed is: (1) a procedure for 
designing aircraft with the desired degree of stall 
warning and for predicting buffet loads; (2) a procedure 
for checking those calculations in the wind tunnel; 
and (3) a procedure for investigating buffeting in flight 
and for correcting such troubles as may arise. 

An analytical solution to a problem requires as its 
starting point a suitable description of the phenome- 
non under investigation. For vibration problems it is 
necessary to be able to specify inertias, dampings, re- 
storing forces, (or their analogies) and external forces. 
It is necessary to understand the kinematics of the prob- 
lem and to visualize the constraints and couplings re- 
lating the variables. 

The masses (inertias) of the buffeting problem are 
those of the aircraft itself. The dampings are both 
structural and aerodynamic in origin, as are the restor- 
ing forces. The external forces result from the fluctu- 
ating loads acting on the stalled body or from the in- 
duced loads created by its wake. It may be of interest 
to discuss how some of these forces may enter into the 
problem. First, suppose a rigid body like a wing is 
supported in a steady stream. The flow configuration 
will depend on the Mach Number, J/, the Reynolds 
Number, R, and the attitude of the body relative to 
the stream, a. The magnitude of the forces will de- 
pend, in addition, on the static pressure, p. These 
variables determine the external force that is the input. 
Since the body is rigid and stationary, there is no damp- 
ing and no inertia. The restoring force supplied by the 
support is then equal and opposite to the aerodynamic 
force. The essentials of this description are shown in 
Fig. 9. The outputs of this system are the force on 
the support and the flow variables in the wake. For 
certain values of the input parameters, the wake be- 
comes large, both the wake and the force fluctuate, and 
buffeting occurs. For the nonrigid case, the displace- 
ments, velocities, and accelerations must also be con- 
sidered. 

Before discussing more complicated cases, it is noted 
that it is not presently possible to estimate the fluctu- 
ations of force and wake in the buffet régime for the 
simple case just discussed, even for elementary con- 
figurations like a two-dimensional cylinder. 

A more realistic case is one in which the wing is at- 
tached to a fuselage and the wake passes near a tail. 
The output we should like to calculate (for example) 
is the vertical acceleration at the pilot’s compartment. 
The relations between the inputs and outputs acting on 
the aircraft components are shown in Fig. 10. It is 
these things that must be calculated in a buffeting 
theory. Items (1), (2), and (3) (numbers correspond 
to those in Fig. 10) have already been discussed. Item 
(4) is the load in the fuselage-tail attachment. Item 
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Fic. 9. Schematic diagram for an elementary case of buffeting. 
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Fic. 10. Possible schematic diagram for buffeting. 











(5) is the motion of the tail resulting from fuselage 
motion, while item (6) is the motion of the wing result- 
ing from fuselage motion. Whether this representa- 
tion of buffeting due to wing stalling is incomplete or 
overly complicated could be investigated in a wind 
tunnel (with suitable attention to the constraint and 
dynamic effects of the model support) using a flutter- 
type model. Tests of a complete configuration com- 
pared with those with the tail off would indicate whether 
the wing or the tail is the principal contributor to the 
fuselage vibrations. Tests with the wing in place to 
provide a wake, but supported separately, would also 
help answer that question. 

The matrix of differential equations that are the 
analog of the schematic diagrams just discussed can- 
not be completely written with presently available 
knowledge. To do so might require an analytical 
representation of some of the following items: 

(1) The fluctuations of force and moment on rigid 
stalled bodies as a function of body geometry, a, J/, 
R, and p. 

(2) The aerodynamic damping of an oscillating 
stalled wing. 

(3) The vertical velocity fluctuations in the wake 
of a stalled body as a function of location relative to 
the wake; body geometry; body a, M/, R, and p; 
and the body’s oscillations. 

(4) The response of a tail to the disturbances in the 
wake. 

(5) The effect of oscillations on the wake of a stalled 
wing. These questions relate mostly to the aerodynamic 
data of buffeting. The use of such data in the theory 
of the vibrations of aircraft will provide a theory of 
buffeting. The last section of this paper discusses sev- 
eral attempts to calculate some of these aerodynamic 


variations. 
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Words and Phrases 


von Karman and Biot,® under the title ‘‘Words and 
Phrases,’’ have a collection of definitions of concepts 
and operations in common use by physicists and engi- 
neers. They also note the existence of a collection of 
legal terms with the same title which includes espe- 
cially the legal interpretation of commonly used ex- 
pressions. A number of words, phrases, concepts, 
and operations have been associated with buffeting. 
Eberle’ has made an interesting survey of some of 
these concepts. An informal review of these is in order 
here to see if they convey any precise information con- 
cerning this phenomenon; also, one new phrase will be 
introduced. 


(1) The Vortex Trai! and the Strouhal Number. 
To obtain the frequency properties, a large part of the 
buffeting literature assumes that the wake of a stalled 
wing is a Karman vortex street (or trail). A single 
characteristic frequency, given nondimensionally by 
the Strouhal Number, carries all the energy of the trail. 
Severe buffeting is supposed to be avoided by keeping 
the empennage natural frequencies far from that of the 
trail, thus avoiding resonance. These ideas are mainly 
based on data from model tests of cylinders and airfoils 
at low Reynolds Numbers and probably cannot be 
carried over to full scale. This was pointed out by 
von Karman!! 
flows, the vortex street pattern only appears for Reyn- 
olds Numbers below a definite limit of the order of 
10°. Above this value, the regular character of the 
wake disappears and is replaced by a completely turbu- 
lent one.’ The upper limit to the Reynolds Number 
at which the trail exists is also discussed by Den Har- 
tog.'” It would seem that a more realistic concept is 


himself: ‘‘Even for two-dimensional 


that the full-scale wake has a continuous frequency spec- 
trum. This implies that resonance with the structure 
cannot be avoided. The solution to buffeting must 
therefore depend largely on reducing the magnitude of 
the fluctuations arising from flow separation and/or 
turbulent wake from the struc- 


keeping the away 


ture. 


(2) The Total Head Wake 


head loss in the wake of a stalled wing are often made 


Measuremeuts of total 
in buffeting investigations. A considerable amount of 
evidence has accumulated which shows that such meas- 
urements are of limited value. Duncan measured a 
““buffeting’’ wake with a detector, which was essentially 
a flexible tail, and found that heavy buffeting can occur 
even outside the total head wake. Eberle mentions a 
case where, in flight, buffeting occurred with the tail 
outside the wing wake. In tests of a low-wing mono- 
plane, White and Hood’ used a fillet to reduce the tail 
buffeting vibrations to one-seventh their original mag- 
nitude. This, however, was accomplished without 
moving the tail from the center of the total head wake. 
In view of these facts, it appears that something more 
meaningful than the total head wake should be found. 


AERONAUTICAL 


SCIENCES—APRIL, 1952 


It is suggested that the “‘vertical velocity fluctuations 
wake may serve that purpose. 

(3) Experiments.—Wind-tunnel tests on buffeting 
have often been made at low Reynolds Numbers, Dun- 
can’s tests, for example, were made at a Reynolds 
Number in the order of 10°. 
to discover the effect of Reynolds Number on wake and 


It would seem worth while 


force fluctuations and then to test only at sufficiently 
Low-speed buffeting might best be stud- 
Other 
questions of testing technique might well be discussed 


large scales. 
ied on actual aircraft in full-scale wind tunnels. 


here. It has already been suggested that vertical ve- 
locity fluctuations may be the important quantity to 
measure in the wake. This may be more useful than 
total head loss measurements or fluctuating wake pres- 
sure measurements. Finally, it is suggested that opti- 
cal methods, such as motion pictures of smoke flow, 
may be useful but must be accompanied by measure- 
ments of fluctuations of force and moment on the stalled 
body and by measurements of fluctuating vertical ve- 
locities (angles of attack) in the wake to provide mag- 
nitudes of the forcing functions. 
(4) Dynamics of Buffeting. 
the vibration characteristics of the airplane into the 


Attempts to bring 


analysis of buffeting have generally considered only the 
This 
is obviously important if the buffet loads on the tail 
Estimates of the magnitude of the 


dynamics of the tail (reference 14, for example). 


are to be studied. 
loads produced by the tail vibrations, however, show 
that they are small compared to the weight of the air- 
craft and cannot be the source of the large c.g. accel- 
erations often measured in buffeting. It would seem in- 
stead that these accelerations must come from the flue- 
tuations of wing lift or from the large amplitude motions 
of fuselage resonance. The wind-tunnel tests suggested 
previously in connection with the schematic diagram of 
buffeting would aid in choosing between the two possi- 
bilities, although it is felt that the fuselage resonance 
possibility is more likely. If this supposition is correct, 
it will be necessary to include the fuselage dynamics in 
analyses of buffeting. 

(5) Stochastic Process.—A 
one in which randomness plays a large part. 


stochastic process 1s 
Games 
In physics, Brownian motion 


Flight through 


of chance are stochastic. 
and molecular diffusion are examples. 
a turbulent atmosphere or taxiing over rough ground 
are examples in aeronautics. Generally, only certain 
mean values of the important quantities enter the anal- 
ysis of stochastic processes, while the concept of phase 
is completely lost. For example, in coin tossing, it 1s 
possible to calculate whether or not a head or tail will 
turn up only on the average and never for a specific 
Liepmann’ has suggested that the flow in the 
wake that methods are 
therefore appropriate for the analysis of buffeting. 

Thus, Strouhal Number, 
head wake, and dynamics of only certain portions of al 


toss. 


is stochastic and statistical 


while vortex trail, total 


aircraft structure may be meaningful concepts in them- 
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selves, they are of dubious significance as related to the 
buffeting problem. It is believed that buffeting is de- 

ndent fundamentally upon fluctuations of a random 
nature and that consequently analytical and experi- 
mental studies must recognize this by application of 
appropriate mathematics and the measurement of more 
directly relevant quantities. 


some Calculations on Buffeting 


It might be helpful to have a way to estimate the 
decay of the disturbances as a tail is raised above a 
wake. This was studied by Graham® by considering 
the Karman vortex trail to represent the wake be- 
hind a two-dimensional wing on which flow separation 
had occurred.t| While it was realized that the vortex 
trail concept did not lead to a correct prediction of the 
frequency characteristics of the wake, it was hoped 
that it might serve here as a rough approximation to 
help study the spread of the fluctuations. It was 
found that, at some distance from the trail, the maxi- 
mum induced vertical velocities occurred on a vertical 
line midway between two adjacent alternate vortices. 
On that line, the amplitude of the velocity was given 
by 
V, = (1'/2b) { sech [(2ry/b) — 0.284] + 

sech [(2ry/b) + 0.28r]} (27 


The distribution of the vertical velocity amplitude given 
by this relation is shown on Fig. 11. If the dimension 
b, the distance between vortex rows, is taken to be the 
width of the total head wake, then this fluctuation dis- 
tribution curve agrees moderately well with the low 
Reynolds Number ‘“‘detector’’ data of Duncan. It is 
sen that a large disturbance can still exist a half total 
head wake width from the edge of the wake. It was 
noticed by Graham that the lateral decay of the ver- 
tical velocity in the flow field of a sine-wave-shaped wall 
approached that given by Eq. (27) for lateral positions 
not too close to the wake. This suggested the use of a 
wave-shaped wall in simulating the field outside the 


tEberle'® suggested the same approach. 
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Lift of a rigid thin airfoil in an incompressible isotropi 
cally turbulent flow. 
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wake. Since the flow about the wave-shaped wall is 
known as a function of Mach Number, a compressibility 
effect is thereby indicated. Graham showed that the 
effect of compressibility was to spread the disturbances 
even farther laterally. The effect was small, since 
the Mach Number involved is that of the vortex trail, 
which is considerably lower than the airplane Mach 
Number. 

Some new ideas on buffeting were recently reported 
by Liepmann.’ He suggests that some well-established 
concepts of turbulence and statistics may be usefully 
carried over to this problem. The wake is considered 
to be turbulent and is described by a continuous power 
The energy of the wake’s fluctuations is dis- 
tributed over a wide frequency range. Hence, reso- 
nance with the tail and fuselage is unavoidable. Flight 
measurements taken during buffeting show that the 
structure does, in fact, vibrate at its natural frequencies. 
Since the fluctuations in the wake are random, buffeting 
The central points in Liep- 


spectrum. 


is a stochastic process. 
mann’s suggestions are that (1) the wake of a stalled 
body is turbulent and is described mathematically by’ 
mean values and by power spectra, and (2) the buffeting 
vibrations of an airplane are a stochastic process. To 
illustrate some of the concepts and methods involved in 
this description of buffeting, Liepmann calculates the 
lift of a rigid two-dimensional thin airfoil in a wake 
which is approximated by isotropic turbulence. Since 
only statistical data are known about the wake, only 
the mean square lift can be calculated. No concept of 
phase even enters the calculation, since the random- 
ness of the fluctuations in the wake means no informa- 
tion on phase is available. The mean square lift co- 
efficient is found to be 


CL? = 42?(w?/U*)fen(mc/L) (28) 


Thus, the mean square response of the airfoil of chord « 
in an isotropically turbulent stream of velocity UL’ de- 
pends on the mean square value of the vertical velocity 
fluctuations, w*®, and the scale of the turbulence, L. 
This relation is plotted in Fig. 12. It is seen that the 
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response to large eddies (Z large) approaches the 
steady-state value given by the lift curve slope 27, while 
the response to small eddies approaches zero. Eq. 
(28) also furnishes useful information concerning ex- 
periments. 

It is necessary that the spectral distribution of 
the fluctuations and the relation between the scale 
of turbulence and the quantity that controls it 
be established. This calculation is a simple example 
of the procedure needed to calculate response from a 
fluctuating wake input. It represents, on Fig. 10, 
wake (3) going into the tail geometry and producing a 
load. The motion of the tail, which brings in the in- 
ertia, damping, and stiffness forces, also has to be con- 
sidered. This consideration will show that the ampli- 
tude of the motion will be controlled by the damping, 
since the vibrations will be resonant. 


CONCLUSIONS 


Analytical approaches to two aircraft dynamics prob- 
lems have been discussed. A scheme for representing 
oscillations of fuel in a tank was developed and used in 
studying the longitudinal motion of an aircraft with 
one partially full fuel tank. The motion of the air- 
craft was found to be composed of three loosely coupled 
oscillations: the phugoid, the airplane short period, and 
the fuel natural mode. The fuel oscillation has a fre- 
quency essentially the same as that of the fuel in a sta- 
tionary tank, and significant damping was obtained only 
by the use of tank baffles. 
for a rational approach to buffeting on the basis of an 


Arguments were advanced 
aeroelastic response problem. It was suggested that 
ideas from turbulence and the mechanics of stochastic 
processes may be useful in setting up a_ buffeting 
theory. 
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The Compressible Boundary Layer 


* 


GEORGE B. W. YOUNG? ann EARL JANSSENt{ 
The Rand Corporation and University of California at Los Angeles 


SUMMARY 


The laminar boundary-layer equations in integral form as ob- 
tained by Moore! are solved on the mechanical differential an- 
alyzer. The method of analysis permits the use of the best 
known variation of the air properties with temperature. Solu- 
tions are obtained for Mach Numbers ranging from incompres- 
sible speeds to the hypersonic régime at ambient temperatures 
of 100°, 400°, and 800°R. for both insulated-plate and heat- 
transfer cases. From these results, an approximate design pro- 
cedure for the calculation of the characteristics of both the lami- 
nar and turbulent boundary layers is obtained and presented as 


an aid to missile designers. 


INTRODUCTION 


Ww" THE EXCEPTION OF MOorRE’S SOLUTION,! 
all existing solutions?~* of the laminar boundary- 
layer equations employe either a power law for the 
viscosity variation with temperature or a constant 
Prandtl Number or both. These assumptions limit 
their usages to the range of Mach Number and ambient 
temperature to which the assumptions are valid. Al- 
though these restrictions are generally not severe (some 
solutions being valid to JJ ~ 5), the present analysis 
removes these limitations and shows the trends of the 
characteristics of the boundary layer at high Mach 
Number and low ambient temperature. The thermo- 
dynamic and physical properties of air, pertinent to the 
boundary-layer solution, employed herein are given in 


the Appendix. 


The laminar boundary-layer equations in integral 
form as obtained by Moore are solved on the mechanical 
differential analyzer. The variation of the flow prop- 
erties through the boundary layer are obtained for 
ambient temperatures of 100°, 400°, and SOO°R. by a 
converging procedure at several wall shear stress con- 
ditions to yield Mach Number variation and at several 
heat-transfer conditions including the insulated plate. 
The ambient temperatures of 100°, 400°, and SOO°R. 
ate chosen in anticipation that the present results will 
be helpful to investigators attempting to correlate re- 
sults from hypersonic tunnels (where the temperatures 
are low) with flight-test data (where the ambient tem- 
peratures are high). 
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EQUATIONS OF THE LAMINAR BOUNDARY LAYER AND 
THEIR DIFFERENTIAL ANALYZER SOLUTIONS 


The laminar boundary layer over a flat plate in two- 
dimensional compressible flow, with steady-state and 
zero pressure gradient conditions, is described by three 


equations: 


Conservation of Momentum 


Ou 4 Ou ra) ( a 
u v = 
eax | Oy dy ay 


Conservation of Mass 


[O(pu)/Ox] + [0(pv)/Oy] = 0 


or ou\?* 
k + pu 
oy Oy 
The following steps and transformations have been 


made by Hantzsche and Wendt and Moore! and are 


simply summarized here: 


Conservation of Energy 
oh oh o 


u + pv = 
F Ox r oy oy 


(a) Transform the equations to total differentia, 
equations, with z as the single independent variablel 


by letting 
z= y/V xX; 02/O0x = 
O2/Oy = 1/V x; v=vVx 
(b) Transform to u as the independent variable by 
introducing G(u) = uw (du/dz) and reduce the system to 
two equations 


dG dh d ( k oI : 
= ( + G 


= 7 
du du du \yc, du 
G(d°?G/du?) = —(1/2)puu 


by combining the momentum and energy equations 
and the momentum and mass equations. 
(c) Put into dimensionless form by letting 


u G 
. ie / 
V (1/2)ppun(he) 


Pr =-pc,/k 


/ 


V hr 
Z — h hp; 


where the subscript B refers to a base temperature 
(taken as 400°R. for this investigation). 

(d) Transform to a new independent variable, /, 
to avoid a singularity which appears in the equations 
rearranged for differential analyzer solution by letting 


9299 
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Fic. 1. Schematic diagram for the differential analyzer solution 
of the laminar boundary-layer equations. 


dn dw pu 


é 7 dt . dt PBHB 7 


The two resulting equations are 


d*j - P | dj 4 (“’) d(1/Pr) 4 ] 
dt? je dt dt dt g 


dg/dn = —w 


since w = 0 when y or fis zero. The equations, along 
with certain associated transformation equations, are 
solved by integrating through the boundary layer from 
the plate (¢ = 0) to the free stream (¢ = ¢;)._ The solu- 
tion consists of j, dj/dt, n, g, and J u/us dt as functions 
of f. 

The complete set of five equations rearranged for 


analyzer solution are 


dj _ 
dt 
“p i| 4 ij + “a(;,) + ‘a| 
- r w y 
| ’ / " le Pr fm 
g=— fwd 
n= JS gdt 
w So a f n dt 
PBHB 
4 IS” (dj/dt) dt 


The sali that relates the original independent 
variables x and y to the independent variable ¢ is 


Vh y 2 
PB B- = = dt 
2uzr Vx LB 


The reciprocal of 1/Pr is obtained by two integra- 
tions, 
log (1/Pr) = f° Pr d(1/Pr) 
Pr — Sf Pr d{log (1/Pr)] 


I 


The schematic diagram for the machine solution is 
shown in Fig. 1. The air properties (1/Pr, pu/paus, 
and u/s), as given in the Appendix are supplied from 
input tables during the course of the solution. Fur- 
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ther details concerning the differential analyzer are 
given in references 6, 7, and 8. 

This is a two-point boundary problem requiring six 
boundary values to determine a solution uniquely, 
The four values at the plate are: 7 = ny, = 0; w = 
Wy» = 0; 7 = Jw, assigned for desired plate temperature: 
and (dj/dt) = (dj/dt), = 0, for the insulated plate 
case; Or g = gy, assigned arbitrary values for the heat- 
transfer case to give a range of Mach Number. 

The two values in the free stream are: 7 = j;, deter- 
mined by the free-stream temperature; and g = g, = 0), 
Note that ¢; itself is not known but is determined when 
g goes to zero. 

The solution* of this two-point boundary problem on 
the differential analyzer involves guessing a value for 
Zw» {or for (dj/dt), in the heat-transfer case], setting it 
in the machine (along with the known values at ¢ = (0), 
If the value of j 





and running the machine to g = 0. 
as g approaches zero is the desired one, then the solu- 
tion is considered complete. If not, a new value for 
£ or (dj/dt), must be guessed and the procedure re- 
This converging procedure is aided by plot- 
(dj/dt), for the heat-transfer 


peated. 
ting 75 VS. £ [or js VS. 
case]. The rapidity of this convergency depends on 
the approximation of the first guess, the linearity of j; 
with g,, or (dj/dt),, and the sensitivity of the machine 
in the particular region of the parameters. The param- 
eter combinations employed are tabulated in Table 1. 


* In certain regions of the parameters the value of 75 obtained 
from the machine was subject to some dispersion—i.e., repeat- 
ing runs with the same initial values gave different values for j;. 
This happened at large values of jy and small values of 75, for 
which there were corresponding large values of gy. The inte- 
grator generating Sg dy [contributing to one of the terms in 
(dj/dt)| was observed to be ‘“‘jumpy.’’ When the machine was 
“‘jumpiness’”’ increased. Correspondingly, the 
By running the machine slowly, the 


run faster, the 
dispersion in 75 increased. 
dispersion decreased to zero. The trouble was believed to be 
due to the slipping of the integrator wheel of this one integrator 
when the acceleration of its moving parts reached a certain 
threshold value. By keeping the speed down, this threshold 
value was never reached. Noruns that could not be repeated or 
which did not give values of 75 and gw or (dj/dt)w which lay on a 


smooth curve were considered valid. 


TABLE 1 
Parameter Combinations 


—Insulated Plate - Heat Transfer 


Assigned Desired Assigned Assigned Desired 
Ju ji Jw gu j6 
l 0.25, } l 2.06 0.25 
2 0.25, 1,2 I 3.27 0.25 
3 0.25, 1,2 1 5.82 0.25 
{ 0.25, 1,2 l 9.94 0.25 

6 0.25, 1,2 3 2.06 0.25, 1,2 

10 0.25, 1,2 3 3.27 0.25, 1,2 

15 0.26, 1,2 3 5.82 0.25, 1,2 

25 3 9.94 0.25, 1,2 

40 Lo 3 12.95 1,2 

6 2.06 1,2 

6 3.27 0.25, 1,2 

6 5.82 0.25, 1,2 

6 9.94 0.25, 1,2 

6 12.95 1,2 
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The Characteristics of the Laminar Boundary Layer at the Nominal Ambient Temperature of 100 


COMPRESSIBLE 


BOUNDARY 


TABLE 2 


to 


LAYER 


R. (Insulated Plate) 


15 z sing 6 6* , g*° 

M (°R.) (CV Re)s (CV Re) ; (°R.) (°R.) x Vv Res x V Res < V Res 

23 99 0.690 0.677 0.851 400 380 13.49 9.86 0.68 
6.42 101 0.651 0.693 0.841 800 642 22.87 19.19 0.64 
8.37 99 0.620 0.702 0.820 1,182 911 31.34 28.20 0.63 
9.57 100 0.593 0.709 0.791 1,550 1,200 37.22 34.24 0.58 
11.67 104 0.555 0.714 0.756 2,246 1,680 19 61 $5.30 0.55 
15.45 102 0.501 0.721 0.711 3,565 2,570 

90.45 90 0.462 0.721 0.670 5,135 3,820 97.09 94.59 0.46 
TABLE 3 
The Characteristics of the Laminar Boundary Layer at the Nominal Ambient Temperature of 400°R. (Insulated Plate) 
Ts T, ag 6 6* - g** - 

VU (*R.) (CV Res (CV Re)u r CR) (He) x V Res x V Res eV Res 
0.22 400 0.674 0.674 400 400 5.02 1.69 0.64 
2.43 102 0.629 0.681 0.842 800 738 7.51 3.94 0.64 
3.47 400 0.601 0.692 0.812 1,182 1,000 10.05 6.10 0.60 
4.28 101 0.575 0.699 0.783 1,550 1,313 11.08 7.83 0.56 
5.52 405 0.541 0.705 0.745 2,246 1,782 14.26 11.15 0.53 
7.49 401 0.499 0.714 0.703 3,565 2,700 19.96 16.81 0.48 
9.59 388 0.458 0.716 0.664 5,135 3,970 26.71 24.07 0.45 
12.44 398 0.417 0.72 0.630 8,160 6,000 36.11 34.08 0.40 
15.95 396 0.383 0.734 0.608 12,645 8,880 50.92 18.81 0.38 

TABLE 4 
The Characteristics of the Laminar Boundary Layer at the Nominal Ambient Temperature of 800°R. (Insulated Plate 
75 7 ig 6 5* g** 

VW (R.) (GWVRe)s (CV Re)u r (R.) (RR) x V Ree V Res ¢ V Reb 
0.26 794 0.671 0.672 800 800 $94 1.67 0.64 
1.74 800 0.638 0.680 0.792 1,182 1,120 6.14 80 0.638 
2.45 808 0.612 0.687 0.767 1,550 1,416 7.3 3.77 0.61 
3.52 797 0.576 0.696 0.733 2,246 1,900 9 02 ».62 0.57 
$.95 810 0.533 0.710 0.693 3,565 2,736 11.88 8.55 0.52 
6.45 792 0.490 0.712 0.659 5,135 $040 15.54 12.24 0.49 
8.54 805 0.448 0.723 0.626 8,160 6,000 20.30 17.82 0.48 
11.33 765 0.405 0.729 0.605 12,645 8,950 28.40 26.57 0.39 

TABLE 5 
The Characteristics of the Laminar Boundary Layer at the Nominal Ambient Temperature of 100°R. (Heat Transfer Plate 
Ts ‘ Nas 6 = s*  ,- s** 

M ("R.) (°R.) (C/V Re)s (GV Redu (°R.) av PrsV Res x © Res x Res x» Res 

8.63 97 400 0.667 0.655 530 0.312 19.62 16.14 0.65 

11.84 101 100 0.639 0.629 720 0.318 26.79 23.35 0.63 
19.19 94 $00 0.584 0.572 1 350 0.331 $15.59 42.61 0.56 
29.64 95 100 0.514 0.504 2,500 0.323 7.28 76.06 0.49 
8.93 101 1,182 0.615 0.696 957 0.307 32.28 28.74 0.62 
12.37 102 1,182 0.594 0.673 1,195 0.329 38.15 34.43 0.59 
18.74 105 1,182 0.553 0.629 1,695 0.339 51.46 18.20 0.53 
32.20 85 1,182 0. 500 0.561 2,925 0.315 104.9 102.9 0.48 
13.39 96 ? , 246 0.553 0.708 1,750 0.313 55.99 52.97 0.54 
22.40 80 0.524 0.668 2,313 0.332 

31.88 Q? 0.468 0.605 $+, 100 0.312 


19d 


Run 

No 
10¢ 
lle 
l6e* 
l7g 

(rep 

22 
13¢ 
L8e 
19e 


The Characteristics of 

Ts 

V (°R 
1.55 $02 
6.31 105 
g 70 101 
14.838 404 
18.44 $10 
6.68 398 
10.20 394 
15.50 394 
18.70 116 


2 
a 


th 





e Laminar 


(°R.) (CoV Re 
1,182 0.591 
1,182 0.572 
1,182 0.538 
1,182 0.481 
1,182 0.450 
2,246 0.53 

2 246 0.508 
2 , 246 0.461 
2,246 0.4387 


r) 


TABLE 6 


Boundary Layer at the Nominal Ambient 1 


‘emperature of 400°R. (Heat Transfer Plate 


7 ( Nu ) 6 5* 5°" 
(GV Re)e (PR VPrVRe)s x V RO y V Res ¢ V Res 

0.680 1,120 0.345 10.11 6.77 0.58 
0.657 1,350 0.354 11.61 8.28 0.56 
0.619 1,834 0.361 15.38 i232 .3oa 0.52 
0.554 3,012 0.348 299 49 19.76 0.45 
0.516 $300 0.333 

0.692 1,975 0.347 15.49 12.35 0.53 
0.657 2,520 0.354 19.52 16.31 0.50 
0.599 3,880 0.339 26.338 23.43 0.45 
0.566 $920 0.327 
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TABLE 7 
The Characteristics of the Laminar Boundary Layer at the Nominal Ambient Temperature of 800°R. (Heat Transfer Plate 
Run T5 : i ( ; Nu ) 6 Res 6* 5** 
No. M ("m.) (7H. (GV Reds (CyV Re)u ("m:) VPrv Re/s x ¥ Res x ~ Res rv V Re; 
10d 3.29 798 1,182 0.621 0.661 1,325 0.349 7.29 3.53 0.62 
11d 4.55 805 1,182 0.602 0.640 1,550 0.353 7.94 1.37 0.60 : 
16d 6.99 792 1,182 0.568 0.605 2,250 0.356 9.94 6.58 0.56 ‘ 
17h =: 10.57)—s«812—s«i1, 182 0.513 0.546 3,340 0.344 13.18 10.30 0.50 3 
22d 13.43 784 1,182 0.476 0.510 4,500 0.339 16.67 14.21 0.46 7 
13d 4.77 801 2,246 0.574 0.678 2,145 0.352 9.78 6.438 0.56 ; 
18d 7.22 804 2,246 0.536 0.647 2,684 0.348 11.33 8.40 ().52 E 
19g 10.98 800 2.246 0.488 0.591 4,100 0.337 15.30 12.31 0.49 c 
23e 13.63 792 =2,246 0.462 0.561 5,100 0.329 18.45 15.62 0.46 
THE CHARACTERISTICS OF THE LAMINAR BOUNDARY (f/f) Various Boundary-Layer Thicknesses 
LAYER . Mee : , . 
From the variations of , j, and J '(u/ug) dt with ¢, 
The characteristics of the laminar boundary layer the variations of u/u; and 77/7 with (y/x)V Re; are 
are determined for both the insulated-plate case and the obtained; also the thickness of u/u; = 0.995, the dis. 
heat-transfer case at nominal ambient temperatures of placement thickness, and the momentum. thickness 
°, 400°, : SOO°R. The results are tabulated it T ; f 
100 » 400", and SUU R. The results are hb, ibulated in are determined. The relation between (y/x) V Re; and st 
lables 2 to 7. The methods employed in the trans-  ;;. ; 
formation of the mathematical parameters into the N 
. P ! Nl 
physical parameters are as follows: y vy fosux | “ pupils fn 
— V Re, = — = 4/2 — / di mi 
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Fic. 3. The film temperature for the insulated-plate case. 


A comparison of the film temperatures with the re- 
sults of Rubesin and Johnson’ are shown in Figs. 3 and 
{. The agreements are excellent at moderate Mach 
Numbers. For Mach Number greater than approxi- 
mately 5, the following equations are recommended: 


(a) Insulated-Plate Case 
T'/T; = (1/18) [13(T,/T;) + 6] 
(b) Heat-Transfer Case 
T’/T; = 0.70 + 0.023 M,? + 0.58 (T,,/T;) 


The results for the nominal ambient temperature of 
100°R. indicate higher values of the skin-friction 
parameter, (C, V Re);. This is due to the existence of 
the maximum in the variation of the density-viscosity 
ratio with the enthalpy parameter (see Appendix) at 
low air temperatures. Higher values are also indicated 
for the temperature recovery factor, inherent temper- 
However, the heat- 
This peculiarity 


ature, and film temperature ratio. 
transfer parameter is little affected. 
would serve as a warning that comparisons of hypersonic 
wind-tunnel data with flight-test data should be made 
with caution. 

Since the flow characteristics for 7; = 400°R. com- 
mand the most interest and the trends are similar to 
those for 7; of 100° and SOO°R., the variations of u/1;, 
1/T;, g, and dj/dt for T; = 400°R. are plotted versus 
(v/x) V Re, in Figs. 5 to 16, with Mach Number as the 


parameter. 


CONCLUDING REMARKS 


It is almost self-evident from the correlation of 7”/T; 
with /;, as shown in Fig. 4 for 7; of 400° and SO0°R., 
that thé ‘approximate design procedure of using the in- 
compressible relations for the laminar skin friction and 
evaluating the air properties at the film temperature, 
I’, should give good results at hypersonic, as well as 
supersonic, speeds in the ambient temperature range of 
missile flights. More important, there are indications 
that this procedure can be applied approximately to 


BOUNDARY LAYER 233 












10 "1 
7,5 140.032 my +058( = ') Ap 
(Ref. 9) 
— o 
Oe 
° ’ . 
6 *Pr ,0.704+0.023M, +0.58 7° 
! 
~e O 7. =1182° R 
: 7,=400°R 
Za © 7, =2246°R 
afs 4 T, *1182° R) > 
a’ 8 =B00°R 
a & 7, =2246°R 
a 
re) 1 i 1 
Ce) 100 200 300 400 
Ms! 
Fic. 4. The film temperature for the heat-transfer case 


















s . 
7," 400°R 
0 | 
) 5 10 15 20 25 30 35 40 
; /Res 
Fic. 5. The velocity distribution in the laminar boundary layer 
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The temperature distribution in the laminar boundary 
layer (insulated plate), 


Fic. 6 


the turbulent boundary layer. That is, use the incom- 
pressible relations for the turbulent heat-transfer and 
skin-friction coefficients and the following relations for 


be 


ie Fe P 7 

— = 1 + 0.032M,? + 0.58 (: -— i) (for M < 5.6) 
T's T; 

a pa ¢ - le - . 
r = ().70 + 0.023M,* + 0.58 r (for M > 5.6) 
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laminar boundary layer (7) = 1,182°R.). 
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In Fig. 17 the application of this approximate method 
to the turbulent skin friction over an insulated plate js 
compared with some data as interpreted by Rubesin 
Maydew, and Varga.'? 


AIR PROPERTIES AS A FUNCTION OF 
TEMPERATURE 


APPENDIX 


| The air properties, u, c,, k, p, and h must be known as 
a function of temperature. The beauty of the present 
method of analysis lies in the fact that the best known 
variation of the air properties can be utilized without 
recourse to approximate expressions, relating the air 
properties to the temperature. This is possible because 
the numerical solutions, employing the mechanical 
The 
air properties in dimensionless parameters as used 


differential analyzer, can absorb graphical data. 


herein are given in Fig. 18. 

The enthalpy information is obtained from Keenan 
and Kaye'!' and Krieger and White.'” 
variation is calculated from Sutherland formula—that 


The viscosity 


1S, 


1+0 
PBERB YT Tp r (0/V T Tp) 
or 


u ( ] ) 1+0 

HB Tp VT T (0/V T T,) 
where @ = 0.522. This variation agrees excellently with 
the experimental data of Fortier!* and Vasilesco'! which 
cover a temperature range from 146.7° to 3,321.2°R. 
The variation of the Prandtl Number with temperature 
possesses a degree of uncertainty because of insufficient 
information on the thermal conductivity of air. How- 
ever, fortunately, the effect of this variation plays a 
minor role in the analysis. The Prandtl Number values 
as used are obtained with values of the thermal con- 
ductivity calculated by Hirschfelder and Curtis’ for 
temperatures up to 500°R. In the temperature region 
between 600° to 2,400°R., the values of Keenan and 
Kaye!!! and Tribus and Boelter! are used. At the 
higher temperatures, the Prandtl Number variation is 
extrapolated using the variations of the air components, 
as calculated with kinetic theory, as a guide. The 
present state of the art on thermal conductivity and 
viscosity is presented in Figs. 19 and 20. In Fig. 21, 
the results of using power laws for the viscosity vari- 
ation are compared with the Sutherland formula 
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Experiments on Interaction of Shock Waves 
and Cylindrical Bodies at Supersonic Speeds’ 


M. V. MORKOVIN,?+ E. MIGOTSKY,+ H. E. BAILEY,** ano R. E. PHINNEY** 
Unwersity of Michigan 


ABSTRACT 


While two-dimensional shock-wave reflections have received 
much attention, in most applications either the shock surface or 
the reflecting surface is curved. Elements of nonviscous theory 
for three-dimensional shock reflections by two of the present 
authors indicated the likelihood of anomalous reflection for con- 
figurations such as that of a plane shock impinging on a cylinder 
aligned with the free stream. This paper presents the high- 
lights of the experimental investigation of this three-dimensional 
shock-cylinder configuration conducted in the University of 
Michigan Supersonic Wind Tunnel at a Mach Number of 
1.90. 

One general result is the striking boundary-layer cross flow 
induced by the sharp pressure gradients across the reflection line 
on body surfaces. For the present configuration, the cross flow, 
in and outside of the boundary layer, results in a separated re- 
gion on the lee side of the body. The effects of boundary-layer 
character and thickness and of shock strength on the overall 
flow and pressure fields are studied by means of body-pressure 
distributions and several visualization techniques. Implica- 
tions of the present findings for practical applications are dis- 


cussed. 
SYMBOLS 

6 = flow deflection angle across incident shock in the vicin- 
ity of the cylindrical body 

w» = meridian angle 

@ = shock-wave inclination angle 

M = free-stream Mach Number 

p = static pressure on surface of cylindrical body 

pf = free-stream stagnation pressure 

Ap = pressure rise on cylindrical body due to presence of 
shock wave 

II’ = trace of shock plane 

é = subscript, indicatiag that the variable is in the plane 
normal to the local direction of shock-cylinder in- 
tersection 

(r:) = boundary layer artificially roughened by paint line 

(rz) = boundary layer artificially roughened by two roughness 
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INTRODUCTION 


i THE LAST 5 YEARS much attention has been 
given to the study of the two-dimensional 
problem of supersonic shock wave reflections from a flat 
surface with boundary layer both analytically'~* and 
experimentally.°~'! Although the problem is far from 
solved, some understanding of the basic variables in- 
fluencing the phenomenon has been achieved, notably 
through references 1, 9, and 10. 


It is generally accepted that the nonviscous theory of 
so-called regular two-dimensional shock reflection’? 
yields the correct values of the pressure jump and the 
shock-wave direction at distances sufficiently far from 
the point of impingment of the shock on the wall. 
Near this point, both upstream and downstream, the 
boundary layer is strongly affected by the shock wave. 
When the boundary layer is laminar, the resulting ad- 
verse pressure gradient is felt far upstream (as much 
as 100 times the undisturbed boundary-layer thick- 
ness,*’'’) and the layer is believed to separate for all but 
the weakest shock waves. The incoming shock re- 
flects as an expansion wave from the free separated sur- 
face and turns it back to the wall, forming a bubble. 
The outgoing compression waves generated by the 
thickening and separation, as well as by the reattach- 
ment of the boundary layer, tend to coalesce with the 
expansion wave and to form a single sharp compres- 
sion wave (the reflected shock) far enough from the 
wall. For a given Mach Number, the interaction be- 
tween the laminar boundary layer and the shock in- 
creases as the shock strength increases and the bound- 


ary layer thickness decreases.'» '° 


Reflection from a wall with a turbulent boundary 
layer is less well understood. For weaker shocks,?~'! 
the boundary layer appears to thicken more or less 
rapidly just upstream of the impingment point, and the 
wall pressure measurements indicate a sharp rise with a 
relatively small upstream influence. For stronger 
shocks,'' local separation and, finally, bifurcation may 
occur with the corresponding shock-expansion-shock 
pattern and increasing upstream influence. The ranges 
in which these régimes occur seem to depend also on 
the Mach Number and the boundary-layer thickness. 
It is believed that the basic difference lies in the fact 
that a given adverse pressure gradient produces a 
much smaller rate of growth when the boundary layer 
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is turbulent than when it is laminar. Reference 9 sional nature of the problem, a rather complex three. In 
demonstrated that the observed difference between the dimensional flow behind the shock waves would ensye es 
laminar and turbulent cases cannot be ascribed to the To solve completely the problem of nonviscous three. ed . 
relative thickness of the subsonic portion of the layers. dimensional shock reflection, one would have to com. ie 
In absence of published data, however, it is believed pute this postreflection flow, which in the two-dimep. shocl 
that the effect of an increase in thickness of a turbulent sional case was merely a uniform parallel flow. id 
boundary layer could be explained by the role of the One of the cases studied analytically in reference 13 we 
increased thickness of the subsonic layer as an avenue was that of a plane shock wave impinging on an infinite - 
for upstream influence. cylinder, which is aligned with the steady free stream. cous 
In practice, the configurations in which a shock wave — It was found that at all supersonic Mach Numbers and Fig 
falls upon a solid boundary are generally three-dimen- at all shock orientations there are regions on the reflect. the f 
sional—.e., either the shock wave or the reflecting (in- ing cylinder where a regular (nonviscous) shock re. to tl 
terfering) surface is curved in such a manner that the flection cannot occur. If the shock did not warp near the 
line of their intersection is not everywhere normal to the these irregular regions, there would also be regions on angl 
oncoming stream lines. It follows that there are at the body where no reflection at all is possible. Clearly, meri 
least two effects not present in the two-dimensional according to the nonviscous theory, the three-dimen- nate 
problem heretofore analyzed: (a) The flow field down- sional nature of the body or shock must induce some sn 
stream of the reflecting wave is not uniform and gener- changes, such as warping, in the incident shock wave to cl 
ally involves a postreflection cross flow (relative to pre- before it reaches the surface. Such a behavior is again to t 
shock stream-line direction) even in the nonviscous basically different from that in the two-dimensional axis. 
case; (b) the (vector) pressure gradient affecting the case. Should the warping be in the form of Mach Y poin 
boundary layer may have an appreciable component type of reflection, pockets of subsonic flow would be onus 
normal to the preshock stream lines; where the created so that the reflection phenomenon would line 
velocities inside the boundary layer are small, one may become more directly dependent on the conditions Nur 
expect a boundary-layer cross flow considerably downstream. Also, any warping whatsoever would A 
stronger than the outer postreflection cross flow (a). make the analytical solution of the postreflection flow eviii 
Two of the present authors established’*® the link be- nearly impossible. regi 
tween the two-dimensional and the three-dimensional It appeared likely that the shock interference with low 
nonviscous reflection theory, the so-called local sweep- the boundary layer, in particular the boundary-layer tion 
back principle. Since in both cases the reflection is cross flow already mentioned, could modify the condi- sibl 
governed by the condition of no velocity normal to the — tions of reflection appreciably and possibly mask the ese 
body and since locally the flow changes (singularity anomalous reflection predicted by the nonviscous the 
type) across incident and reflected shocks literally theory. Since the plane-shock cylinder configuration be | 
dwarf any flow changes due to continuous curvature of also represents the simplest case of the important wing- F 
stream lines or shock waves, one needs merely to con- body interference problem (wing influencing body pro 
sider the local two-dimensional reflection problem in without reciprocal influence of body on wing), it was de- is t 
the effective plane, defined as the plane normal to the cided to explore this case of three-dimensional shock ain 
intersection of the incident shock with the body. reflection experimentally. inte 
The result of such point-by-point solutions along the ncl 
intersection line would yield initial values of velocity EXPERIMENTAL OBJECTIVES AND TECHNIQUES The 
on the surface of the body; since these initial values of , a 
velocity are not constant because of the three-dimen- The experiments were conducted in the 1.90 Mach lar 
Number channel of the University of Michigan Super- a 
- sonic Wind Tunnel. This tunnel, which is intermittent on 
auiibhens ye and has an 8- by 13-in. test section, is described in ; 
. | yr reference 14, and its detailed calibration is reported im 7” 
we | nt reference 15. The model configuration can be seen in ae 
)7 idl Figs. 2 and 8; problems connected with model design 
6 Pe and instrumentation are discussed in detail in refer- -_ 
° 1" AZ A 7 ences 16 and 17. . 
o 14 24 Fi Broadly speaking, the objectives of the experiments a 
vs Kg were, by various visualization techniques and pressure 
e Z SAS | | , 2.0 measurements: (1) to explore the regions on the cyl- 
fe JS 5 — inder where anomalous reflection is indicated by the 
Z Pog Ls arc a nonviscous theory;'® (2) to obtain the main charac- 
10 2 1.4 1.6 I 2.0 22 24 26 teristics of the boundary-layer and postreflection cross 
EFFECTIVE MACH NUMBER, My ag OI eae mn ; 
a a ls Rated flows; (3) to assess various effects on (1) and (2) due to 
upon nonviscous theory. changes of the boundary layer and of the shock strength. 
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INTERACTION OF SHOCK 


In pursuing objective (1), a special technique was de- 
signed to check optically a point in the region of 
anomalous reflection which otherwise cannot be seen. 
This was accomplished by simultaneously rotating the 
shock system by 25° around the axis of the cylinder 


25° from its normal 


and yawing the schlieren beam by 25 
position (see Fig. 9a). The choice of these angles can 
be understood in the light of the results of the nonvis- 
cous theory of reference 13 which are summarized in 
Fig. 1. The abscissa, M,, represents the projection of 
the free-stream Mach Number, J/, on the plane normal 
to the line of intersection between the shock plane and 
The angle 0, is the effective shock wave 

Thus, at the w = 90° 
i.e., the ordi- 


the cylinder. 
angle in this normal plane. 
meridian, //, = M sin 6 and @, = 90° 
nate csc 06, = 1; as free-stream Mach Number  var- 
ies or as the strength of the incident shock varies due 
to change of wedge angle 6, the points Py) corresponding 
to the 90° meridian are represented on the horizontal 
axis. In the same sense, shock conditions at any 
point on the incident side of the cylinder will be repre- 
sented by a point in the sector between the straight 
line VM, sin 6, = 1 and the horizontal axis for all Mach 
Numbers and shock strengths. 

According to the nonviscous theory, any region on the 
cylinder corresponding to points above the curve A is a 
region where regular reflection is possible. Points be- 
low curve A correspond to regions where regular reflec- 
tion of a truly plane incident shock becomes impos- 
sible. Finally, points between curve B and the axis 
esc 6, = 1 refer to the region on the cylinder where, if 
the shock could remain plane, no reflection at all would 
be possible. 

For a given Mach Number and shock strength, the 
projection of the W/ vector normally to the shock plane 
is the same at all points and is equal to /, sin #é, = M 
sin 6. Thus, as one proceeds along the shock-cylinder 
intersection, one covers a segment of one of the family of 
inclined straight lines 7, sin 0, = constant in Fig. 1. 
The starting point of the segment corresponding to 
w = ( may be anywhere, say above curve A where regu- 
lar reflection is possible; the segment must then always 
cross curves A and B and end on the horizontal axis, 
which corresponds to w = 90 

The segments corresponding to three of the cases 











studied by pressure measurements at JJ = 1.90 
namely, 0 = 35.2°, 41.8°, and 42.9° (effective wedge 
- 7.00" - 831 - 4.19°—+ 
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Model configuration 
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Fic. 3a. Pressure distribution on cylindrical body. w = 0°; 
6 = 10.2° 


angles 4.0°, 10.2°, and 11.2°)—are shown dotted in 


Fig. 1. In the case 6 = 10.2°, one sees from points cor- 
responding to 0° and 25° meridian, Py) and P»5, that the 
shock can have regular reflection at the normally ob- 
served meridian w = 0° but that the nonviscous 
theory indicates anomalous reflection at the meridian 
w = 25°, which was made visible by the yawed-schlieren 
system. 

Ordinary schlieren photographs with different knife- 
edge positions were taken regularly. In addition, the 
flow was studied by shadowgraphs and the light-scat- 
tering technique recently used by N.A.C.A. (so-called 
vapor technique). In pursuing objective (2) a special 
mixture of china-clay, oil of wintergreen, and iso-safrole 
was developed by Murphy,' of the tunnel staff. A 
thin film of this colorless mixture applied to the body 
evaporates at the proper rate during the 20-sec. tunnel 
runs so that the flow lines in the inner portion of the 
boundary layer become visible (Figs. 11b, 12b, 13b, 15b, 
and 16). The tuft technique was also used for indica- 
tion of average flow directions near the surface. 

The location of static pressure orifices (0.025 in. in 
diameter) is shown in Fig. 2. The cylinder was moved 
fore and aft and was rotated about its axis in the cases 
6 = 4.0° and 11.2° in order to obtain as complete a 
coverage of the shock neighborhood as desired. The 
model axial position for which orifice No. 1. fell 0.29 in 
forward of the centerline of the tunnel window was 


arbitrarily designated as standard axtal position, and all 
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Fic. 3b. Pressure distribution on cylindrical body. w = 90 
6 = 10.2°. 


axial displacements were measured and recorded with 
respect to this standard. At first, a rotatable wedge, 
W,, 3.50 in. wide and 3.96 in. long, was used as the 
shock generating unit. As the blocking characteristics 
of this particular configuration in the tunnel became 
clearer, a second wedge, W (5.50 in. wide, 5.46 in. long, 
and 0.40 in. closer to the body), furnished a larger test- 
ing region in the vicinity of the shock not appreciably 
affected by the outside interference such as that due to 
wedge tips, model nose and shoulder, etc. The tunnel 
and cylinder interference effects are essentially elimi- 
nated by considering the pressure increments, Ap, of 
body in presence of wedge over the body alone con- 
figuration (see, for instance, Fig. 3a). 

In order to bring out the overall effects of large 
changes in body cross section, some experimental in- 
vestigations of interaction between the identical plane 
shock wave and a long square body were made. The 
length of the bedy, shoulder to base, was 12.50 in., and 
the depth and width were 1 in. The pyramidal nose 
was 7 in. long and had a half angle of 4°5’. Fig. 15a 
shows the square body and a typical interaction pat- 
tern. 

The Reynolds Number at the standard plane on the 
body was 5.43 X& 10°. 
obtainable at this station. 
to 5 in. aft of the standard position, optical studies 


Laminar boundary layer was not 
By shifting the cylinder 3 


(without pressure distributions) of the shock-wave inter- 
ference with a laminar boundary layer were made (Figs. 
9a and 9b). The condition of turbulent boundary 
layer near the pressure orifices was designated as 
“smooth” when the transition point was not controlled 
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artificially, rough (7;) when a narrow line was painted 
on the body 1 in. downstream of the cone-cylinder 
shoulder, and rough (72) when two rings (wire diameter 
0.022 in.) were slipped on the cone (Fig. 8). For runs 
in the smooth condition there was an unpredictable* 
variation of the transition point (as much as 3! > in. 
and a consequent variation in the boundary-layer thick. 
ness which in turn was responsible for a larger scatter of 
pressure measurements in the vicinity of shock jm. 
pingement (Fig. 3a). The paint line in the rough (r,) 
condition appeared to bring about a uniform boundary 
layer judging by the elimination of scatter in the total 
The total 
head pressure probes were made by flattening 0.065-in, 


head measurements on the body surface. 


outside diameter stainless-steel tubes to give an open- 
ing having an inside dimension of 0.007 in. The probes 
were soldered to the rear side of the body (Fig. 11b) 
so that they would not influence the static pressure 
measurements on the front side. 

In order to obtain the effective shock deflection angles 
6, the negatives of the various schlieren photographs 
were projected onto a ground glass, and the shock in- 
clination near the cylinder relative to the trace w = 0) 
was read several times by different people on a special 
protractor. From the average value of 0, the deflection 
angle 6 was deduced by using graphs from reference 19. 
Such a procedure may lead to errors of as much as 
+(0.30° in the present case. 

In the preceding description of the experimental 
techniques, the likely sources of errors were brought 
out. It is believed that the errors actually present do 
not influence the understanding of the basic effects that 
were selected for study in objectives (1), (2), and (3). 


THE THREE-DIMENSIONAL FLOW FIELD 


The problem of shock-wave diffraction around the 
cylinder can also be visualized as a pressure and flow 
field response to a sudden (step) application of a con- 
stant angle of attack to the portion of the cylinder 
downstream from the inclined shock plane. According 
to the linearized view of the problem, there would be 
local distortions of the flow and pressure fields in the 
vicinity of the “‘loading’”’ line on the body; these distor- 
tions would disappear with increasing distance down- 
stream, until ultimately the flow would settle down to 
the (doublet) flow around an infinite circular cylinder 
at the given angle of attack. In the nonlinear view, the 
distortions of the field near the loading line may well 
be greater and may include the effects of Mach Y forma- 
tion and warping of the shock wave out of the incident 
plane. Viscosity is likely to modify the field distortion 
by allowing distortion of the effective flow boundary by 


virtue of boundary-layer changes. It is likely that 





* It is likely that the phenomenon is connected with variable 
heat transfer from the body caused by differences in the tem- 
perature of the tunnel air supply and the body at the beginning 
of the run. Despite much effort, no conclusive correlation has 


as yet been obtained. 
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INTERACTION OF 


even in the nonlinear viscous case the local field distor- 
tions will die down with increasing distance downstream 
and the flow will settle down to the viscous flow around 
an infinite cylinder at the same angle of attack (yaw). 
The distortions of the field due to the boundary 
laver are twofold. First, there is the essentially two- 
dimensional effect of the upstream propagation of the 
pressure rise associated with the shock as discussed in 
the Introduction. This effect, of course, modifies 
basically the boundary layer as a whole (see for in- 
stance the boundary-layer edge in Fig. 9b). Second, a 
boundary-layer flow across meridional planes is gene- 
rated, thus further modifying the effective boundary 
surface. The inclination (relative to the free stream) 
of the flow lines in the boundary layer (Figs. 11b, 12b, 
and 13b) on the side of the body appears to be always 
somewhat less than the corresponding shock-wave angle 
but is several times as large as the angle of deflection 
across the shock outside the boundary layer. This 
flow-line inclination decreases with distance downstreain 
from the shock but the boundary-layer cross flow per- 
sists. The corresponding partial disappearance of the 
boundary layer from the w = O° meridian and its ac- 
cumulation on the w = 180° meridian are clearly visible 
in all schlieren photographs with horizontal knife edge. 
It is notable that the overall features of the field 
just described are all present in the case of shock dif- 
fraction around a square body (Fig. 15) and are there- 
fore likely to occur for most convex elongated bodies. 
Fig. 15b of the flow lines on the side of the square body 
shows that while these are inclined at slightly over 40 
(corresponding to a shock angle of 42.9°), the line along 
which the flow lines start sharply curving is inclined 
nearly 60° to the free stream. This illustrates the ef- 
fect of the cross flow as an “‘increasing obstacle’’ to the 
oncoming flow as one proceeds along the shock-body in- 
tersection line and of the observable forward warping of 
the shock wave (Fig. 15a). Both of these effects can be 
observed on schlieren photographs of the cylindrical 
body such as Fig. 12a. The shock continuation be- 
yond the body on the lee side indicates some forward 
warping; the continuous compressive wave, which 
coalesces into a shock upstream of the main shock con- 
tinuation on the lee side, indicates the boundary-layer 
effect. This forward effect at an angle greater than the 
shock angle is naturally evidenced by the higher in 
clination of the line along which the pressure rise be- 
gins on the composite pressure distributions of Figs. 4 


and 5. 


LEE-SIDE EFFECTS 


The phenomenon of shock reflection on the lee side 
(90°< w < 270°) is fundamentally different from that 
on the incident side ( w < 90°) even in the nonviscous 
case. Thus, on the incident side, the oncoming shock 
wave induces velocities that are directed into the sur- 
face of the body, and it can be considered that their 
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Composite arrangement of pressure distribution 
6 = 4.0". 


Fic. 4. 


required cancellation brings about the reflected shock 
wave. As the shock wave spreads toward the lee side 
after passing the w = +90 
duced velocities have a component away from the body. 
Even if one knew the line along which the shock inter- 
sects the body on the lee side, it is not clear physically 
that these outward radial velocity components need be 
completely canceled and that the flow must continue 
to follow the surface, especially in the viscous case. 
The cancellation of the outward component or any part 
of it would have to be effected by an expansion wave 
The expansion 


meridians, the shock-in- 


originating along the intersection line. 
wave would tend to weaken the shock wave which 
generated it and form a ‘‘shadow”’ region (partly or 
completely shielded from the oncoming shock). 
Shadowgraphs (Figs. 9b and 10), schlieren photo- 
graphs (Figs. 8, 12a, and 14), and irreproducible pic- 
tures obtained by the vapor technique show that the 
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incident shock is indeed broken up and does not reform 
in the w = 180° plane for a distance of one to three 
diameters from the body. The regeneration of the 
shock on the lee side is undoubtedly connected with the 
thickening and accumulation of the boundary layer. 
In the laminar case the boundary layer can be clearly 
seen to curve away from the body at some distance up- 
stream from the plane of the shock in Fig. 9b. Simi- 
larly, in the turbulent case, schlieren photographs 
(Figs. 8, lla, 12a, and 13a) indicate that the outside 
flow in the vertical plane w = 180° deflects away from 
the cylinder surface at an angle somewhat less than the 
deflection angle across the oncoming shock (but nearly 
proportional to it). The wedgelike space thus created 
is apparently occupied by slower moving air made up of 
the original boundary layer on the lee side and from air 
swept in from the bottom and sides of the cylinder by 
the previously mentioned cross flow of the boundary 
layer. 

The separation on the lee side thus seems to be akin 
to that which may occur on bodies with sweepback.” 
Since the shock wave induces a crosswise component of 
flow, this similarity is not surprising. As a matter of 
fact, the flow resembles that over bodies of revolution 
at moderate angles of attack (except in the immediate 
vicinity of the shock wave), and there is some evidence 
of corresponding vortex formation discovered by Al- 
len.”! 
ary-layer flow lines on the lee side of the cylinder (6 = 
14.3°), and Fig. 15b shows a similar pattern on top of 
the square body (6 = 11.2°). The flow lines neck in 
where the shock influence is felt first and then diverge 


Fig. 13b exhibits a typical formation of bound- 


It is not clear 
whether the two regularly appearing flow lines, which 


from the w = 180° meridian again. 


seem to absorb others and become rather strong, mark 
the line of separation on the side of the body or not. 


Once the lee-side separation appears it seems to per- 
sist, even when an expansion wave with the opposite 
pressure gradient affects it (Fig. 8). In general, an ex- 
pansion wave falling upon the body produces practi- 
cally no boundary-layer cross flow in the opposite 
direction, Since the local cross component of the 
pressure gradient in an expansion wave some distance 
from the expansion corner is a mere fraction of the cor- 
responding contponent due to a shock wave, this is 
readily understandable. 

From the preceding observations it appears likely 
that there will be a tendency toward lee-side separation 
on bodies of revolution (and elongated convex bodies in 
general), even when the oblique shock is generated by a 
Whether the 
strength of the boundary-layer cross flow is materially 


wedge abutting on the side of the body. 


changed when the ratio of the body diameter to the 
boundary-layer thickness is increased as on full-scale 
prototypes remains to be seen. The design of winged 
models for tunnel investigation of those aspects of super- 
interference associated with the 


sonic wing-body 
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boundary-layer build-up and its effect in the region of 
tail surfaces is now in progress at the University of 
Michigan. 


PRESSURE FIELD 


Figs. 3a and 3b exemplify the effects of boundary. 
layer thickness on pressure distributions at w = 0° and 
90° meridional stations for a moderate shock strength 
(6 = 10.2°). For the thicker boundary layer, the net 
pressure rise in Fig. 3a is smaller, but the shock is felt 
This latter behavior is in qualitative 
agreement with the concept of the subsonic sublayer 
(which also thickens) as avenue of upstream propaga- 
tion. On the side of the body (w = 90°), the smooth 
and rough (72) boundary-layer conditions unexpectedly 


farther upstream. 


lead to different pressure levels (Fig. 3b); this effect is 
probably caused by the difference in the boundary- 
layer cross flow in the two cases. Since it is known 
that the pressure distributions at w = 90° rise again 
farther downstream (Fig. 5), the pressure may ullti- 
mately reach the same level in the two cases. 

The considerations concerning the pressures at w = 
90° have implications for pressure measurements by 
long-needle probes and static pitot tubes. It has been 
generally agreed that near shock waves the pressure 
orifices on the probes should be located at w = 90 
rather than w = 0° or 180°. 
strong influence of the nature of the probe boundary 


Reference 9 reports a 


layer on the pressure measurements near shock waves 
and recommends (with caution) making the boundary 
layer turbulent. 
be extended to the small probe scale, even these latter 


If the present larger scale results can 


measurements are subject to a considerable degree of 
distortion. Fig. 14 exhibits the similarity between 
the shock-interaction patterns for a probe and the pres- 
ent model. In particular, the reflected wave at w = 0 
and the lee-side separation are indicated. 

The pressure field around the cylinder at 6 = 4.0 
and 11.2°, together with the shock trace //’, is repre- 
sented in the spatial composite pressure distributions 
of Figs. 4 and 5. The reader can interpret most of the 
features of these distributions in terms of the concepts 
of the boundary-layer cross flow and the postreflection 
The latter flow tends to relieve the over- 
As al- 


ready discussed, details of these pressure distributions 


outer flow. 
pressure on the incident side rather rapidly. 


depend on the nature and thickness of the boundary 
layer and on shock strength. The maximal values of 
Ap/ po are shown irrespective of their longitudinal post- 
tion, in Fig. 7, as a function of w. For comparison, 
pressure jumps based on nonlinear nonviscous theory"’ 
are shown in Fig. 7 for 0 < w <w, and for w = 90 
(point C). 


well as weak shocks, the maximum pressure jumps re- 


In the case of relatively strong shocks, as 


corded fall short of the values predicted by the theory 
even at w = 0°. The pressure jump on the w = 90 


meridian is also much smaller than the single jump 
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across the oncoming shock wave (point C). This 
smallness of the wall pressure jumps can probably be 
ascribed to the effect of a relatively thick boundary 
laver that could accommodate sizable radial pressure 
gradients and to the three-dimensional relieving effect 
of the cross flow. The latter effect is probably re- 
sponsible for shifting the smallest values of pressure 
jumps to the vicinity of w = 120° meridian. Two 
aspects of the overall pressure variation deserve closer 
attention. 

(1) The forward bump in the pressure distributions 
on the incident side of Fig. 5 is associated with the 
boundary-layer cross flow and has some characteristics 
of a dead-air bubble.* In Fig. 6, the static pressure 
variations at w = 55° and 63° are matched with the 
surface total-head variations at w = —50° and —59°. 
It would appear that at the location of the bump the 
component of flow velocity against the total-head probe 
opening has vanished so that the static and total pres- 
sure values are essentially the same. For the rela- 
tively thick boundary layer and relatively strong shock 
of 6 = 11.2°, the same behavior seems to take place 
over most of the incident side; for thinner boundary 
layers and weaker shocks, the static pressure bump and 
the total-head dip develop only for higher w’s on the in- 
cident side of the body. For instance, at 6 = 10.2' 
with smooth turbulent boundary layer (Fig. 3a), the 
bump is absent and develops only somewhere between 










































w = 25° and w = 55°. For the smooth turbulent 
* See also next section. 
44 7 
40}- —--—+—_--—} 
36} jt 
& 
— ae — 
a 
° / | | STATIC 
_ | / | 
ae eee ee ee ee 
}: 
W 
« 
- 
a re 
. “TOTAL HEAD 
2 (BODY ALONE) 
63° -59° 
55° -50° 
> 
AZ 
1.6 1.2 0.8 0.4 STD 0.4 0.8 


AXIAL POSITION , INCHES 
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Fic. 7. Experimental and theoretical variation of pressure ratio 
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Fic. 8. Overall view of shock reflection pattern from circular 
body. Horizontal knife edge; 6 = 10.2°; rough (re) condi- 
tion. 
boundary layer at 6 = 4.0° the effect is not discernible. 


It is interesting to note (Fig. 5) how this bump, asso- 
ciated with the boundary-layer cross flow, gradually 
becomes the main pressure peak as w increases past 
80°—1i.e., while the lee side is partly shielded from the 
direct shock, a relatively forward pressure rise occurs 
there due to the boundary-layer cross flow. 

(2) Comparison of pressure variations on the lee 
side shows that, irrespective of the boundary-layer 
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thickness and shock strength, the pressure starts rising 
again slowly as one crosses the meridional planes s 
ward w = 180°. This is consistent with the settling 
down of the flow to that around the cylinder at an 
angle of attack. This latter flow exhibits a similar ad. 
verse pressure gradient in both the potential and yis. 
cous case.”! 


COMPARISON WITH THEORY AND WITH 
Two-DIMENSIONAL EXPERIMENTS 


The optical data obtained in this investigation show 
that the anomalous effects predicted by the nonviscous 
theory'® are masked by the boundary layer so far as 
can be observed at w = 25° for the critical case § = FIG 
10.2° (see Fig. 1), irrespective of the laminar or turbu- 
lent nature of the boundary layer. In the laminar 
case no evidence of a Mach Y type of reflection was ob- 
tained (Figs. 9a and 9b). This is not surprising inas- 
much as the protuberance-type geometry of the 
laminar boundary layer in the neighborhood of the in- 
tersection changes the local requirements on the de- 
flection angles which were the critical conditions for the 
existence of a regular reflection in the nonviscous theory. 
In the turbulent case there was also no significant dif- 
ference in the basic shock reflection pattern between 
w = 0° andw = 25°. It had been conjectured" that a 
Mach Y type of reflection might occur in the predicted 
region of anomalous reflection; however, for a given 
boundary layer, no Mach Y configuration could be as- 


cribed to the change from w = 0° tow = 25° location.* 
As a matter of fact, for a given shock strength any ob- 
servable Mach Y shocks appear to be related primarily - 
to increase in thickness of the boundary layer. [There 20) 


seems to be no clear-cut criteria for distinguishing be- 
tween bifurcated regular reflections and Mach Y re- 
flections from either schlieren photographs or wall 
pressure measurements probably because the transi- 
tion is rather gradual (see Figs. 13c, 13d, 23, and 32 of 
reference 9) ]. 


It is interesting to note that, while there is boundary- 
layer cross flow and therefore a modification of the ef- 
fective boundary of the three-dimensional body at all 
shock strengths, the pressure bump associated with the 
cross flow does develop only at higher w's on the inci- 
dent side when the shocks are weaker. It is difficult to F 
determine the demarcation line between the peak type 
(Fig. 3a) and the bump-plateau type of pressure distri- 


butions, especially since it is likely to depend on bound- ss 
Vv 
* The schlieren system could not be yawed more than 25° t 
because of tunnel window buckets, and thus, unfortunately, no t! 
meridians beyond w = 25° could be observed v 
Pp 
0 
a 
Fic.9. Typical shock reflection pattern from circular body with 1 
P wee . ° = a 9° 
laminar boundary layer viewed in w = 25° plane; 6 = 10.2°. 


Top (a) Vertical knife-edge schlieren. Bottom (b) Shadow graph. a 
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(Figs. 9a and 9b) is similar to that obtained for the two- 
dimensional case*~' as described in the “ Introduction.”’ 
However, in the present case, the edge of the boundary 
layer has a greater curvature corresponding to a 
stronger envelope-type reflected shock upstream of the 
impingement line. Whether the explanation is to be 
sought in the three-dimensionality of the configuration 
or in the greater shock strength is not clear. 

In the case of the turbulent boundary layer, the re- 
flection patterns obtained are in fair agreement with 
those of reference 11. For weak shocks (Fig. lla) the 
reflected shock seems to be followed by a weak expan- 
sion wave and another faint compression wave. This 
behavior, as well as the sharpness in the corresponding 
Fic. 10. Typical eo reflection pattern from pressure rise (Fig. 4, w = 0°), agrees also with the re- 

laa sults of references 9 and 10. As shock strength is in- 
creased, the upstream influence (which also increases 
with boundary-layer thickness) grows and a configura- 
tion of an evelope-type shock, an expansion fan, and a 





second envelope-type compression results (Figs. 8, 
12a, 13a, 14). In reference 11 (where the optical re- 
solution was unimpaired by the finite width of the re- 
flecting body) this configuration was associated with 





Fic. lla. Typical shock reflection pattern from circular 
body with turbulent boundary layer; smooth condition; hori- 
zontal knife-edge schlieren; 6 = 4.0°. 








Fic. 12a. Typical shock reflection pattern from circular body 
with turbulent boundary layer; rough (7;) condition; horizontal 
knife-edge schlieren; 6 = 11,2°. 





Fic. 1lb. Surface stream-line patterns on circular body; 6 = 
4.0°, 


ary-layer thickness. From the present experiments it 
would appear to lie in the region of anomalous reflec- 
tion'® between the curves A and B of Fig. 1. Whether 
the development of the bump is actually associated 
with the impossibility of regular shock reflection (in 
presence of the boundary layer) on the side of the body 
or not will remain for speculation. 

The reflection pattern of the shock wave (6 = 


10.2°) from the incident side of the body (both w = 0 





and w = 25°) in the case of the laminar boundary layer Fic. 12b. Surface stream-line patterns; 6 = 11.2 
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Fic. 13a. Typical shock reflection pattern from circular body 
with turbulent boundary layer; rough (7) condition; horizontal 
knife-edge schlieren; 6 = 14.5°. 





13b. Surface stream-line pattern on lee side of circular 
body; 6 = 14.5°. 





Fic. 14. Shock reflection patterns from circular body and 
probe; smooth condition; horizontal knife-edge schlieren; 6 = 
10.2°. 


local separation (bubble) of the turbulent boundary 
layer. In the case of the square body (6 = 11.2°) this 
bubble can also be seen (Fig. 15a); the corresponding 
boundary-layer flow lines (irreproducible in print) cor- 
roborate its existence. 

There is a suggestion of a turbulent boundary-layer 
bubble on the cylinder in most schlieren photographs 
with horizontal knife edge (Figs. 8, 12a, 13a, 14, etc.) 
for 6 = 10°.2. 


to the presence of a bubble at w = 0° or w = 25° on the 


The evidence of the flow lines relative 


cylinder is inconclusive at all 6's observed. 
Shadowgraphs such as Fig. 10 indicate that the second 

reflected gradual. A 

similar conclusion can be reached from the second re- 


compressive wave is rather 
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flection of the reflected shocks from the rear of the 
wedge (Fig. 12b). 

It must be remembered that three-dimensional] ¢f. 
fects, in particular, the cross flow that thins the bound. 
ary layer downstream from the point of incidence at 
w = 0° may influence the comparison (including the 
two-shock configurations) with other investigations. 
The fact that a clear-cut bifurcation causing separation 
without reattachment!! was not obtained at the strong- 
est shocks (6 = 14.3°) may be thereby explained. (Op 
the other hand, bifurcated reflections in two-dimen. 
sional passages are strongly influenced by the down- 
stream characteristics of the diffusion process and may 
be no more representative of the local shock-b« yundary- 
layer conditions. ) 


APPLICATIONS TO WIND-TUNNEL TESTING 


In Fig. 11a one can discern several weak shocks (prac- 
tically Mach lines) which appear to issue from the sur- 
face of the cylinder upstream of the shock incidence. 
Three pairs of these shocklets appear more or less dis- 
tinctly in all the schlieren photographs in the same rela- 
tive position with respect to the tunnel, irrespective of 
[A fourth pair of Mach lines, 
which move with the body, is at times observable (Fig. 


body axial position. 


13a); these are caused by the occasionally loosened 
juncture between parts of the model.] The stationary 
Mach lines are connected with the very problem of re- 
flection of (weak) plane shocks from a cylinder studied 
in this report. construction of 


the reflected wave reveals that the reflected rays form a 


A simple, “‘ray-by-ray”’ 


three-dimensional (nearly conical) surface; the shape 
of stronger reflected waves observed by the vapor tech- 
nique and from their interception on wedge (Fig. 12b) 
confirms the construction. The three pairs of station- 
ary Mach lines are really profiles of three-dimensional 
reflected Mach surface brought out by the schlieren 
light beam normal to the sidewalls. The plane inci- 
dent waves and their reflections from the square body 
travel in vertical planes from sidewall to sidewall and 
cannot therefore be picked up directly by the schlieren 
system (Fig. 15a). 

The fact that weak nonobservable sidewall disturb- 
ances can be located by their visible reflections from 
cylindrical surfaces can be used as an effective tool in 
The 
actual waves whose reflections are observed in Fig. 8 


tunnel calibration and tunnel-flow improvement. 


are extremely weak; they have not been identified dur- 
(In- 
cidentally, as in most calibrations, much use was made 


ing the original calibration of the wind-tunnel. 


of static pressure probing, which is subject to additional 
error near the weak waves, as previously discussed, and 
may mask the presence of the disturbances. ) 

The visualization technique used here to establish 
the cross-flow characteristics can be employed effec- 
tively for qualitative studies of most test configura- 


tions. As an example, it simultaneously gives typical 








info 
wed 
circt 
Fig. 
gra} 
(6 = 
face 
stre 
tion 
the 
(3) 
bou 
the 
tral 
tio1 
con 
side 
anc 
by 
for’ 
teri 
sho 
ary 
bul 
7 
to | 
tes 
un 
pre 
mo 
str 
gra 
the 
cre 
fie! 
tio 
me 


th 


bo 
sh 
ad 
sic 
fa 


Mt the 


al ef. 
ound- 
ice at 
g the 
tions, 
ation 
rong- 

(On 
men- 
own- 
may 
lary- 


n 
n 











INTERACTION OF SHOCK WAVES & CYLINDRICAL BODIES 247 


information about the flow around the tips of a finite 
wedge, the interaction of a 7° expansion wave with a 
circular cylinder, and tunnel blocking characteristics in 
Fig. 16. With the aid of corresponding schlieren photo- 
graph, one can conjecture that (1) the high pressure 
(§ = 16°) from the lower surface affected the upper sur- 
face at angles from tips slightly higher than the free- 
stream Mach angle; (2) there was some flow separa- 
tion at the squarely cut tips of the wedge (which raises 
the question of proper tip design at supersonic speeds) ; 
(3) the expansion wave produced an insignificant 
boundary-layer cross flow on the body as expected from 
the weak pressure gradients; (4) it was primarily the 
trailing-edge shock that (unexpectedly) caused separa- 
tion on the body and contributed to tunnel blocking (in 
conjunction with similar separation on the floor and 
sidewalls), see vortex formation at aft stations of body 
and floor; (5) the juncture between the floor, covered 
by relatively thick boundary layer, and the swept- 
forward support strut exhibited undesirable charac- 
teristics (probably a A shock on floor and a detached 
shock on wedge near the floor) judging by the bound- 
ary-layer outflow near wedge and the local vortex 
bubble at juncture of floor with sidewall. 

The preceding example is intended to draw attention 
to the known but often neglected truism in wind-tunnel 
testing that qualitative techniques may lead to better 
understanding of the flow than quantitative force and 
pressure measurements and may save much time (and 
money) especially in exploratory investigations. The 
strong boundary-layer effects due to oblique pressure 
gradients demonstrated in this report (and the fact that 
the resolution of probes and pressure orifices is de- 
creased because of the smaller scale of models and flow 
fields) suggest a more frequent coupling of visualiza- 
tion techniques with the standard pressure and force 
measurements in supersonic wind-tunnel testing. 


CONCLUDING REMARKS 


From reference 13 and the present study it is clear 
that to the previously recognized basic parameters 
governing the reflection of shock waves from solid 
boundaries (boundary-layer character and thickness; 
shock strength and nature-step or impulse) should be 
added the factor of the geometry of the three-dimen- 
sional shock-body configuration. Three complicating 
factors arise for truly three-dimensional shock-wave re- 
flections. First, in the cases considered, the exten- 
sion of the nonviscous singularity type of theory leads 
to an impasse: either the nonviscous reflection is 
highly anomalous and even impossible in some regions 
of the reflecting body, or the body in some manner 
warps the oncoming plane shock wave toward the up- 
stream direction. Second, a complicated post re- 
flection flow occurs even in the nonviscous case, the 
initial values for which should come from the non- 
viscous singularity theory. In the two-dimensional 





Fic. 15a. Shock reflection pattern from square body. 





Fic. 15b. Surface stream lines on side and top of square body 





le: 


Fic. 16. Surface stream lines on body wedge and tunnel floor 
for interaction of expansion wave and circular body. 


case gratifying progress in understanding the reflection 
of shock waves from laminar boundary layer on a wall 
has been achieved by studying simplified models of the 
viscous flow which approach asymptotically the non- 
viscous solution. In the three-dimensional cases stud- 
ied, even the nonviscous solution appears to be elusive, 
if it exists. Third, pressure gradients across stream 
lines due to the shock generate a boundary-layer cross 
flow that further modifies the effective boundary sur- 
face. 

These three factors appear to give rise to a general 
interaction pattern (flow and pressure fields) between 
plane shocks and elongated bodies with widely different 
cross sections. Its main feature is a strong boundary- 
layer cross flow and accumulation (separation) on the 
lee side roughly proportional to the shock strength (yaw 
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velocity component) and a rapid relieving of pressures 
on the incident side. The flow field gives indications 
of settling down to an asymptotic flow at shorter down- 
stream distances from the shock plane than expected 
from nonviscous considerations. In contrast with the 
two-dimensional case, this apparent asymptotic flow is 
intrinsically viscous, since the lee-side separation basi- 
cally modifies the effective boundary of the body. 
Thus, even if a quantitative understanding of two-di- 
mensional shock only a 
qualitative appreciation of the three-dimensional case 


reflection were achieved, 
based on experiments is likely. 

Inasmuch as the boundary-layer cross flow generated 
by oblique shock waves is likely to be important in 
many practical applications, it may be worth while to 
conjecture on how the findings of this report are likely 
to be affected by size. Increasing the size of a body or 
wing will generally not alter the strength or inclination 
of the shock relative to the oncoming stream lines, but 
the boundary layer will thicken. In accordance with 
the present results (Fig. 3a), the thicker boundary layer 
should accommodate larger pressure gradients normal 
to the surface so that the pressure gradients across 
stream lines near the wall should be decreased. On the 
other hand, the actual mass of low-velocity fluid in the 
inner portion of the boundary layer which responds to 
the pressure gradient will be increased by the thicken- 
ing of the boundary layer. It is therefore unlikely 
that the cross-flow effects will be less important on full- 
scale prototypes than on wind-tunnel models. 


REFERENCES 


1 Lees, L., Interaction Between the Laminar Boundary Layer 
Over a Plane Surface and an Incident Oblique Shock Wave, Prince- 
ton University, Aero. Engrg. Lab., Report No. 143, January, 
1949. 

2? Tsien, H. S., and Finston, M., Interaction Between Parallel 
Streams of Subsonic and Supersonic Velocities, Journal of the 
Aeronautical Sciences, Vol. 16, No. 9, pp. 515-528, September, 
1949. 

3 Robinson, A., Wace Reflection Near a Wall, Cranfield College 
of Aeronautics, Report No. 37, May, 1950. 

* Lighthill, M. J., Reflection at a Laminar Boundary Layer of 
Weak Steady Disturbance to a Supersonic Stream Neglecting 
Viscosity and Heat Conduction, Quart. Jour. Mech. & Appl. Math. 
Vol. II, Pt. 3, pp. 303-325, September, 1950. 

5 Fage, A., and Sargent, R. F., Shock Wave and Boundary 
Layer Phenomena Near a Flat Surface, Proceedings Royal Society 
A., Vol. 190, 1947. 


THE AERONAUTICAL SCIENCES 


APRIL, 1962 


6 Ashkenas, H., Dhawan, S., Liepmann, H. W., Marble F 
and Roshko, A., Problems in Shock Reflection, GALCI1 Report ‘ 
A.A.F., 1948. 

7 Shapiro, A. H., W., Some 
Experiments on the Interaction of Shock Waves with Boundary 
Layers on a Flat Plate, M.1.T. Gas Turbine Laboratory, Meteor 
Report No. 29, October, 1948; also, Jour. Appl. Mech., Vol, 17 
No. 2, pp. 126-131, June, 1950. 

8 Liepmann, H. W., 
pp. 39-66, Symposium on Experimental Compressible Floy 
U.S. Naval Ordnance Lab. Memorandum, No. 1133, June, 1949 

® Liepmann, H. W., Roshko, A., and Dhawan, S., On the 
Reflection of Shock Waves from Boundary Layers, GALCIT 
Report, 1950; N.A.C.A. T.N. No. 2334, April, 1951. 

Barry, F. W., Shapirv, A. H., and Neumann, E. P., The 
Interaction of Shock Waves with Boundary Layers on a Flat Sur. 
face, M.I1.T. Gas Turbine Laboratory, Meteor Report No. 52, 
March, 1950; also Journal of the Aeronautical Sciences, Vol. 18 
No. 4, pp. 229-238, April, 1951. 

1 Bardsley, O., and Mair, W. A., The Interaction between an 
Oblique Shock-wave and a Turbulent Boundary-layer, Phil. Mag,, 
Ser. 7., Vol. 42, pp. 29-36, January, 1951. 

12 Courant, R., and Friedrichs, K. O., Supersonic Flow and 
Shock Waves, pp. 318-350; Interscience Publishers, New York, 
1948. 

18 Migotsky, E., and Morkovin, M. V., Three-Dimensional 
Shock-Wave Reflection, University of Michigan Supersonic Wind 
Tunnel Memorandum No. 118, August, 1949; also, Journal of the 
Aeronautical Sciences, Vol. 18, No. 7, p. 484, July, 1951. 

14 Garby, L. C., and Nelson, W. C., University of Michigan 
8 X 13 inch Intermittent-flow Supersonic Wind Tunnel, Univer- 
sity of Michigan Memorandum No. 59, Engineering Research In- 
stitute, June, 1950. 


Neumann, E. P., and Barry, F 


Boundary-Layer Shock-Wave Interaction 


' 


% Culbertson, P., Calibration Report on the University of Michi- 
gan Supersonic Wind Tunnel, University of Michigan Memo- 
randum No. 36, Engineering Research Institute. 

16 Morkovin, M. V., Migotsky, E., and Bailey, H. E., On 
Three-Dimensional Shock-Wave Reflection, University of Michi- 
gan Supersonic Wind Tunnel Memorandum No. 191, January, 
1951. 

17 Migotsky, E., On the Reflection of Shock Waves in Three 
Dimensions, Ph.D. Thesis, University of Michigan, April, 1961. 

18 Murphy, J.S., and Phinney, R. E., Visualization of Boundary- 
Layer Flow, Readers’ Forum, Journal of the Aeronautical Sci- 
ences, Vol. 18, No. 11, p. 771, November, 1951. 

19 Dailey, C. L., and Wood, F. C., Computation Curves for 
Compressible Fluid Problems; John Wiley & Sons, Inc., New 
York, 1949. 

20 Jones, R. T., Effects of Sweepback on Boundary Layer and 
Separation, N.A.C.A. T.N. No. 1402, July, 1947. 

21 Allen, H. J., Pressure Distribution and Some Effects of Vis- 
cosity on Slender Inclined Bodies of Revolution, N.A.C A. T.N. 
No. 2044, March, 1950. 








rf) 








irble, F 
Xeport to 


V., Some 
Soundary 
’ Meteor 
Vol 17, 


‘eraction 
le Flow 
le, 1949 

On the 
‘ALCIT 


P. 1% 
lat Sur. 
No. 52, 


Vol. 18, 


veen an 
. Mag, 


mw and 


York, 


1sional 
Wind 
of the 


chigan 
niver- 
ch In- 


Tichi- 


[emo- 


.» On 
lichi- 
uary, 
lhree 
1951. 
lary- 


Sci- 


for 


New 
and 


Vis- 





Shock-Tube Studies of Transonic Flow over 


Wedge Profiles’ 


WAYLAND GRIFFITH? 


Princeton Unwersity 


SUMMARY 


Steady flow fields around two-dimensional wedges have been 
studied for Mach Numbers between 0.85 and 1.80 using the 
shock tube. It has been found that the location and shape of 
detached shock waves are independent of the wedge angle as 


jong as the shock is not near to attachment. 


The position of the sonic line and the pressure distribution over 
the wedge have been calculated from the density fields measured 
with a Mach-Zehnder interferometer. For a given wedge, identi- 
cal pressure distributions are found over a considerable range of 
supersonic flows. These results suggest simple formulas giving 
the drag coefficient for each wedge angle which are expected to 
hold in most of the transonic speed range. 

Experimental data on the drag of a thin wedge through the 
entire transonic range are available for the first time and are in 
excellent agreement with recent theoretical calculations based on 


transonic small-disturbance theory. 


SYMBOLS 


a = velocity of sound 

c = wedge chord 

Ca = drag coefficient based on D and ¢ 

Ci(t/2) = drag coefficient based on D and t/2 

Cp = pressure coefficient = (p — ~i)/(1/2)pi Vi? 

D = pressure drag of half wedge 

K = transonic similarity parameter = (1 VW) 
[ay +1)]*” 

L = slant length of wedge 

M = Mach Number 

P = pressure 

§ = detached shock distance measured from shoulder 
of wedge 

t = wedge thickness 

u’, v’ = perturbation velocities in x and y directious, re- 
spectively 

Vi = free-stream velocity 

X = fractional distance along incline of wedge from tip 
to shoulder 

¥ = specific heat ratio 

6 = wedge half angle 

p = density 

Superscript 
Oi = critical condition, state where local flow velocity = 


local velocity of sound 
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Subscript 
.” = stagnation value 
( ) = free-stream value 
( )e = refers to region behind detached bow shock 


(1) INTRODUCTION 


— PROBLEM OF determining flow patterns in the 
transonic range presents considerably greater diffi- 
culty than cases when the flow is purely subsonic or 
supersonic. From the theoretical standpoint special 
troubles arise first because the differential equations of 
flow are nonlinear and second because the regions of 
subsonic flow where the equations are elliptic are sepa- 
rated by an unknown boundary from supersonic regions 
having equations of hyperbolic type. Wind-tunnel data 
on flows near the speed of sound have so far been unob- 
tainable because tunnel wall interference results in 
blocking of the flow. 

This paper will describe a method by which the shock 
tube has been used as a sort of intermittent wind tunnel 
in which flow patterns over the entire transonic range 
are determined using a Mach-Zehnder interferometer. 
The results may then be compared with theoretical 
predictions and wind-tunnel measurements where these 
are available. 

Two rather distinct methods of approaching prob- 
lems involving mixed flows have been followed: one that 
deals primarily with the shape and location of detached 
shock waves in supersonic streams and the other that 
makes use of the extremely powerful transonic small 
disturbance theory! for treating flows in which the ve- 
locity is never far from the speed of sound. 

Maccoll and Codd? and Drebinger* have carried out 
numerical solutions for the flow around two-dimensional 
wedges at particular Mach Numbers, while some ex- 
perimental data on detached shocks have been ob- 
tained in a wind tunnel by Pack and Groth‘ and with 
the hydraulic analogy by Laitone.® Approximate 
methods for predicting the position of detached shocks 
have been developed by Laitone and Pardee® and 
Moeckel’ and for axially symmetric flows by Dugundji® 
and Laitone.? Experiments on axially symmetric 
flows have been carried out by Ladenburg, van Voor- 
his, and Winkler" and Heberle, Wood, and Gooderum.!! 
Lin and Rubinov'? and Busemann!* have discussed 
some general features of the problem. 

The transonic similarity approximation gives an 
equation that is linear in the hodograph plane, From 
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Fic. 1. The location of the shock reflected from a 90° wedge 
as a function of time after impact showing how a steady position 
is reached. 


this, properties of flows near sonic velocity having theo- 
retical interest are discussed by Guderley.'' Tsien 
and Baron’ have expanded the oblique shock and 
Prandtl-Meyer relations in terms of the similarity pa- 
rameter to find lift and drag coefficients for flat plates 
and asymmetric wedge airfoils in pure supersonic flows. 
Vincenti and Wagoner,'® Cole,” and Guderley and 
Yoshihara" have obtained solutions giving the pressure 
distribution over thin wedges. Recent wind-tunnel 
experiments of Liepmann and Bryson” have extended 
experimental data up to \J = 0.860 and down to V/ = 


mire ® 


(2) SUPERSONIC FLOW IN THE SHOCK TUBE 


Since the theory and operation of the shock tube have 
been well described already,” we shall include only a 
few special considerations relevant to the range of 
strong shocks used in these experiments. The Prince- 
ton shock tube is 4 by 18 in. in cross section and 38 ft. 
long. The tube is divided into regions of high and low 
pressures by a partition that shatters quickly when it is 
punctured. A shock whose strength depends on the 
original pressure ratio then advances through the low- 
pressure region accelerating the air to a certain Mach 
Number. The region of air flowing along behind the 
shock moves at uniform velocity and has the special 
features that no boundary layers or turbulence exist at 
the start. Following this region of uniform flow comes 
the ‘‘cold front,’”’ the boundary originally occupied by 
the diaphragm that now separates the compressed gas 
from the gas originally in the high-pressure end which 
has expanded and cooled. To accelerate air to a Mach 
Number of 1, a shock wave having a pressure ratio of 


very nearly 5:1 must pass through it. To generate a 
shock of this strength, an initial pressure ratio of 50: 1 
is needed. If helium is substituted as the high-pressure 
gas, the higher velocity of sound in helium makes pos- 
12: 1,21 


sible the more favorable ratio of Stronger 
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shocks will produce flows of Mach Number up to 1,89 
All of the wedges were '/, in. thick except for the 15° 
wedge which was 0.080 in. thick. As a result of the 
small size of the models and low density used in these 
experiments the Reynolds Number is small, the order 
of 10°. 

The fact that steady supersonic flow can be observed 
in the shock tube is shown by results similar to those 
given in Fig. 1. When a shock wave strikes an ob. 
stacle, a reflected wave travels out into the stream tel]. 
ing the oncoming particles of the objects presence, 
This wave grows for a time (about 500 microsec. for a 
model '/, in. thick) and then reaches a steady position, 
It has thus become the detached bow wave of the wedge 
in a supersonic stream. The steady flow is ended when 
(a) reflections from 
the tube walls return, (b) the incident shock comes 
back up from the end of the tube, (c) the rarefaction 


any one of several things occur: 


that started into the high-pressure gas when the dia- 
phragm burst reaches the bow wave, or (d) the ‘“‘cold” 
gas that was originally behind the diaphragm arrives, 
In the last case, a steady flow may again be set up before 
anything else happens, and no theoretical limit is placed 
on the Mach Number attainable. The highest in this 
tube to date is 3.7. A special advantage of studies 
made immediately behind the incident shock lies in the 
complete absence of a boundary layer prior to the pas- 
sage of the shock. It is therefore possible to observe 
nonviscous aerodynamic effects at the start and then 
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The two calculated values of Maccoll are marked M, and Dre- 
binger’s result is marked D. The solid line is the experimental 
curve, and the dashed line is the assumed transonic approximation 
for the detachment Mach Number. 
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see what modifications take place as steady-state 
boundary layers are formed on the models. 

All of the experiments reported in this paper were 
made in the ‘‘hot’’ gas behind the incident shock. With 
the existing arrangement alternatives (b) or (c) of the 
vious paragraph were always observed to spoil the 


pre . f i 
The reappearance of the incident shock reflected 


flow. 
from the end of the tube was clearly visible, the catchup 
of the rarefaction being somewhat obscure. No mech- 
anism like the Mach reflection for shocks acts to 
straighten up a rarefaction wave, so it usually arrived 
with a considerable slant and changed the flow pattern 
unsymmetrically. In another investigation, Weimer’ 
has shown that the ‘‘cold front’’ has a thickness about 
twice the tube width, or 8 in., and is clearly observ- 
able. The only remaining means whereby the flow 
could be disturbed without detection is from top and 
bottom wall interference. Pictures taken with an 
artificial wall placed near the model show that this too 
would be unmistakable, so the region of possible steady 


flow is reasonably clearly known. 


(3) DETACHED SHOCKS 


First, a set of experiments were made to determine 
the location and shape of detached bow waves on wedges 
of 10°, 15°, 20°, 30°, 40°, 60°, 90°, and 180° total angle. 
The bow-wave distance from the shoulder of the wedge, 
measured in terms of wedge thickness, is plotted as a 
function of Mach Number in Fig. 2. As has been 
suggested by several workers,® * '! the bow-wave posi- 
tion is practically independent of the wedge angle as 
long as the bow wave is not near to attachment. Under 
these conditions the approximate method for predict- 
ing bow-wave locations’ is in fair agreement with ex- 
periment. Data from the Michigan shock tube taken 
by Geiger and Mautz*’ are found to be taken mainly at 
times before equilibrium is established. With increas- 
ing Mach Number, bow waves for blunter wedges suc- 
cessively leave the family and reach a constant distance 
when attachment occurs. These points are indicated 
on the graph. As the flow Mach Number approaches 
1, the time required for the reflected wave to move out 
to an equilibrium position increases rapidly so that no 
data were obtained in this investigation on bow-wave 
positions for 17 < 1.16. In a later section a method 
will be described whereby this difficulty is sidestepped. 
In view of the fact that the experimental curve passes 
nearly through the detachment points it seems reason- 
able to assume tentatively that the locus of attachment 
points gives a good approximation to'the bow-wave 
distance curve in the interval 1 < M < 1.2. From 
transonic similarity, the formula for this is 


s/t = 1.558/(M? — 1)” (1) 


Bow-wave shapes for Mach Numbers near 1.37 are 
shown for various wedges in Fig. 3. Those for wedge 


angles above 40° are the same within the accuracy of 











re are 


14° 20° 40° 60° 90° 180° 
Fic. 3. The bow shocks for various wedges at Mach Numbers 
near 1.37. When the shock is far from attachment the wave 
assumes a “‘standard”’ form. The symbols on the horizontal axis 
locate the tips of the various wedges according to the notation in 
Fig. 2. 


observation, the 20° case is seen to have a different 
shape even though the foremost point is the same, and 
and the bow shock is attached and straight for a 14.4° 
diamond. According to the shock polar theory, there 
is a narrow range of Mach Numbers in which a shock is 
attached but the flow behind it subsonic. A few pre- 
liminary experiments bear out Guderley’s conclusion 
that the shock is curved under these conditions.'4 A 
detailed survey of the attachment process was not 
made because viscous forces near the surface introduce 
additional unknown effects, but some conclusions may 
be drawn with the information available. For a Mach 
Number well below that at which the bow wave will 
become attached to a wedge, the shape matches the 
“standard” shape, Fig. 4. As the Mach Number in- 
creases, the shock curvature increases, but its center 
position still follows the standard formula. When the 
shock just attaches, it is curved, and the subsonic field 
between it and the shoulder is nonuniform. Accord- 
ing to theory, the Mach Number range for transition 
from a detached shock to an attached shock with uni- 
form supersonic flow along the incline is only about 
0.01, but experiments indicate that the interval of 
nonuniform flow is much greater than this because of 
the boundary layer. Hence, the decrease in drag upon 
shock attachment in a real flow is probably not nearly 
so abrupt as potential flow theory predicts. 
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The “‘standard”’ shape of detached shocks for various 
Mach Numbers. 


(4) MertTHODS OF MEASUREMENT 


In previous studies with the shock tube, two meth- 
ods have been available for computing the incident 
shock strength: measurement of its velocity and ob- 
servation of the density change across it with the inter- 
ferometer. For shocks weaker than 6:1, these differ 
by less than | per cent, and it has always been assumed 
that the flow Mach Number may be found from the 
computed density ratio. In the supersonic range, wave 
angles for attached shocks may be measured after 
steady flow has been established. Mach Numbers 
computed from these values are larger than those de- 
termined from measurement of the incident shock ve- 
locity, while the fringe shifts observed across steady- 
state bow waves lead to agreement with the wave angle 


measurements. A number of special experiments were 





Fic. 5. 
a steady flow picture reveal the lines of constant fringe shift and 
therefore lines of constant density as distinct gray bands. Since 
the fringes are exactly matched ahead of the bow wave, the lines 
observed represent half-integral fringe shifts. The upstream in- 
fluence of the shock-boundary-layer interaction on the observa- 
tion windows is clearly visible. 


Superposition of the undisturbed fringe pattern and 


conducted to determine wave angles and shock attach. 
ment points for wedges of 10°, 14.4°, 20°, and 39° 
(M = 1.24, 1.32, 1.43, and 1.60, respectively). When 
the Mach Number so computed was plotted against the 
Mach Number determined from the shock velocity, a 
straight line made a satisfactory fit, extrapolating be no 
correction of the velocity Mach Number at J/ = 1,07 
and a correction of +10.5 percent at IJ = 1.6, the upper 
limit in these experiments. 

No tests have been devised as yet to enable one to 
choose between various conceivable explanations for 
this shift in flow Mach Number, but careful consider. 
ation points to the most likely cause as a change in 
specific heat ratio y at the high temperatures produced 
by these strong shocks or the effective narrowing of the 
channel due to wall friction. 

The theory of the Mach-Zehnder interferometer re- 
lates a density change pz — p; in two-dimensional flow 
to the fringe shift 7.2. by a proportionality constant that 
may be evaluated for any particular experimental ar- 
rangement. 

When an interferometer is adjusted to have equal 
optical paths everywhere, a uniform field is observed. 
If a disturbance is introduced, any fringes that appear 
will then be lines of constant density, but one will be 
unable to determine the order of a fringe and cannot 
compute the corresponding densities. For quantita- 
tive work, it is therefore necessary to observe the shift 
When a “‘blank’”’ (the undis- 
turbed fringe pattern) is superimposed on that showing 


of finely spaced fringes. 


the fringes in a disturbed condition, those regions where 
fringes are just 180° out of phase become readily visible 
as gray strips, and these are lines of constant density 
(see Fig. 5). By moving the blank normal to its 
fringes, the density contours may be followed over the 
entire field. The surprising feature of this method is 
the clarity of the gray region, its width appearing to be 
only about one-tenth of that between successive fringe 
crossings. Thus under good conditions, lines of con- 
stant shift may be located with an accuracy of about 
+ 1/20 fringe. 
as that using direct measurement but becomes less 


This method is, in general, as precise 


satisfactory in regions where fringes on the blank and 
the picture are nearly parallel or have widely different 
spacing. Numbers can be attached to the curves by 
taking some undisturbed region as the zero condition 
or by making absolute measurements of fringe location 
from marking pins placed in the field of view. 

When discontinuities in density are present, as is the 
case with shock waves, two special techniques are avail- 
able: setting the light source off axis and using the white 
With the first method, the effect of the 
discontinuity is distributed over adjacent rays so that 
If white 
light is used and the interferometer is in good adjust- 


light fringes. 
one can follow individual fringes across it. 
ment, only about three fringes will be clear, because 


different wave lengths quickly get out of step. The 
central fringe is easily located throughout the field and 
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ige shift computed in terms of the original fringe 


its frit 


spacing. 
(5) THe Sonic LINE 


The subsonic region immersed in supersonic flow is 
hounded on the rear by the sonic line, so the location of 
this line is of considerable theoretical interest. Since 
the density field is known, it is only necessary to com- 
pute the density ratio appropriate to sonic velocity. 
This will vary for each stream line, depending on the 
inclination of the shock wave. Two points may be 
found immediately. The air near the body has gone 
through a normal shock so its state is known exactly, 
and the point where the sonic line reaches the wedge 
Also the point where the sonic line 


may be computed. 
Since the density ratios 


meets the shock is easily found. 
for these points differ only slightly, the entire curve may 
be filled in fairly accurately. Representative examples 
are shown in Figs. 6-9. In cases where bow waves are 
the same, sonic lines are the same. 

Theoretical considerations regarding the intersection 
of the sonic line with the wedge due largely to Buse- 
mann!* and Guderley'* may be outlined as follows. In 
a supersonic region bordering subsonic flow downstream 
characteristics approaching the sonic line are shown to 
represent rarefactions, while the other family of char- 
acteristics represent compression waves reflected from 
the sonic line. If these are to be reflected from a 
boundary as expansion waves, again so that the flow 
i.e., without shocks—the boundary 
must be convex. A supersonic region embedded in 
subsonic potential flow therefore has to be adjacent to a 
solid body, since there is no other way a suitable bound- 
A second important conclusion may be 


remains smooth 


ary can exist. 
reached from an extension of the incompressible poten- 
tial flow theory result that the velocity at a convex 


corner must be infinite. In compressible flow then, 
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Fic. 6. Density field for wedge of 180° angle in supersonic 


flow. In Figs. 6-9, the fact that the sonic line reaches the body 
slightly ahead of the shoulder may be ascribed to the boundary 
layer. 
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Fic. 7. Density field for wedge of 60° angle in supersonic flow. 
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Fic. 8. Density field for wedge of 40° angle in supersonie flow. 
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Fic. 9. Density field for wedge of 20° angle in supersonic flow 


the velocity at a convex corner will certainly become 
supersonic. 

These two results—i.e., that flow over a convex corner 
has to be supersonic and that the entrance to a super- 
sonic region must be adjacent to a convex body—lead 
to the conclusion that the sonic line must be attached 
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Fic. 10. Density field for flow at WM = 1.50 past a plate with 
a half-circular nose. The sonic line meets the surface 59.5° around 
from the front of the body. 
to the shoulder of the wedge. Furthermore, it will be 
perpendicular to the incline, since the case of a sharp 
corner will just correspond to a Prandtl-Meyer expan- 
sion in a limited region. The figures show that in a 
real flow these considerations are slightly modified by 
viscous effects near the surface for the sonic line actually 
approaches the wedge a little ahead of the shoulder. 
It is important to note that the sonic line is not neces- 
sarily normal to stream lines everywhere. Consider- 
ing each stream tube as a case of one-dimensional flow, 
sonic velocity will clearly be reached at the minimum 
cross section, but there is no requirement that these be 
adjacent. 

When there is no sharp corner, the location of the 
sonic line is considerably harder to predict theoretically 
but may be computed from experimental results in the 
manner already described. For a body whose nose is a 
half-cylinder, the location of the sonic point will probably 
depend on M. For M = 1.50, it is located 59.5° 
around from the front as shown in Fig. 10. 

Maccoll and Codd? have computed numerically the 
density field for flow at M = 0.7 and 1.5 past a 40° 
wedge and at M = 1.5 for a square-ended body. Vor- 
ticity produced by the shock is shown to have little in- 
fluence and so is neglected. What appears to be a more 
serious assumption is that the sonic line is perpendicular 
to stream lines everywhere. The subsonic region pre- 
dicted is considerably smaller than observed in expert- 
ment but it will be seen in the next section how this has 
little effect on the accuracy of computed pressures on 
the wedge. 

A more accurate numerical solution, in the sense that 
no restriction was placed on the angle between the sonic 
line and stream lines, was worked out by Drebinger*® 
for a finite incline of 1:2 (53.2° wedge) at WV = 1.44. 
The resulting shock and sonic line are compared with 
experiment in Fig. 11. Except for the vicinity of the 
shoulder already discussed and the fact that the bow 
wave is slightly farther out in the computed case, the 
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two results are taken to be in agreement within the 


experimental uncertainty. 


(6) PRESSURE DISTRIBUTION 


Pressure distributions on the front surface of the 
wedges have been computed and are shown in Fig, |9 
for the examples presented in Fig. 3. Here the ratio 
of the pressure to free-stream pressure p//, is plotted 
as a function of the fractional distance XY from the tip 
to the shoulder. Since there is always zero velocity at 
the tip and (nearly) sonic velocity at the shoulder with 
a detached shock, the end points are the same for the 
same Mach Number. When the bow wave is attached, 
as for the 14.4° diamond at M = 1.39, the high stagna- 
tion pressure is not attained at the tip, and the pres- 
sure drag is much reduced. 

A comparison of the pressure distributions on the 
same wedge at different Mach Numbers may be made 
by plotting p/p,* where p* is the pressure correspond- 
ing to sonic velocity. According to theory / is equal 
to p* at the shoulder and [(1/2) (y + 1)]”"~"p* = 
1.8929p* at the tip. In Fig. 13, data for the 60° and 
180° wedges at several Mach Numbers have been 
plotted in this way. The points lie along the same 
curve within the accuracy of the measurements. The 
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Fic. 11. Comparison of experimental and theoretical shock 
and sonic line locations for flow at J = 1.44 past a 53.2° wedge. 
The dashed curves are Drebinger’s calculated result. 
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Fic. 12. Pressure distribution from the tip to the shoulder 
in terms of free-stream pressure for several wedges at Mach 
Numbers near 1.37 and a 14.4° diamond at M = 1.39. The 
shock is attached to the diamond. 


possibility of this was first suggested by Maccoll and 
Codd, who found that the pressure distribution curve 
for a 40° wedge had the same shape at M/ = 0.7 and 
1.5. Their result for a 180° wedge is indistinguishable 
from the experimental curve shown in Fig. 13, as is the 
case with the 40° wedge on a curve not shown. Since 
the pressure at the tip and shoulder are independent of 
the computed shock position, it is reasonable that the 
computed pressure distribution along the incline is 
relatively insensitive to variations in shock position as 
long as the shock is at least far enough away to have 
the standard shape. 

Using this result, a drag coefficient for the transonic 


range may be found as follows: 


b/Pi = (b/p*) (p*/P) (bx»/pr) = 
[2 (y + OT ‘Y= “(p p*) ( pe0/ Pi) 


The last term is the familiar Rayleigh pressure ratio, 
which is a function of Mach Number only. The pres- 


sure coefficient 


Cy = (p — pr)/(/2)ei Vi? = 2(p/pi — 1)/yM? (2) 
thus becomes 


4 2 inl i 2 po 

© * {( ) ; 1 s) 
yMPL\y + 1 p” pi 

The drag coefficient based on wedge chord as used 


by Cole, Liepmann, and Vincenti is defined as 


Ca = D/(1/2)piVi'e (4) 


OVER WEDGE 


bo 
wt 
or 


PROFILES 


D is the pressure drag of the half wedge per unit width. 
In order to obtain a meaningful value for the flat- 
nosed model (180° wedge) and to facilitate comparison 
with Maccoll’s computed results, a drag coefficient based 
on the half thickness is defined as 


af 
Ca(t/2) = D/(1/2)pi.V12(t/2) - C,dX (5) 


where X is the fractional distance along the incline from 
the tip to the shoulder. 
For the 60° wedge this gives 
C,(t/2) = 1.429(0.773 p/p. — 1)/ Me (6) 
and for the 180° wedge 


C,(t/2) = 1.429(0.917 p/p, — 1)/M2 (7) 


This is in close agreement with Maccoll’s result, 


C,(t/2) = 1429(0.9277 p/p: — 1)/My? 


From what has been said these formulas should yield 
the front pressure drag for Mach Numbers in the in- 
terval from about 0.7 to at least 1.7. 


(7) TRANSONIC SMALL-DISTURBANCE THEORY 


The transonic similarity theory is the counterpart of 
subsonic approximations leading to the Prandtl- 
Glauert formulas and Ackeret’s rule in supersonic flow. 
Here perturbations from the critical velocity a* are 
assumed to be small compared to a*, so that the x and 
y velocity components are 


, 


u=at*+u’, v=v', u'v' <<a* (8) 


If irrotational flow is assumed so that a potential func- 
tion & may be introduced for the perturbation veloci- 


ties the equations of gas dynamics reduce to 


O° eb, is. i. Of O° 2 OP O*b (9) 
Ox Ox? a* Oy Oxdy 


Oy" ge” 
By means of an affine transformation, the appropriate 


similarity parameter is found to be 


K = (1 — M)/[i(y + 1)]” = 
(1 — M?)/2[6(y + 1)]** (10) 
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pressure p* for the 60° and 180° wedges at several Mach Numbers 
showing that the distribution is independent of Mach Number. 
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where 6 is the wedge half angle. Each solution of the top and bottom walls of the tube. The incline Was sec 
equation therefore provides an infinite number of pat- followed by a 0.6-in. straight section and was held from nut 
terns for similar bodies at appropriate Mach Numbers. the back by a sting extending to the end plate on the 10 
In the hodograph plane, a linear equation is obtained shock tube. The model was mounted nearer the top hig! 
known as the Tricomi equation, of the tube than the bottom so that any wall interfer. sup 
(02y/Ou,2) — my, (02y/du,2) = 0 ence would cause an unsymmetric flow. No such effect exp 
uy = (y + 1)u’/a,* 1% = (y + 1)v’/a*f sad saertneepeieane ; mot 
lhe fact was mentioned earlier that for Mach Num. 
It is an equation of this general form for which several _ bers close to 1 sufficient time is not available for the re- 
solutions have been obtained. Guderley and Yoshi- flected bow wave to reach an equilibrium position 
hara'’’ have computed the flow over a wedge-shaped Obviously, in subsonic flow it would never do so but 
profile in a free stream at a Mach Number of | and would continue back up the tube indefinitely and be. we 
obtain for the drag coefficient come vanishingly weak. As it turns out experimentally, Ma 
Co* = 1.7588 /(-y + 1)” (12) this makes no difference to the establishment ofa steady nh 
pressure distribution on the wedge for after about 350, 8 
which for a 15° wedge in air gives Cp* = 0.0441. — sec. equilibrium is observed, which lasts for 1,600x sec., sm 
Cole” uses a doublet to represent the singularity in the — until the original shock returns from the end of the tube. _ 
hodograph plane corresponding to the flow at infinity The advancing bow wave becomes nearly normal and the 
and gets a series representation of the drag coefficient just replaces the compression extending to infinity by a - 
for 7 <1. AtM=1, sudden jump in pressure. Thus it is possible to obtain - 
Co* = 1.6787/(y + 1)” (13) steady flow patterns in the neighborhood of a model 
i even though conditions are not steady everywhere in 4 
in good agreement with Guderley and Yoshihara. The _ the tube. Bm 
slope of the drag curve at 1/ = 1 is found to be 46/ The drag coefficients on a 15° wedge computed from “a 
(y + 1). This has special theoretical significance, experiments over the range of times just mentioned 598 
which will be discussed later. for 0.86 < M < 1.11 are plotted in Fig. 14. Results i 
Vincenti and Wagoner'® solved the Tricomi equation from measurement of the flow past the 15° wedge and ma 
by numerical relaxation methods for four values of the a 14.4° diamond of 3-in. chord with steady bow waves ai, 
similarity parameter chosen to fill in as well as possible for 1.15 < MW < 1.40 are also plotted. The circles are cli 
the region between 1/ = 1 and the pure supersonic wind-tunnel data of Liepmann and Bryson." 5 
régime. In terms of A defined by Eq. (10) the values Within the approximations used in transonic theory, 7" 
used were —0.242, —0.352, —0.460, and —0.529. A the pressure coefficient may be written alternatively as -” 
discussion of the proper way to apply results of this ee hla atin val 
small disturbance theory to real airfoils for Mach Num- Cr = (b — pi/(1/2 pr | sca Maa on 
: ees Eh aa eee —2u’/V, = 2(M,? — m*)/(y + 1) (14 eat 
bers not equal to | is included in the next section. 
where m is the local Mach Number. If one uses the . 
(8) EXPERIMENTAL RESULTS FOR A 15° WEDGE AND last form and applies the conclusions of Section (6) ag 
COMPARISON WITH THEORY that the local Mach Number distribution is constant a 
The 15° wedge used for experiments in the sonic '™ nee tear 2 ee cl>indeenmmpee: Aane! OM) a, =1 is of 
range was made much smaller (¢ = 0.080 in.) than pre- Oand Ge Rope nt Che drag curve at 0 = Iie 
vious shock-tube models in order to increase the rela- dCp/d Al, = 46/(y + 1) (15) . 
tive flow duration and cut down the influence of : : : ae ; : - 
This was pointed out by Liepmann and is identical with for 
,  Cole’s theory. From Fig. 14 this is seen to be in ex- ins 
m cellent agreement with experiment. Fig. 15 shows th 
experimental pressure distributions for Mach Numbers to 
: PS ie near 1. The vertical marks denote points where J/ = of 
F : 1 and 0.815 (the Mach Number at mid-chord according Vi 
We = to Guderley’s calculation). 
4 .* Vincenti and Liepmann discuss thoroughly the diffi- 
Yi. culties that arise in applying small perturbation theory 
ie. .% to real flows. As they point out, there is no theoretical Ty 
; | means of choosing among the various forms for C, in 
; Eq. (14), sitice they are all equivalent in the approxt- - 
i F . r qo ation used. Similar uncertainties exist in the value 
M of the Mach Number which should be associated with U1 
een ie ee ae eas Rally Papers each value of AK for a given wedge. Accordingly, 
curve in terms of experiment and theory : two numbers are computed for C, corresponding to the Fi 
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second and third forms for C, in Eq. (14), and the Mach 
number 1s found from the first expression for A in Eq. 
(0). The results are plotted in Fig. 14. Clearly, the 
higher values in subsonic flow and the lower ones in 
supersonic flow correspond to taking —2u’/V;. The 
experimental data enable us to choose this form as the 
more reliable one to use. 


(9) CONCLUSION 


The flow around wedges in the transonic range has 
been studied experimentally using a shock tube and 
Mach-Zehnder interferometer. Detached shocks for 
wedges of 10°, 15°, 20°, 30°, 40°, 60°, 90°, and 180° 
total included angle have been observed for 1.16 < AJ < 
1.80. The shape and location of detached shocks, which 
are measured appropriately in terms of distances from 
the shoulder where the sonic line reaches the wedge, 
are found to be a function of Mach Number only and 
independent of the wedge angle as long as the shock is 
not near to attachment. The numerical solution for 
a 53.2° wedge at MJ = 1.44 obtained by Drebinger is in 
good agreement with experiment, while the shock and 
sonic line predicted for a 40° wedge by Maccoll and 
Codd are somewhat in error, presumably because their 
assumption that the sonic line crosses stream lines nor- 
mally everywhere is too severe. A measurement of 
experimental pressure distributions over the wedge in- 
cline shows that Maccoll’s values at the surface are cor- 
rect, however, and that identical distributions exist 
for a considerable range of supersonic flows. This re- 
sult and Maccoll’s observation that the computed 
values for a 40° wedge at AJ = 0.7 and 1.50 are the same 
suggest a simple formula for the drag coefficient for 
each wedge angle which is expected to hold in most of 
the transonic range. Such extends to 
wedges of any angle the predictions of transonic small 
disturbance theory that the pressure distribution on the 
inclined portion of a thin wedge is constant over a range 


a conclusion 


of Mach Numbers near |. 


For a 15 
tions may be observed in the neighborhood of the model 
for all Mach Numbers between 0.860 and 1.16, thus clos- 
ing for the first time the gap in experimental data around 
the speed of sound. Measured drag coefficients are found 
to match wind-tunnel data of Liepmann at either end 
of this interval and confirm the theoretical predictions of 


Vincenti, Cole, and Guderley. 


wedge it has been found that steady condi- 
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Remarks on Dynamic Loads in Landing 


ADAM ZAHORSKI* 


North American Aviation, Inc. 


SUMMARY 


Engineering aspects of the solution of the equation of motion 
for the inertia forces on an aircraft wing during the impact of 
landing are presented. A new and simplified, yet rigorous, ap- 
proach to the analysis of dynamic loads which can be easily 
applied to an arbitrary forcing function is described. Working 
limits are established for the number of modes of vibration which 
are effective in contributing to the dynamic responses. Ex- 
amples are given to show how the inertia forces are affected by 
different modes of vibration, by variation in mass distribution, 
and by shapes of the forcing function. 


STATIC AND DyNAMIcC LOADS 


QO” of the more difficult problems of structural 
engineering is to determine the magnitude of the 
external loads acting on a structure. Most loads fixed 
by higher authorities or codes are usually specified as 
static loads, subdivided into dead and live loads. How- 
ever, live loads are basically static loads—i.e., for an- 
alysis purposes, live loads may be considered to be 
static loads that are varying relatively slowly in a pre- 
scribed sequence. Static loads as such are usually not 
difficult to evaluate and can be checked experimentally 
with relative ease. But, on the other hand, it has been 
known for a long time that structures must be able to 
sustain considerably higher loads than those calcu- 
lated from the principles of statics. In many cases 
these higher loads are caused by forces moving over the 
structures and producing dynamic effects. Dynamic 
forces are much more difficult to calculate, much harder 
to interpret physically, and, because they act for only 
a short time, are difficult to measure. Only recently 
have reliable instruments been developed which can 
record transient forces for engineering analysis. An 
airplane is subjected to dynamic loading by flight 
through gusts, by rapid maneuvers, by the firing of guns, 
and by the impact of landing. This paper considers 
the dynamic loads developed over a wing structure 
during and just after the impact of landing. 


Basic EQUATIONS 


Fundamentally, the problem is similar to the case of 
a bar that vibrates when struck with a sharp blow. 
Vibration implies that part of the vibrating body is 
undergoing accelerated motion and that inertia forces 
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are opposing the acceleration. These inertia forces are 
the dynamic forces that must be determined in the 
problem. 

There are two kinds of the inertia forces: (1) those 
that are associated with linear accelerations and are 
expressed by formulas of the type F = ma, and (2) 
those that are associated with angular accelerations and 
are expressed by 7’ = Ja. The former force usually 
causes the bending or flexure of a wing; the latter, the 
torsions or twisting of a wing. 

The problem is so broad in scope and so complex 
mathematically that it is necessary to establish certain 
reasonable limitations. The first limitation is that 
only the dynamic loads that produce bending of the 
wing will be considered. (However, many statements 
also could be extended to include torsional loads.) 
The second limitation is that the analysis is confined 
to landing in which the main gear wheels touch the 
ground simultaneously, thereby producing identical 
force vs. time relationships at each wheel and only 
symmetrical modes of vibration in the wing. 

The mathematical tools are well developed for han- 
dling the response of an elastic system to transient 
forces. Physicists have derived mathematical equa- 
tions that represent the natural laws; mathematicians 
have developed many ways of solving these equations. 
The engineering task is threefold: (1) to choose mathe- 
matical constants that appear in the equations so that 
they realistically represent the physical system; (2) 
to select and to simplify the mathematical methods for 
the problem in hand so that a solution can be obtained 
in the shortest possible time with a sufficient degree of 
accuracy; and (3) to present the final results in such a 
manner that additional facts may be deduced concern- 
ing the laws of nature. 

The wing is represented mathematically by a line 
possessing certain elastic properties, known as EJ. For 
each element of the line an equation of motion is set up 


O° (: 7 a) n o*y fae, £) (1 
=i m = f(x, 
Ox? Ox? of? = 


This partial fourth-order differential equation repre- 
sents the summation of forces in the vertical plane: 


as follows: 


elastic forces of deformation, inertia forces, and exter- 
nally applied forces. The equation is similar to the 
one used by electrical engineers where they desire to 
find the response of a circuit to suddenly applied voltage. 
We are obliged to them for various methods of solution 
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Fic. 1. Actual and approximated drop curves showing the forcing 
functions for maximum and normal landing weights. 


of those types of equations known as “boundary prob- 
lems.”’ 

The method of solution adopted in this paper is that 
of the natural modes of vibration, sometimes called the 
“method of characteristic or eigen-functions.”’ 


EXTERNAL TRANSIENT FORCES 


External transient forces, which cause vibration in the 
wing, are represented in Eq. (1) by f(x, ¢) and are as- 
Actually, they are not known and 
Sev- 


sumed to be known. 
thus must be determined as a separate problem. 
eral attempts have been made to determine these forces 
analytically, showing the variation of force vs. time 
graphically. Some methods proposed are too difficult 
to be of practical value, since they involve the solution 
of nonlinear differential equations for any particular 
case. Others, such as those proposed by Wright Air 
Development Center,’ are simple in calculation but are 
based on arbitrary assumptions and seldom agree with 
the experimental results. However, one W.A.D.C. 
assumption can be adopted without serious error—.e., 
the trapezoidal shape of the forcing function—which 
considerably simplifies the lengthy calculations for 
dynamic response. 

Many dynamic analyses have been made on the 
basis of forcing functions that were obtained from drop 
testing of the whole gear assembly. Results of such 
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analyses may be extremely misleading. The technique 
of drop testing was developed primarily for testing the 
landing gear at ultimate load and not for obtaining dy- 
namic load data for design of the wing. Because the 
dynamic response, both in magnitude and shape, is not 
linearly proportional to the maximum of the applied 
load, information obtained from an arbitrary free drop 
(e.g., 27 in.) of the landing gear is practically useless so 
far as the dynamic loads are concerned. 

The four drop curves shown in Figs. la and 1b were 
obtained experimentally for the same landing gear at 
two different weights and at two different drop heights. 
They are supposed to represent conditions of maximum 
weight and normal weight landings at limit and ultimate 
landing-gear load factors. Dynamic responses are cal- 
culated for each condition, thus showing the effects of 
variation of two parameters. For the purpose of sim- 
plifying computations, each drop curve was approxi- 
mated by a trapezoid or a triangle, as shown by the 


dotted lines. 


SOLUTION OF EQUATION OF MOTION 


The solution of equation of motion [Eq. (1)] by the 
method of characteristic values assumes that the de- 
flection y(x, ¢) can be represented by a sum of two func- 


tions, such as 
y(x, t) = > Vj(x)p,(t) (2) 
j=0 


where J(x) is the function of the span coordinate (x) 
only and is known as the shape of the jth natural mode 
of vibration and where p(t) is the generalized time 
coordinate associated with the jth mode. 
It can be shown? that for all natural modes of vibra- 
tion the following relations are valid: 
tl. Ff; 


4 


dx? dx? 


b . 
ry 7 _ JOfor] # k 
f mY,;Y,dx = \N, for j = & (4) 


where w, is the frequency of the jth natural mode and 


= wm; (3) 


N;isa constant. 

It will suffice here to state that Eq. (3) can be solved 
for Y,, the shape of natural modes, and that, without 
loss of generality, Y; at the wing tip can be set equal to 
unity. In that case, p,(t) in Eq. (2) can be interpreted 
also as the tip deflection of the wing caused by the ex- 
ternal disturbing forces f(x, ¢). 

If Eq. (2) is substituted into Eq. (1) and modes of vi- 
bration Y, eliminated by Eqs. (3) and (4) then the dif- 
ferential equation of motion Eq. (1) will be transformed 


into 
P(t) + w;*pj(t) = A;F(t) (5) 
where 


A, = RY,(L)/N; (6) 








260 JOURNAL OF THE 

R = the maximum vertical force exerted by the 
landing gear on the wing 

Y,(L) = the ordinate of jth mode at landing gear 

N; = defined by Eq. (4) 

F(t) = the forcing function, a 
quantity that describes the variation of 


nondimensional 


transient forces applied to the wing 


The maximum value of the forcing function is taken as 
unity. Thus, the product RF(t) is the history of the 
transient. 

The solution of Eq. (5) is known in the form# 


p(t) = (Aj/w))D,(t) 


“J 
— 


where D(¢) is Duhamel’s integral 
t 
Dt) = [ F(r) sin wj(t — 7) dr (8) 
0 


Integration of Eq. (8) is usually difficult. If F(7) is 
an experimental curve and is given by a set of ordinates, 
then the numerical integration becomes tedious and 
must be done by a modern high-speed computing ma- 
chine. Even when high-speed machines are available, 
a prohibitive length of time is usually required to de- 
termine the dynamic response of the wing for many 
different forcing functions that may arise for various 
landing conditions. 

It is possible to simplify calculations if we note (1) 
that for the dynamic analysis the chief interest is in the 
distribution of accelerations along the wing and not in 
the distribution of deflections and (2) that certain types 
of superposition are permissible in calculation of Duha- 
mel’s integrals and one simple integration will be suffi- 
cient to cover vast varieties of forcing functions. 

First, note from Eq. (2) that the acceleration at any 
point along the wing is 


Hx, 1) = YS V(x) AO (9) 
j=0 
and from Eqs. (5) and (7) that the wing tip acceleration 
is 


p(t) = A;[F(O) — w,D;(0)] (10) 


GRADIENT FORCING FUNCTIONS 


Let F,(t) be a forcing function that starts at 4; = 0 
and uniformly increases with time. This function 
can be represented by a straight line passing through 


the origin. Since the slope is constant, we have 


F(t) = F,{(t)-t; (11) 


where F;,'/(/) is a constant and subscript 7 indicates 
that time ¢; is measured at the origin of F;. The func- 
tion represented by Eq. (11) will be called a “gradient 
forcing function.’’ Resolution of landing force vs. 
time history into four components, each of which is a 
linear function of time, was proposed by Midwest Re- 
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Trapezoidal forcing function (a) and its equivalent in 
straight line functions (b). 


Fic. 2. 


Substituting Eq. (11) into Eq. (8) and integrating, 
we get Duhamel’s integral for the gradient functions 


Dij(t) = Fi ((ti/w;) — (1/a,;?) sin wyt;] (12) 
and, substituting the above into Eq. (10), we get 
bib) sad A,(F,’ /«;) sin wy; (13) 


In general, the forcing functions of any shape F(!) 
can be approximated by a number of straight lines form- 
ing a polygon, and the integration of Duhamel’s in- 
tegral may be performed along several broken lines. 
However, this integration can be greatly simplified by 
superposition of Eq. (13), which will be demonstrated 
on a trapezoidal forcing function. 

Let F(t) be a trapezoidal forcing function whose 
T;, T2, and 


This function is replaced by one, two, 


shape is defined by three time intervals: 
T; (Fig. 2a). 
three, or four straight-line functions Fo, Fi, F2, and F;3 
(Fig. 2b), with slope so adjusted that the summation 


of all forces at any time ¢ is equal to F(t). Thus, 


F(t) = >F; (t) 

Because the superposition of forces is permissible in 
calculations of responses (reference 4, page 405) and 
if the time coordinates are measured from the proper 
origins at intersections of Fy, Fi, Fs, and F3 with the 
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abscissa aS shown on Fig. 2b, then 
D,(t) = D(to) + Dj) + ete. = DED ji(t,) 
and, consequently, from Eq. (7), 
p(t) = pte) + p(t) + ete. = Depylti) 
and 
p(t) = Pilto) + bt) + ete. = Wpjilts) 
where Pj;(t;) is given by Eq. (13). 
For example, by superposition of response, we can 
say that the acceleration in the region between 7) and 


T; (Fig. 2a) is 


] . l ; 
pt) = A; (=) Sin wjlg — () Sin wf; 


or, since? = fg and?, = ¢ — 7), 


l l 
A ( ) sin wt — ( --) sin w(t — 7) | 
wT wl) 


This can be written 


} = 
p(t) 


sin (1/2)w,;T; 


(1/2)w,T, cos w,; [t — (1/2)7T;] (14) 


L 


p(t) = A; 
This equation is of a particular interest. The form 
of the equation is very similar to that derived by 
Wasserman.’ However, it is much more general and 
reveals that both the lag angle and the amplitude of the 
acceleration are functions of the frequency of natural 
vibration and of the slope of forcing function. In 
fact, the amplitude of ,(¢) is an oscillatory function of 
the product w;7}. If (1/2)w;71 = a, 27, 32, etc., the 
amplitude of acceleration will be zero and therefore 
there will be no dynamic response in that particular 
mode in this region. 

Fig. 3 shows the results of calculation of the dynamic 
response in one mode only due to three different forcing 
functions that vary only in the initial rate (slope) of 
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Fic. 3 
of forcing function. 
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the application of the disturbing force. For computa- 
tion of these curves, Eq. (14) was modified somewhat, 
as it will be shown later, by introduction of the struc- 
tural damping that changed a simple harmonic motion 
to a motion with an exponential decay of amplitude 
after the break in the forcing function. The dynamic 
magnification factor in that figure is defined as a num- 
ber by which the design load factor must be multiplied 
to obtain vertical acceleration at the wing tip in units 
of gravity due to dynamic response of the elastic mode 
only. However, a precise definition is of no importance 
here, since the curves are presented only as an illustra- 
tion that the initial slope of the forcing functions may 
be an important parameter in the calculation of dy- 
namic responses and that an arbitrary drawn line may 
give misleading results. The figure shows that, if the 
force attains its maximum in 0.04 sec., then the mag- 
nification factor for wing tip acceleration will be 2.12 
but that, if the same force attains the same maximum 
in 0.07 sec., then the magnification factor will be only 
0.34. 


NUMBER OF EFFECTIVE MODES 


How many natural modes of vibration should be 
used in the dynamic analysis is the question of practical 
importance. Too few or too many modes may lead 
to erroneous conclusions. 

Two different criteria may be simultaneously used 
in order to limit the number of effective modes. The 
first criterion is based on the consideration of the mag- 
nitude of wing tip deflections of various modes in rela- 
tion to overall dimensions of the wing; the second con- 
siders the effects of the structural damping on the 
amplitudes of deflection and of acceleration. 


Trp DEFLECTIONS 


The amplitude of the wing tip acceleration /,() 
will be, in general, proportional to the ratio V;(L)/Njw;. 
In this ratio, only w; increases for higher modes, but 
there is no prescribed pattern of values for both Y,(L) 
and N;. Therefore, it might be expected that high 
acclerations may occur at some higher frequencies 
than were considered in the analysis. On the other 
hand, it can be seen from Eq. (5) that the ratio p,/p, 
is in all cases inversely proportional to w;*._ This means 
that deflection amplitudes will in all cases diminish 
rapidly for higher modes in contrast to the case for ac- 
celerations. 

To arrive at a rational conclusion as to how many 
modes are necessary in an analysis, the order of magni- 
tude of the wing tip deflection of each consecutive mode 
must be considered in relation to the overall dimensions 
of the wing. For the wing of large size, the tip de- 
flection of a few hundredths of an inch, which was 
theoretically derived for a line beam with uniformly 
distributed masses and stiffnesses, is not realistic. It 
can never be detected experimentally and, in fact, 
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does not exist on the actual wing. If the computed de- 
flections are not realistic, then it is reasonable to deduce 
that the accelerations derived from those deflections 
also are not realistic. 

Figs. 4a and 4b show wing tip deflections calculated 
for the first three elastic modes at two different load 
conditions. These deflections demonstrate that for a 
particular airplane the modes of beyond the third 


should not be considered in the dynamic analysis. 


STRUCTURAL DAMPING 


The second approach to limiting the number of modes 
is the consideration of the effects of structural damping 
upon deflection and acceleration amplitudes. 

There are some questions as to what the properties of 
structural damping are, how damping should be used 
in the differential equations, and if the mode of vibra- 
tion still remains orthogonal. But, if the damping 
coefficient is small, say of the order of 5 per cent, then 
there should be no objection to using basically the same 
method as in the case when damping is neglected. 

If it is assumed that the modes of vibration are 
damped out proportionally to exp(—yot), then a simple 
diagram showing the variation of amplitude with differ- 
ent frequencies will be helpful in deciding how many 
modes should be included in the analysis. 

If damping is introduced into the equation of motion, 


then Duhamel’s integral will be, instead of Eq. (8), 
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D,(t) = / Firje peo ( 
/0 


This will be integrated for a damping coefficient 
“« = 0.05 (assuming that u << 1) and a gradient forcing 
function F(r) = F’(r)r, where F’(r) is a constant. 


” sinw(t — r) dr (15) 


Dt) = (wt — 0.10) — e~°™ “i sin (wt — 0.10) 


(16) 
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Substituting the above into Eq. (10), we get 


pit) = Aj(Fi'/w;)dy (17) 
wihiere 


di; = 0.10 + e~° *i sin (wt; — 0.10) (18) 


A new function d; is introducéd because of its use- 
fulness in the calculation of the dynamic response and, 
in addition, because it can be generated quickly (about 
35 per min.) on IBM equipment. The dynamic re- 
sponse to any arbitrary forcing function which can be 
approximated by a number of straight lines need be 
obtained by the generating function d; only once for any 
given frequency w. At this stage, proper application 
of superposition means only automatic counting and 
stacking the cards by the machines. 

When ),(t) has been determined by the 
method, all that remains is to substitute its value into 
Eq. (9) to obtain the spanwise distribution of acceler- 
ation in the form of a family of curves for various values 


above 


of time. 


COMPARATIVE RESULTS 


The method outlined above for the solution of the 
dynamic loads is simple in application and may be 
handled by relatively inexperienced engineers. 

Thus, it is practicable to calculate the wing tip ac- 
celerations (/;) due to four different forcing functions 
(Fig. 1) from ¢ = 0 to ¢ = 0.500 sec. after the impact 
for each 0.001-sec. interval (Figs. 5 to 8). 

It is entirely possible that in many engineering offices 
only one dynamic analysis is made on the bases of only 
one drop curve and one weight condition. Any one of the 
diagrams of Figs. 5 to 8 could be accepted as typical 
for what could be expected in the behavior of the wing 
during the impact of landing. These four diagrams 
show that the variations in p,(¢) are very large, and any 
conclusions derived from any one diagram may be er- 
roneous if projected on any other landing conditions. 
















006 
60 t 
40 
xz 0 05 SECONDS 
r ; : 
4 0 04 SEGONDS 
= 30 
z 008 SECONDS 
° 
4 , 
$ 20 O11 SECONDS 
. 
z 
o 
a 
z 10 
° 
E o —<$$_+_—_—> NY 
1o 
CENTER WING 
LINE TIF 
<i j 
oO 10 20 30 a 50 60 70 8 30 100 
PERCENT OF SPAN 
Fic. 12. Dynamic magnification factor in bending vs. span 


for various times for normal landing weight (ultimate drop curve) 
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This possibility of error can be further emphasized by 
showing, through the application of Eq. (9), the span- 
wise distribution of the dynamic load factors (Figs. 9 
to 12). 

Load factor or, in this case, the dynamic magnifica- 
tion factor (d.m.f.) is defined as the acceleration in 
gravity units per unit landing load factor at the gear. 
This permits comparison of all curves on the basis of 
one common denominator that is the unit of the maxi- 
mum reaction at the wheel during landing and permits 
extrapolation of the effects on dynamic loading due to 
It is 
the author's contention that inasmuch as accelerations 


small changes in weights and intensity of landing. 


due to dynamic responses are not linear functions of 
landing load factors, it is much more realistic to use for 


all design purposes only those forcing functions that 


actually can occur in landing and to apply the desired 
safety factor to the dynamic accelerations only afte, 
the calculations are computed. 
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A Navy Symposium on Aeroballistics, under the sponsorship of the Bureau of Ordnance, will 
be held at the Huntington Hotel in Pasadena, Calif., on May 13 and 14, 1952. 
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A Study of Transonic Gas Dynamics 
by the Hydraulic Analogy 


E. V. LAIFONE* 
Unwersity of California at Berkeley 


SUMMARY 


The theory governing the propagation of surface water waves 
js analyzed to show that an analogy with the two-dimensional 
flow of a perfect gas exists only if the water depth is approxi- 
mately '/, in., the model is fairly large, and the shock waves or 
hydraulic jumps present are of the weak type having a negligible 
increase in entropy. 

Photographs of the water surface waves corresponding to 
shock waves in front of simple wedges and circular cylinders are 
presented to illustrate these restrictions. It is pointed out that 
the surface-wave group velocity, which could be considerably 
less than the so-called wave velocity, provided the only correct 
analogy to the speed of sound in two-dimensional gas flow over 
aclosed body. 

Motion pictures were taken in order to study the effect of con 
stant acceleration on bow shock waves. It was found that at 
high rates of acceleration the location of a detached shock wave 
was dependent mainly on a nondimensional acceleration param- 
eter and was practically independent of the instantaneous Mach 


Number. 


THE ANALOGY AND THE EFFECT OF THE 


GROUP VELOCITY 


ANALYSIS OF 


[' ONE ASSUMES THAT SURFACE WAVES are propagated 
across the free surface of a large body of shallow 
water over a horizontal plane bottom in such a manner 
that the vertical accelerations are negligible, then the 
horizontal motion is the same at all depths. In this 
case all infinitesimal disturbances move at a constant 
velocity of propagation given by +/gh, where h is the 
water depth. 

This motion is directly analogous to a two-dimen- 
sional compressible gas flow where V gh corresponds 
directly to the speed of sound, a, and the infinitesimal 
change in /: is directly proportional to the isentropic 
variation of either density (p), or temperature (7), in 
a perfect gas having a specific heat ratio of y = 2 (see 


reference 1) so that 


hp * 5) y—- 1 - (6) 7) 
eo -. s he | (1) 


V ygRT = V yp/p 


II 


sh=a= 


5 


% 


for y=2 only. 

This analogy is valid for either stationary or non- 
stationary isentropic flow. However, the analogy is an 
approximation that is valid only so long as the vertical 
acceleration is negligible and Eqs. (1) are true for all 
surface waves of any wave length. 


Presented at the Annual Summer Meeting, I.A.S., Los Angeles, 
June 27-28, 1951 


* Associate Professor. 


The first difficulty in the analogy between the shal 
low-water surface waves and the two-dimensional isen- 
tropic gas flow is the fact that the actual velocity of the 
surface-wave propagation is not given exactly in Eqs. 
(1), which are only approximations, but is given more 
nearly by (see reference 2) 


/gn QnT\ Oh 
+ tanh - = 


2 pi 


| 


nr?T h? 
\ eh + 2 : ( | (2) 
gpr A* 

This indicates that the vertical acceleration is not 
negligible, even when the frictional drag on the bottom is 
neglected, unless the wave length \ > h or unless h ~ 0 
when X finite. In the latter case, as  ~ 0, the second 
term containing the surface tension 7° becomes pre- 
dominant for small }. These shorter wave lengths 
correspond to the so-called capillary ripples that are 
entirely a surface phenomena and do not follow the 
analogy defined in Eq. (1). 

Fig. 1 is a plot of c vs. \ for ordinary tap water with a 
surface tension 7 = 0.005 Ibs. per ft. and various water 
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Hydraulic tank and towing carriage. 


Fic. 2. 


c>a= V gh 


independent of A, only for small water depths of the 
The 
waves corresponding to \ values smaller than this are 
entirely capillary ripples that are always unavoidably 
present in the analogy. Their effect is almost negligible 
on the formation of the hydraulic jump, analogous to 
shock waves, but they do appear in the photographs of 
the surface waves, although their relative effect may 
be practically eliminated by using a large model. 


order of '/, in. and even then only for X 2 | in. 


If the ratio of h/d is not sufficiently small, then the 
analogy fails, since there is no constant corresponding 
to the speed of sound in isentropic gas flow. Asa mat- 
ter of fact, if c varies with A, not only does the analogy 
break down but even the concept of the wave velocity 
c, Eq. (2), itself loses all physical significance. In this 
case c should actually be termed a phase velocity rather 
than a wave velocity (see reference 2), since it is neither 
the rate of propagation of the energy nor the wave 
form. Then, no given wave form can be propagated 
without a continuous change in shape, although the 
overall motion proceeds at a constant velocity called 
the group velocity. From a physical point of view, 
this so-called group velocity is always more important 
than the wave velocity, since the group velocity is di- 
rectly related to the rate of propagation of the energy 
itself. The wave velocity is only of importance when 
it coincides with the group velocity. The individual 
waves of a given wave length forming a grouped disturb- 
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ance move more rapidly than the group velocity of the 
disturbance itself, but these individual waves die out 
as they approach the front of the group, their energy 
being continually transferred to the disturbance itself 
It is only when the group velocity coincides with the 
wave velocity that a given disturbance can be prop- 
agated without a change in form, as in two-dimensional 
supersonic flow, at a constant velocity. 

It is shown in reference 2 that the group velocity jg 
given by j 


[c — X(dc/dd) | 


Therefore, the group velocity for any finite combination 
of waves with the same wave length dX is given by the 
intercept, on the c axis of Fig. 1, made by the tangent 
to the curve of c vs. \ at the given X._ It is then ob- 
vious from Fig. | that the capillary ripples, having a 
small \, travel with a group velocity greater than the 
wave velocity c. On the other hand, gravity waves, 
having a larger \, must have a group velocity less than 
c. Their group velocity is exactly one-half of the wave 
velocity in deep water (k > ) no matter how large 
the model itself may be. Only as the water becomes 
more shallow does the group velocity approach c, which 
Vv gh, 
which is independent of wave length for all \ except 


in turn approaches the constant value, a = 


those corresponding to the extremely small capillary 
ripples. 

It is also seen from Fig. 1 that decreasing the water 
depth below 0.2 in. will make matters worse, at least 
for the standard value of surface tension for pure water. 
Likewise, a 50 per cent decrease in the surface tension, 
as shown by the dotted line in Fig. 1, adversely affects 
the analogy for h = 0.2 in. 

The second important restriction on the analogy 
arises in the investigation of shock waves. The local 
water depth directly corresponds to the local value of 
the enthalpy or temperature in the two-dimensional 
compressible gas flow. Ina gas flow there is no change 
in the stagnation enthalpy (or stagnation temperature 
across an ordinary shock wave, since the flow is adiabatic 
(although nonisentropic). However, for the hydraulic 
jump or surge formed in water, there exists an actual 
loss in the stagnation water depth due to the formation 
of eddies that eventually increase the water temper- 
ature. This is an additional loss of available energy 
or increase of entropy which has no counterpart in the 
analogous ordinary shock wave in gas flow. Conse- 
quently, the shock-wave analogy is limited to that 
range of Mach Numbers and flow deflections wherein 
this loss in stagnation depth is negligible. It is ob- 
vious that this range must include the shock waves 
wherein the change in entropy can be neglected, as de- 
fined by (1)? — 1)* > 0 for normal shock waves or by 
the flow deflection 6° — 0 for oblique shock waves (see 
reference 3). Therefore, the shock-wave analogy 1s 
valid in the limiting case either as /) — | or as the 


wedge angle 6 —~ 0. However, a comparison of the 
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computed values for the pressure changes through a 
normal shock show that the analogy is well within 1 per 

cent if My < 1.5. This range is, of course, satisfactory 
ior the transonic investigation of any normal or oblique 
shock. This has been verified experimentally in refer- 
ence 4, where it is shown that for the shock-wave anal- 
ogy the water depth, as would be expected, must be 
of the order of '/, in. 

Table 1 gives the two-dimensional oblique shock- 
wave relations for the maximum flow deflection (@nz.) 
which can be produced at a given free-stream Mach 
Number .V/ in a perfect gas for y = 2, while Table 2 
gives the same relations for y = 1.4 corresponding to 
air (Table 2 is taken directly from reference 3, and Table 
|iscomputed in the same manner). 

It is immediately noted that the increase in y pro- 
duces a large increase in the strength of the shock wave, 
especially at the higher Mach Numbers. It is impor- 
tant to note that the values given in Table 1 do not cor- 
respond to the hydraulic relations for the hydraulic 
jump or surge actually occurring in the water channel. 
The hydraulic-jump relations approach the perfect gas 
relation for y = 2 (Table 1) as My— 1 and give nearly 
the same values for My < 1.5. However, for Mo > 1.5, 
some of the hydraulic-jump relations depart consider- 
ably from the perfect gas relations for y = 2, and the 
analogy therefore cannot be used for any finite model be- 
yond this range, since all of the isentropic flow regions 
are still governed by Eq. (1) with y = 2. This limita- 
tion, of course, does not restrict investigation in the 
transonic range wherein M < 4/2. 

The oblique shock-wave analogy may, of course, be 
carried to higher values of My by limiting the values of 
§so that the entropy change, which varies as 6*, is neg- 
ligible. In all cases, however, the analogy is valid only 
for y = 2, since all the flow expansion and isentropic 
compression regions are governed by Eq. (2) and not 
by the hydraulic-jump relations. The hydraulic- 
jump relations may be used alone, ignoring Eq. (2) 
and y = 2, only for the trivial case of an infinite wedge 
in an infinite stream. 

Although no real gas has y = 2, still the analogy is 
useful as an experimental verification of any theoretical 
equation developed from a study of two-dimensional 
gas flow. Attempts have been made to use other than 

rectangular water channels in order to obtain different 
For example, an equilateral triangular 
However, this is 


values of y. 
cross section gives a y value of 1.5. 
based on the a priori assumption that the hydrostatic 
pressure relation may be used or that the vertical ac- 
celeration is zero. Fig. 1 shows that this is true for 
ordinary model sizes only for water depths of the order 
of 1/, in., even in an unrestricted infinite body of 
water of uniform depth; consequently, the analogy must 
fail for the deeper water depth required in such a chan- 
nel. 

This restriction to small water depths also makes the 
use of a flowing water channel difficult, since the flow 
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TABLE 1 
Detachment Angle for Perfect Gas Relations (y = 2) 
Omaz p p H On 
( Deg.) Mo M DP p Ho ( Deg.) 
0 l ] l l l 90 
2! Bry 0.957 1.34 1.156 1.000 42.4 
5) 1.30 0.9438 1.61 1.27 0.987 68.9 
10 1.56 0.925 2.33 1.50 0.970 64.8 
15 1.90 0.932 3.44 1.76 0.900 62.3 
20 2.42 0.941 5.70 2.08 0.755 61.0 
25 3.51 0.965 12.4 2.48 0.521 60.2 
30 x l x 3 0 60 
TABLE 2 
Detachment Angle for Perfect Gas Relations (y = 1.40) 
Rene pP p H Ou 
(Deg.) My M Po po Hy ( Deg.) 
0 ] l I l l 90 
2'/s 1.14 0.961 1.24 1.17 0.999 74.1 
5 1.24 0.944 1.43 1.29 0.995 70.8 
10 1.42 0.925 1.84 1.54 0.972 67.5 
15 1.61 0.919 2.36 1.81 0.939 65.7 
20 1.84 0.920 3.07 2.14 0.870 64.9 
25 2.13 0.928 4.13 2.55 0.758 64.6 
30 2.52 0.941 5.90 3.06 0.591 64.8 
x 6 0 67.8 


over the fixed bottom cannot be sufficiently uniform 
when the optimum depth is only '/, in. Therefore, 
the following investigation was made in a rectangular 
channel 4 ft. wide and 20 ft. long, where the model was 
towed by means of a carriage. The section of the tank 
bottom where the pictures were taken was made of 
A description of the apparatus used, as shown 
The photo- 


glass. 
in Fig. 2, is given in references 5 and 6. 
graphs were taken by having the light source below the 
glass bottom of the tank and the camera above the 
water surface. The camera was focused on a ground- 
glass screen supported above the model and | in. above 
the water surface. 

Fig. 3 was obtained by using a high-intensity carbon 
arc-light source, while Figs. 4, 5, and 6 were obtained 
by using a Kodatron speed lamp. In Fig. 3 the grid was 
below the glass bottom, while in the later pictures the 
grid was marked on the ground-glass screen, which was 
above the water. 

The present tests were taken for steady transonic flow 
near M, = 1 and for accelerated flow, since these are 
the most difficult to test experimentally in ordinary 
wind tunnels but present no difficulty with a towing car- 
riage in a still water channel. 

It is again stressed that these cannot be considered 
as model tests but only as an experimental check of any 
theoretical compressible gas-flow equation wherein 
the numerical value of y is taken as 2. There is no 
necessity of attempting to correlate any parameter 
other than the Mach Number for steady flow, since 
there is none in this analogy. For example, correlating 
the Reynolds Number would have no physical signifi- 
cance because the boundary conditions are entirely 
different. The only logical basis of comparison is the 
one available from neglecting all viscosity effects. 
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DISCUSSION OF THE EXPERIMENTAL ResvuLts 


Fig. 3 shows the supersonic flow pattern about a 
*/, in. thick wedge having a nose angle of 20° (9 ~ 
10°) in a still-water depth of 0.33 in. The equivalent 
to the free-stream Mach Number (1/9) is defined by 
the carriage speed (I) and the still-water depth (h) as 


My = V/V gh 


The theoretical minimum value of JJ, necessary for an 
attached shock wave to exist at the nose of a 20° wedge 
is 1.56 for a compressible gas having y = 2 (Table 2), 
while the theory for a hydraulic jump or surge requires 
My = 1.51. : 

Fig. 3a shows the equivalent of a detached shock 
wave in steady flow for 1J) = 1.29, while Fig. 3b shows 
the shock wave already at the nose but still curved, 
proving that it has just become attached (see reference 
3) at My = 1.48. 
agrees more closely with the hydraulic-jump theory 


As would be expected, this result 


than it does with the compressible gas theory (y = 2), 
However, this difference does not exist in the transonic 
range defined by My < 7/2. For example, a 10° 
wedge has the same minimum JJ, of 1.30 for an at- 
tached shock to exist either in a compressible gas 
(y = 2) orasa hydraulic jump. The experimental re- 
sults for a wedge having 6 = 5° closely checked this 
value (see references 6 or 7). 

Fig. 3c presents a completely supersonic attached 
shock wave at My = 1.52 and points out a further diffi- 
culty in the analogy—namely, the relatively large 
thickness of the hydraulic jump or surge as compared 
to that of the ordinary shock wave in air. This thick- 
ness effect invalidates the analogy unless the shock- 
front radius of curvature is always large. Fortunately, 
this is generally the case in transonic flow, especially 
with a detached shock wave. 

Fig. 3d illustrates the effect of acceleration on the 
shape of the detached shock wave, which now con- 
tinually The 
celeration parameter is defined by the acceleration (A) 
multiplied by the nose length (C) and divided by the 


varies in shape. nondimensional ac- 


speed of sound (a) squared, where 


(AC/a*) = M)?(AC/V?) 


The larger this parameter, the more the distortion of the 
flow pattern (see reference 7). 
Fig. 4 demonstrates the effect of water depth upon 
The pic- 
y-in. thick 


the wave formation and the group velocity. 
1.2 with a 
wedge having a nose angle of 30° (6 = 15 
that for water depths greater than the optimum of 
'/, in., the main surge front representing the detached 


tures were all taken at JJ) = 
It is seen 


shock wave is pushed further from the nose and loses 
the original smooth curvature. Fig. 4b indicates the 
formation of the secondary wave fronts that alter the 
shock 


These secondary wave fronts correspond to shock waves 


shape and location of the detached analog. 
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at an equivalent free-stream Mach Number 


formed 
it is actually greater than that given by My = V 


The expression for J) only applies to the wave 
gths that have a horizontal tangent (dc/dX = Q) in 
Fig. |, so that the wave velocity is equal to the more 
For the other wave lengths 


the 


V gh. 
Jen 


important group velocity. 
the wave velocity itself has no direct effect, and, as 
previously discussed, the group velocity must be used. 
Therefore, the equivalent free-stream Mach Number 
for any group of these wave lengths would be given by 


M, = V/[e — X(dc/dy)] 


which is greater than ./, for the shorter wave-length 
gravity waves and less than Mp for the capillary ripples. 
As previously pointed out the group velocity for any 
given Wave length is obtained by the intercept on the 
vertical c axis of the tangent to the curve in Fig. 1. 
Consequently, different wave lengths, all having dif- 
ferent wave velocities, would still have the identical 
group velocity if they had a common tangent. Then, 
if a whole series of wave lengths had approximately the 
same group velocity, any small disturbance composed 
of them would move at this constant group velocity, 
which would be equivalent to the speed of sound, al- 
though the shape of the disturbance would vary with 
time unless the wave velocities were also equal, as pre- 
viously explained. A strong wave front or surge, cor- 
responding to a steady-state shock wave, could be 
formed only if a considerable portion of the curve of ¢ 
vs. \ (Fig. 1) formed a nearly straight line. This is so 
because any finite body essentially produces an in- 
finite number of infinitesimal disturbances, all traveling 
with their particular group velocity. These infini- 
tesimal disturbances can produce a finite surge moving 
at the constant speed V > +/ gh and corresponding to a 
steady-state shock wave only if a finite range of wave 
lengths have approximately the same group velocity. 
If the water depth is greater than '/, in., it is therefore 
possible to have several common group velocities, cor- 
responding to several different speeds of sound. If the 
capillary ripples are neglected, then the longer gravity 
waves would have the highest equivalent speed of 
sound given by +/gh, while the shorter gravity wave 
lengths would have an equivalent speed of sound con- 
siderably less than the minimum wave velocity. Pre- 
vious analyses of the analogy (e.g., reference 8) have 
been in error because they did not realize the impor- 
tance of the group velocity and because they did not 
stress the fact that a surge is a finite disturbance that 
travels at a constant speed much greater than that of 
an infinitesimal disturbance. 

As the depth is increased above '/» in., the main surge 
becomes more distorted until, finally, because of the 
large variations in the group velocities, no steady-state 
surge can be maintained. This unsteady state was 
evident when Figs. 4c and 4d were taken. 

Fig. 5 depicts the effect of model size and surface 


tension. The same geometric shape is used at the 
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FIG. Sb: h=0.388 





























FIG. Sd: h* 0.255, T reduced 50% 


FIG. 5c: h* 0.570 


@r16*, M,*1.20, t#3in. 


same Wp as in Fig. 4, the wedge now being 3 in. thick. 
It is seen that the larger model size has produced a rela- 
tive decrease in the appearance of the capillary ripples, 
since they correspond to extremely small values of X. 
Similarly, the water depth must now be slightly greater 
before the secondary wave fronts become as strong. 

Fig. 5d shows the effect of decreasing the surface ten- 
sion from 0.005 to 0.0025 Ibs. per ft. by adding a wetting 
agent (Oronite). Strong secondary wave fronts are 
again formed, even though the water depth remained 
at its previous optimum value. This is explained by 
referring to the dotted line in Fig. | which corresponds 
to the 50 per cent decrease in surface tension for a water 
depth of 0.2 in. The dotted’line shows that all \ from 
approximately '/; to 1 in. and from 1 to 2 in. have almost 
a set of common group velocities and could produce a 
fairly strong secondary surge. 

Fig. 6 shows exactly the same effect of water depth, 
model size, and surface tension on a circular cylinder. 
The conclusion is that the optimum water depth for 
pure water is about '/, in. and that the model size 
should be as large as possible in order to have the larger 
wave lengths predominate and thereby decrease the 
relative effect of the shortest gravity waves and the 
However, the model size is limited 


capillary ripples. 
For example, a 7!/>-in. 


by the width of the tank. 
diameter cylinder was found to ‘“‘choke”’ this 4-ft. wide 
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FIG. 6b: h= 0.500, D#1.25 
































FIG. Gc: h=0.255, D*2.5 FIG. 6d: h*0.255, D*2.5, 


T reduced 50% 


CIRCULAR CYLINDERS AT Mo*!.20 _ 
tank for all Mo < 1.5 (see reference 5). This choking 
phenomena was easily observed as the cylinder started 
moving. First, the detached shock wave became 
strong enough to form a Mach reflection at the walls. 
Then the normal shock portion of the Mach reflection 
increased in length until it completely spanned the 
channel moving upstream, ahead of the cylinders. 


COMPARISON OF DETACHED SHOCK WAVE LOCATION 
WITH THEORY 


Because the detached shock wave is in the transonic 
region and especially because it does not have a sharp 
radius of curvature, it should provide the most reliable 
data to test the analogy to a perfect gas with y = 2. 
Figs. 7 and 8 present the experimental measurements of 
the detachment distance along the line of symmetry 
for the wedges and the cylinder at various values of 
My. It was found in every case that increasing the 
still-water depth moved the detached shock consider- 
ably further out. 
noticeable effect of model size on the detached shock 


In the transonic range there was no 


location, at least for water depths near '/, in. How 
ever, the 50 per cent decrease in surface tension seemed 
to move the detached shock slightly nearer to the circu 
lar cylinder (see reference 5). 

Unfortunately, there is no satisfactory theory yet 
available for predicting the detached shock location. 
Two approximate theories, given in references 9 and 10, 
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are included for comparison. The simplest approxima. 
tion is the so-called geometric method of Moecke|s 


which is based on the assumption that 


y,/L = tan Omar. (9 


” 


where Ona, is the maximum stream-tube deflection that 
can occur at a given My (Table | or 2). y, is the vertj. 
cal distance from the sonic point on the body surface 
to the axis of symmetry of the body. The normal 
shock or nose of the detached shock wave is then at 
the distance L from the location of y,. For a wedge the 
sonic line appears at the shoulder, so y, is one-half the 
thickness (see Fig. 7). For an arbitrary body the 
sonic line can be assumed to occur at the point where 
the surface slope is equal to On7,; consequently, for a 
circular cylinder the distance 1; from the centerline 
to the detached shock is given by 


(D/2)/1, = sm 0... (4 


The predicted distances given by Eq. (4) are pre- 
sented in Fig. 8, using the values of 092, for a com- 
pressible gas with y = 2 (Table 1) and those corte. 
sponding to the hydraulic-jump theory. As would be 
expected, the two curves approach one another for 
My < 1.5 but diverge considerably at the higher Mach 
Numbers. As previously mentioned, the values of 
Omar. for a hydraulic jump exactly approach those of 
Table | for My < 1.5 but are more closely approximated 
by those of Table 2 for slightly higher Mo. 

As pointed out by Moeckel,’ Eq. (3) is more nearly 
valid for blunt’ bodies and for transonic Mach Numbers 
not too near unity. This seems to be the case in Figs. 


7 and 8. However, the theory of reference 9 predicts 
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Fic. 8. Effect of water depth and model size on the shock detach- 
ment for circular cylinders. 


no effect of the nose shape forward of the sonic line, or 
wedge shoulder, and Fig. 7 proves there is a considerable 
effect for slender wedges near a Mach Number of unity. 
The theory presented by the author in reference 10 is 
more applicable in this case, since it is an approximate 
first-order solution that is valid for extremely slender 
wedges near a Mach Number of unity. This theory (see 
reference 7 or 10) yields the expression 
Vs 


ee ees (M ad (5) 
‘a = = meq ie o 
7 a alll eS 0 


for the wedge shown in Fig. 7. Then, for the analogy 


with a compressible gas having y = 2, 


: = (1 + ) (6) 
y, 6 exp {5.92 [(My — 1)7/0]} — | 


As would be expected, Fig. 7 shows that Eq. (6) 1s in 
better agreement with the experimental data nearer a 
Mach Number of unity, especially for the slenderest 
wedge. 


THE EFrect OF CONSTANT ACCELERATION ON THE 
DETACHED SHOCK WAVE 


The main effect of acceleration, as indicated in Fig. 
3d, was to produce a nonsteady detached shock-wave 
shape that seemed to vary almost periodically from a 
rounded nose to a flattened, normal shock nose. In 
order to study this phenomena, a series of photographs 


was taken of the strong detached shock wave existing 
in front of a rectangular block 2 in. wide (Fig. 9). The 
model was started from rest at constant acceleration. 
‘Lhe pictures were obtained by mounting a motion-pic- 
ture camera on the towing carriage itself. 

Fig. 9 shows a sequence of photographs at intervals 
of approximately 0.078 sec. The complete series of six 
pictures was approximately the cycle representing the 
periodical change in shape of the detached shock. As 
the Mach Number increased, at constant acceleration, 
periodic secondary wave fronts formed between the 
body and the detached shock. These secondary surges 
gradually overtook the detached shock, flattening its 
nose and keeping it at almost a constant distance ahead 
of the body, independent of the increasing Mo. As 
seen in Figs. 7 and 8, the detached two-dimensional 
shock wave, in steady flow at a constant velocity, is 
nearer the nose at the higher Mach Numbers. On 
the other hand, it is known that for the one-dimensional 
flow ahead of a moving piston the shock waves are fur- 
ther ahead of the piston at the higher speeds and still 
further with a finite constant acceleration. In the 
present tests it was found that the detached shock 
wave ahead of a two-dimensional rectangular block, 
starting from rest at a constant acceleration, approxi- 
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mately followed the steady-state detachment distance 
for small values (less than 0.01) of the nondimensional 
acceleration parameter Ay,/a*. However, as the mag- 
nitude of this parameter was increased, the detachment 
distance became more nearly a constant, practically 
independent of the instantaneous Mach Number. The 
tendency of the detached shock wave to move toward 
the body as the Mach Number increased seemed to be 
counteracted by the energy supplied to it by the periodic 
formation of the secondary surges formed downstream 
and continually overtaking and joining the main de- 
tached shock front. 

Further studies are now in progress to determine the 
reason for the apparently periodic fluctuation and to ex- 
tend the validity of the analogy to other nonsteady flow 
problems. This seems to be the most useful service of 
the analogy, if valid, since it is a field that is difficult to 
investigate in ordinary two-dimensional gas flow. 
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Analysis of Nonuniform Columns and Beams 
by a Simple D.C. Network Analyzer 


GEORGE W. SWENSON, JR.* 


Washington Unwersity 


SUMMARY 


The analysis of beams, columns, and beam-columns whose 
flexural rigidities and/or lateral loadings are arbitrary functions 
of position can be solved with the aid of a simple electrical 
network analyzer. In the case of the column, the lowest critical 
load can be determined with accuracy. The 
electrical network is of ladder configuration, with displacement 
represented by node-pair voltages and moments by currents. 


considerable 


The column problem requires the use of a number of ‘‘negative 


resistors,’ which are simulated by subnetworks of positive 


resistors and adjustable voltage sources. 
Problems are discussed which illustrate the application of the 


analyzer. In the case of a column with a single abrupt transi- 
tion in rigidity, the value of the critical load is determined to 


within 0.5 per cent of the analytically calculated value. 


INTRODUCTION 


— DETERMINATION OF THE DEFLECTION of a non- 
uniform beam or beam-column or of the critical 
load of a nonuniform column can be simplified and 
expedited by the use of an electrical network analyzer. 
The general method of attack is to replace the differen- 
tial equation of bending by a set of approximately 
equivalent finite-difference equations, a set of linear 
algebraic equations that must then be solved simul- 
taneously to yield the deflections at a number of 
equally spaced stations along the member. For ac- 
ceptable accuracy, one must write a relatively large 
number of equations; therefore, the computational 
labor involved in their solution is great enough to 
warrant the use of some mechanical aid. The similarity 
in form of the difference equations to the current equa- 
tions for a lumped-parameter electrical network suggests 
that such a network be used as an analog to the problem 
at hand and that the solution be obtained by measure- 
ment of the electrical quantities. When such a possi- 
bility is further investigated, it is found that both 
positive and negative network parameters are required. 
This suggests, in turn, the use of alternating currents, 
with inductive reactances as the positive parameters 
and capacitive reactances as the negative parameters 
or vice versa. However, for reasons of accuracy and 
economy, it seems desirable to use an all-resistance, 
direct-current network, if possible, and it has proved 
feasible to synthesize the required ‘“‘negative resistors”’ 
from positive resistors and voltage or current sources. 
Received August 10, 1951 
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The result is a network analyzer capable of solving 
the general slender beam-column problem rapidly and 
accurately by an iterative process of adjustment of the 
network. The saving in labor conventional 
numerical techniques is considerable, and the device 


over 
is useful either for synthesis or analysis 


NOMENCLATURE 


The following symbols are used: 

E = Young’s modulus, in general, a function of x 

h = length of finite-difference interval 

I = second moment of area, in general, a function of x 

tn = forcing current inserted at mth node 

l = length of beam or column 

M = bending moment, in general, a function of x 

P = axial load (force) on a column 

Por. = critical load of a column 

R- = base resistance of network analyzer 

Rox resistance connected between mth node and ground 

V, = node-to-ground voltage at mth node 

x = space variable in direction along member 

y = lateral deflection of member 

Bn = Ph?/(EI), — 2, coefficient of y, in the nth difference 
equation 

on = (M/EI),, normalized bending moment in the nth 


difference equation 


THE DIFFERENCE EQUATIONS AND THEIR ELECTRICAL 
ANALOG 


The differential equation of bending of a slender 


beam-column is 


(d?y/dx?) + (P/E]Dy = M/EI (1) 


where the coefficient of y and the term on the right- 
hand side may vary with x. Referring to Fig. 1, the 
approximation to the second dexivative at point 0, 
taken from the calculus of finite differences, ' is 


(d?y/dx")o ~ (Ya + VY» — 2¥0)/h? (2 
Thus, at point 0, Eq. (1) is approximated by 
Va t+ M t+ volhPP/(ED)o — 2) = W(M/ED)) (8 


A similar equation can be written at each of the equally 
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The finite-difference lattice 
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the rheostat R, is adjusted manually until the voltmeter is large 
4 V reads zero. Under this condition, a current V,/R. it is p! 
V; ba flows out of the circuit, info the node of the analyzer Whe 
so that the circuit simulates the desired negative colum! 
R resistor. negati 
R In the present analyzer, the resistors R,, are arranged found 
22 ' : 
to be plugged into jacks on the analyzer panel. |p. js to 
se expensive volume control potentiometers are used grounc 
a their values being adjusted on a Wheatstone bridge that p 
Fic. 2. The d.c. network analyzer. during the setting-up process. Though 15 or » the ad 
junction points would probably be sufficient for mog a curt 
spaced stations on the member. As / becomes vanish- problems, the present analyzer has provisions for soly- repres 
ingly small, these equations become exact. After a ‘8 100 simultaneous opuncrens, having been con. =. 
few trial problems have been solved, experience structed originally for application to the two-dimen. = ’ 
indicates the number of intervals necessary to assure sional wave equation.” Because adjustment of one the ar 
the desired accuracy. As an example of the complete ° the resistors R, affects the balance of each of the o . 
setup of a problem, the following matrix equation other negative resistance circuits, an iterative adjust. rent 1 
represents a member divided into four intervals, with ™€"t process must be applied until all the voltages | -— . 
both ends (stations | and 5) assumed to be fixed: are zero. The process converges for all axial loads Pers 3 
less than, and also for a small range just above, critical the bi 
B, 1 O ye pe load. detert 
I Bs | ys | = | 3 (4) The results described below could be duplic ited on produ 
O 1 Pst LM pi a nine-equation analyzer whose fixed resistors are of trary 
The electrical analog of this system is most easily | per cent tolerance. For measuring the voltages, an the vi 
set up on the basis of Kirchhoff’s current law, which electronic d.c. voltmeter is needed, capable of indicating 
states that the sum of the currents entering a junction O° millivolt. A value of R of about 100 ohms is 
point, or node, must be zero. Fig. 2 illustrates the ComyEMneNs. Because all the per ene: volt 
network with the desired properties. The current-law 8° supplies have one common terminal, they _ be A] 
equation that applies at node 0 is derived from the same source, which should preferably term 
be a battery of storage cells for constancy of voltage. 
¥o—~ Ve 4 Vy — Vo Vo _ in (5) About 30 volts is adequate. It is estimated that such 
R R Rw an analyzer could be built at a cost for materials of was i 
or less than $100, exclusive of battery, voltmeter, and and t 
V, + V, — Vo(R/Ro + 2) = oR Wheatstone bridge. 
Thus, V’,, = Yq, provided that DETERMINATION OF CRITICAL LOAD This 
Ron = —R(ET),/(h°P) (6) The difference equations of a slender column with -_ 
and no lateral loads are those of Eq. (4) with the ¢’s equal a 
1, = h?(M/ED),/R (7) to zero. Such a system has a nontrivial solution only eg 
“fp . e ° — . ) 
An exception to Eq. (6) occurs at the station adjacent if its determinant vanishes, and this condition might ee 
to a station whose displacement is assumed to be zero. wes oe oe ites eteeennae tages tle Nera how 
In this case, the network at the “adjacent” node has numerical labor is great when the number of stations | 
the configuration of node 2 in Fig. 2, corresponding to 
the first row of the matrix in Eq. (4), and Tt 
Rx» = R/{1 — W®P/(ED)»2] (8) scrib 
THE NETWORK ANALYZER | 
In Eq. (8), Ry» is positive for values of h which give y 
reasonable accuracy. R,,, however, representing the “a 
node-to-ground resistance at a point in the interior | _ 
of the network, is negative for any value of . Such = 
a negative resistor can Le simulated by the circuit of thay 
Fig. 3, in which R, is a resistance whose value is equal | the 
to the magnitude of the desired negative resistance, 
R,», and in which the resistors R’ have arbitrary but This 





equal values. With the circuit in place in the analyzer, Fic. 3. The negative resistor circuit 
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ANALYSIS BY A D.C. 


enough to yield acceptable accuracy, so that 


is large 
it is preferable to get the results from the analyzer. 


When the analyzer is set up to represent such a 
column, the value of P being arbitrary, and when the 
negative resistances are systematically adjusted, it is 

The remedy 


found that a trivial solution is obtained. 
is to apply a battery between some one node and 
ground to force the network to assume a solution at 
that point conforming to the applied voltage. When 
the adjustment process is completed, it is found that 
a current exists in the battery circuit. This current 
represents the lateral load at that point required to 
cause a deflection corresponding to the battery voltage. 
This is true for P either less than or greater than P,,.; 
the analyzer does not concern itself with the fact that 
the deflection is unstable in the latter case. The cur- 
rent is positive (toward the node) for P less than P.,. 
and negative for P greater than P.,, For P equal to 
P.,., a finite deflection requires no lateral force, so that 
the battery current is zero. Thus, the critical load is 
determined by graphing the required to 
produce a certain node-to-ground voltage at an arbi- 
trary station in the network and by interpolating to 
the value of P which corresponds to zero current. 


current 


EXAMPLE I. CRITICAL LOAD OF A 
TAPERED COLUMN 


A pin-ended column whose flexural rigidity has the 
form 
(9) 


was investigated by the above method, using h = //10, 
and the critical load was found to be 


Po, = 14.3(EI)o/I? (10) 


This is 3.4 per cent lower than that of a uniform column 
with the same average flexural rigidity. The time 
required for this solution, including calculation of 
coefficients, setting up the network, and measurement 
of deflections for several different lateral loads, was 


approximately two hours. 


EXAMPLE II. CRITICAL LOAD OF A 
STEPPED COLUMN 
This pin-ended column has a flexural rigidity de- 
scribed by the equations 
EI = (EI), 0<¢x<l1/2)\ 
EI = 1/8(EI)p, W2<xK<lf 


(11) 


Because one station was located at the point of transi- 
tion, the rigidity at this point was arbitrarily made 
equal to the average of the other two values. Again, 
using ten finite-difference intervals (nine equations), 
the critical load was determined to be 

Per. = 0.196(ED)o/?? (12) 


This value checks within 0.5 per cent of the previously 


NETWORK ANALYZER 
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Deflection of the stepped column. 


published analytical solution,* obtained with the aid of 
the Laplace transformation, and within 1.7 per cent of 
an experimental value* for a Douglass fir column. 


HIGHER MODES OF BUCKLING 


Should it be desired to know the critical loads for 
higher modes of buckling, they can be determined by 
a trial-and-error process. The column is divided into 
sections corresponding to an initial estimate of the 
locations of the points of zero displacement. Then 
the critical loads of the individual sections are deter- 
mined by means of the analyzer, and the locations 
of the zero points adjusted to make the respective 
critical loads more nearly equal. This process is 
continued until the desired accuracy is attained. The 
critical load for the mode in question is then equal to 
the critical load of any one of the sections. 

This scheme, and indeed the entire technique of the 
analyzer herein described, can be applied equally well 
to problems in the oscillations of nonuniform strings, 
resonators, coaxial electro- 


plane-wave acoustical 


magnetic resonators, and so on. 


BEAM-COLUMNS AND BEAMS 


It is necessary to apply a simulated lateral load in 
determining the critical load of a column; thus the 
solution of such a problem also includes the solutions 
of one or more beam-column problems. For example, 
Fig. 4 illustrates the deflection of the stepped column 
under the influence of a lateral load applied at the 
center and an axial load of slightly less than the critical 
value. Deflections produced by arbitrary distributions 
of lateral loading are easily determined by inserting 
currents at the several node points of the network, 
according to Eq. (7). 

The beam problem in the absence of axial load is 
particularly simple. In this case the term in y of 
Eq. (1) is absent, so that no node-to-ground resistors, 
either positive or negative, are required. The solution 
is available immediately when the appropriate currents, 
as determined from Eq. (7), are inserted into the 
appropriate nodes. These currents are most con- 
veniently obtained from a common high-voltage power 
supply through dropping resistors whose values are 
high compared with the base resistance of the net- 
work. 

(Continued on page 278) 





Reply to Discussions of the Wright Brothers 
Lecture’ 


P. B. WALKERT 
Royal Aircraft Establishment 


HE DISCUSSIONS FOLLOWING the Wright Brothers 

Lecture have added much valuable material to 
the original paper. Though there are minor differences 
of opinion on some matters, the discussions as a whole 
are remarkable for the general measure of agreement on 
fundamentals. The comments which follow are mainly 
directed towards answering specific questions or in 
response to requests for further information. 

Mr. A.B. Callender asks how the fatigue test criterion 

an alternating load of 7'/2 per cent of the design 
ultimate, applied two million times—was obtained. 
It is, of course, difficult to trace all the tentative proc- 
esses leading to an empirical criterion, but there are 
a few points which may be of interest. 

The single fatigue test has to represent a wide range 
of loading conditions, corresponding to an extensive 
working portion of the endurance (SN) curve. From 
this standpoint, the 7'/. per cent alternating load is 
chosen as being well placed in a sensitive region of the 
typical curve and as not calling for undue length of 
testing time. In addition, the criterion has to cover 
in a somewhat indefinable way the cumulative effect 
of mixed cycles. In this connection it is of interest 
to record that recent British calculations indicate that 
7'/2 per cent for the alternating load is a little on the 
high side theoretically, but not more so in practice 
than is required to give a reasonable working margin 
of fatigue strength. 
Mr. Callender’s figures, though it is endorsed by a 
later speaker (Mr. Hitchcock). 

It is intended that the complete criterion be modified 


This conclusion is at variance with 


as necessary with increasing experience or to meet 
special conditions. In general, it is simpler to alter the 
required number of cycles than the alternating load. 
The two million’ cycles stipulated is regarded as cor- 
responding roughly to a life of 30,000 flying hours or 
about 10 years. Within reasonable limits, however, 
proportionate reduction is permissible when a shorter 
life is expected for reasons not concerned with fatigue. 
Adjustments are also to be expected when fatigue load- 
ing data are available for particular aircraft operating 
over particular routes. The original criterion was in- 
tended to set a conservative standard at which designers 
were advised to aim. Later experience, obtained since 


* Walker, P. B., The Experimental Approach to Aircraft Struc- 
tural Research, Journal of the Aeronautical Sciences, Vol. 19, No. 
3, p. 145, March, 1952. 

t Head, Structures Department. 


the lecture, indicates that the criterion may not always 
err on the side of safety. 
Professor N. J. Hoff’s ccmments are most encouraging, 
In a deliberate presentation of a case for greater recog- 
nition of experimental work in the structural field, 
there was a danger of appearing to minimise the value of 
theoretical work. 
course, no conflict between the two branches of stryc- 


mathematical and There is, of 
tural research, but it is good to have this view endorsed 
by an outstanding leader in theory of structures. 

Mr. Paul Kuhn, whilst recognising the value of small- 
scale work, expresses doubts as to its accuracy. I do 
not share these doubts. As stated in the original 
paper, the model work is not regarded as a method for 
absolute prediction of stresses in real structures, but 
British experience leads to the conclusion that it is 
reliable for comparative work such as the examples 
given in the paper. 

On the question of comparison between “‘aerody- 
namic’ and “‘structural”’ ribs, Mr. Kuhn gives N.A.C.A. 
results conflicting with those of R.A.E., and suggests 
that this discrepancy is an indication of the margin of 
I would suggest some other expla- 
nation, though I cannot speak with certainty in the ab- 


experimental error. 


sence of detailed knowledge of the N.A.C.A. tests. 
The R.A.E. models had ribs designed to represent 
typical ribs as used on actual aircraft. Heavier ribs, 
especially ribs with heavy flanges, might behave differ- 
ently; but in this event a direct comparison becomes 
difficult and results are open to different interpreta- 
tions. According to R.A.E. calculations, moreover, 
the participation stresses mentioned by Mr. Kuhn are 
not significant for ordinary light ribs, and for heavy 
ribs become significant only at angles of sweepback 
exceeding 60°. The R.A.E. tests, incidentally, did not 
go beyond 45°. 

Mr. Kuhn's suggestion that the engine loading effect 
could have been predicted from theoretical work al- 
ready published in England is not strictly correct. 
The relevant portion of this work dealt only with a 
The effect of the offset engine load 
is not so simple, and it is much greater than could have 


change of torque. 


been accounted for by the then existing theory. 

In reply to Mr. Kuhn's query about the future use 
of the counting accelerometer, it can be stated that this 
will be installed on as many aircraft as operators will 
tolerate. With increasing realisation of the fatigue 
danger, it may well be on the way to becoming a piece 
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REPLY TO DISCUSSIONS 


wndard equipment for fitting to all civil and trans- 


of ste 
port aircraft. The point about speed being signifi- 
as well as acceleration is fully appreciated, and is 


cant 
covered in my reply to Mr. Hitchcock. 

Mr. Walter Ramberg raises some vital points about 
design and test interpretation which unfortunately 
cannot be treated adequately in brief comments. As 
I understand the term, ‘general instability’ did not 
occur in any of the tests analysed. Nevertheless, it 
does not follow that the only alternative was local 
weakness which “‘could and should have been caught 
before the test of the complete structure.’ Stress 
diffusion, as mentioned in the paper, and other allied 
effects often lead to premature failure of a structure 
whose components would have been regarded as satis- 
factory when treated individually. 

As regards Mr. Ramberg’s next point, it is not gener- 
ally possible to make a clear-cut distinction between 
failures which could have been prevented by careful 
design and construction and those which could not 
In the first place, the economics 
The prospects of 


have been predicted. 
of design comes into the picture. 
predicting any form of failure depend upon the time, 
money, and effort that can be expended in design re- 
search and theoretical study. Without setting precise 
limits to these, it is impossible to divide premature 
failures generally into predictable and unpredictable 
classes. In the second place, it has to be borne in 
mind that the designer normally makes every endeavour 
to produce a satisfactory structure, and does not hesi- 
tate to make laborious calculations when these are 
The result is that—apart 
premature failure, 


thought to be necessary. 
from the few unqualified mistakes 
when it occurs, frequently admits of no simple explana- 
tion and is attributable to a complexity of causes, 
Faced with such a situation, it becomes difficult for 
the structural specialist to say exactly what it is that 
the designer is supposed to have predicted. 

Professor E. EF. Sechler does a great service in giving 
a warning against unrestrained enthusiasm for analog 
and digital computing machines. Development in this 
field is extremely important, but must not be allowed 
to reduce interest in the direct experimental approach. 
The two approaches are not alternatives, and both 
increase in value when they are followed in parallel. 
Incidentally, in applying the Xylonite-model tech- 
nique there are times when a model can almost be re- 
garded as an analog machine rather than as a physi- 
cal representation of a structure. 

Mr. J. F. McBrearty appears to be seriously con- 
cerned at what he regards as the over-ready acceptance 
of the design philosophy by which structures are tested, 
broken, modified and re-tested as a matter of deliberate 
policy. I entirely agree with him in pointing out the 
dangers of unrestrained use of the principle which under- 
lies this philosophy. As stated in my paper, it is 
possible to produce a structure that is altogether too 
weak, and beyond remedial action, whilst the time 
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expended in repairs and re-tests might well become 
I do not agree, however, that the principle 
In Britain, 


intolerable. 
is being misused to any appreciable extent. 
ihe evidence indicates that a good practical compro- 
mise is being achieved, and I would have thought that 
much the same applies in America. 

I entirely agree with Mr. McBrearty that testing is 
no substitute for a proper stress analysis. In my view, 
testing of a structure that has not been fully investi- 
gated theoretically is largely a waste of time, money, 
and effort. Testing is a safeguard against direct 
mistakes remaining undetected, but this is only part 
of its purpose. It should be regarded as a check on 
theory, of a kind which will benefit not ony the parti- 
cular aircraft tested but many others which follow. 

Mr. L. M. Hitchcock's remarks on ordinary static 
strength testing are of great interest, and a welcome 
additon to my own records. There is only one point 
on which comment appears appropriate, and that con- 
cerns the advisability of having a central testing unit 
available to the whole country. The situation in 
U.S.A. is different from that in England owing to the 
much greater distances involved. Even in Britain, 
however, it is considered that the central unit should 
not operate in competition with industry, and that 
its contribution must be additional and complementary 
to industry's own effort. It is considered that there 
is a purpose to be served by a central advisory body in 
the strength testing field, but that such a body cannot 
function realistically without first-hand testing expe- 
rience. Furthermore, there are long-term and uneco- 
nomic research projects which industry cannot normally 
be expected to undertake. 

On the question of testing pressure cabins in water, 
I think Mr. Hitchcock overestimates the difficulties. 
Pressure testing equipment, whether permanent or 
improvised, is much more easily produced than might 
at first sight be expected. It is, moreover, not at all 
difficult to arrange fatigue tests to reproduce the rais- 
ing and lowering of cabin pressure as occuring on air- 
craft in operational use. Successful fatigue experi- 
ments, in fact, have been made in the water tank illus- 
trated in the paper. 

The value of information additional to that given by 
the counting accelerometer is appreciated. There is 
a more elaborate version which records photographi- 
cally speed, time, and altitude, as well as acceleration 
counts. This is necessarily heavier and more expen- 
sive than the simple one, but does not need to be fitted 
so frequently. It should be understood that, as Mr. 
Paul Kuhn stated, this additional information is re- 
quired for the benefit of subsequent designs or for the 
same aircraft operating under different conditions. 
The simple version covers adequately the normal 
fatigue hazards of the particular aircraft on which it 
is fitted. 

There are a number of direct queries put by Mr. 
Hitchcock on fatigue which can be answered briefly. 
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(a) It is intended eventually to cover a range of both life, but then the life increases rapidly to something 
high and low temperatures in the testing programme of the order of 100 times. 

for the Meteor tail plane. (b) As to the choice of Dr. Karl Arnstein’s comments reached me later 
materials, I prefer not to comment on the particular than the others, and there is only one specific point 





materials mentioned, but feel justified in saying that on which I would like to comment. It is true tha 
in future more consideration will have to be given to scale-model work, just as mathematical analysis | 
the fatigue properties of materials and less to ultimate cannot when directly applied be more accurate than 


static strength—at least for civil and transport air- the simplifying assumptions. The main object jn 
craft: (c) There is a nickel-chromium steel (British using the models for comparative work as described in | 
Specification D.T.D.247) which has about the same my paper, however, was to nullify most of the errors 


coefficient of expansion as aluminium alloy. Itappears introduced by simplifications. With reasonable care 
to have been used at one time for dealing with a similar these can usually be made to affect equally the two 
problem for engines, but is probably unsuitable for models that are being compared, and so bring out 
bolt manufacture. Research to discover a suitable — thedifferences relating to the one particular phenomenon 
steel might well be justified. (d) The interference fit or structural principle that is being studied. 








found necessary to produce increased life of joints I am grateful to Dr. Arnstein for his useful summary Sec 
having steel bolts loaded in shear is approximately of the whole situation for present-day structural work. ra 
three thousandths of an inch for '/2-in. and */,-in. It is true that ‘‘many roads lead to Rome,” and, while | Deps 
diameter steel bolts. It is interesting to remark that we should all study the methods of others, success js Ur 
as the negative tolerance is progressively increased possible only with those methods in which we have om 
from zero there is at first no significant change in fatigue gained confidence through direct experience. | 
T 
flow 
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ond: 
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tion 
Analysis of Nonuniform Columns and*Beams by a Simple D.C. oat 
Network Analyzer shot 
(Continued from page 275) a 
CONCLUSIONS adaptable to any statically determinate problem in- — 
The network analyzer technique for nonuniform volving a slender beam or column. 
beams and columns compares favorably with con- = 
ventional techniques of numerical analysis, especially REFERENCES | main 
with regard to time consumption, mental effort, and 1 Southwell, R. V., Relaxation Methods in Theoretical Physics, vort 
accuracy. As compared with other types of analog 1st Ed., p. 19; Oxford University Press, London, 1946. | we 
computing equipment, the d.c. analyzer offers out- * Swenson, G. W., Jr., A Direct-Current Network Analyzer for | Vo 
standing advantages in economy and accuracy. Solving Wave-Equation Boundary. Value Problems, Ph.D. Thesis, " 
‘ : University of Wisconsin, 1951. | =m 
Although the examples presented involve only pin- ? Thomson, W. T., Critical Load of Columns of Varying Cross | abo 
ended members, any other boundary conditions can Section, Journal of Applied Mechanics, Vol. 2, No. 2, pp. 132-134, 
be accommodated as well. The analyzer is therefore June, 1950. | 
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RIEF REPORTS of investigations in the aeronautical sciences and discussions of papers pub- 


lished in the JOURNAL are presented in this special department. 
completed approximately 8 to 10 weeks after receipt of the material. 
not hold itself responsible for the opinions expressed by the correspondents. 


The publication will be 
The Editoral Committee does 
Contributions sohuld not 


exceed 800 words in length. 


Secondary Flow in Cascades 


Hsuan Yeh 

Department of Mechanical Engineering, The Johns Hopkins 
University, Baltimore 

December 20, 1951 


HE PROBLEM OF SECONDARY FLOW IN CASCADES received much 
‘aan recently because of its close relationship with the 
flow in a turbomachinery. In a recent article! in the Readers’ 
Forum of this JouRNAL an attempt was made to derive the sec- 
ondary flow in cascades from the principle of conservation of cir- 
culation. The interpretation of this principle (and, hence, the 
conclusion reached in reference 1), however, needs some clarifica- 
tion 

An approximate solution to this problem has been obtained by 
Squire and Winter? and Hawthorne.* However, the following 
shows that essentially the same results can be obtained in a much 


more direct manner. 
For incompressible nonviscous steady flows, the principle of 
permanence of vorticity can be expressed by* 


(Q-V)V 


(V-V)a = (1) 


dQ/dt = 

where Q is the vorticity vector and V is the velocity vector. 
If our interest is to find the variation along a stream line of the 
vorticity component which is parallel to the stream line (this 


vorticity component gives rise to the so-called secondary flow), 
Now, 


we seek to find the value of 7-V(Q- V7). 


Y=(P-vyat+(@avr+0xK(VxXa)+2Xx(9 X D) 


via- 7 


is a vector identity. Therefore, the dot product of V with the 
above gives 


V-V(a:- V7) = V-((V-V)a + (a°V)V) 


(2) 


Let — be the vorticity component parallel to the stream line, q¢ 


be the absolute magnitude of the velocity vector V(i.e., g = 
Hence, Eq. (2) can be written, with the 


(Q- V)/q. 
help of Eq. (1), as 


V|), then ¢ = 


i> a > 


V-V(éq) = 27-|(a-V)T) (3) 


Integrating along a stream line, we obtain, 

"2 7 [(Q- 
(tq) = 2 ds (4) 
1 q 


tq )o 
line. A 


However, 


where ds is an elementary distance along the stream 


similar expression appears as Eq. (10) of reference 3. 
in evaluating the integral of Eq. (10) of reference 3, one experi 
ences some difficulty when there is a stagnation point along the 
stream line. This difficulty seems to have been removed when 
the above equation |Eq. (4)] is used instead, since the integrand 
of Eq. (4) remains finite when g goes to zero. Although this shows 


that £, contrary to some earlier theories, does not become infinite 


downstream of a stagnation point, the calcuiation is not made any 
easier because the linearization process (below) fails when € is 
not very small. 

To obtain a solution to Eq. (4), it is desirable that some ap 
proximation be made. The most common one is the assumption 
that the value of — is so small that the main 7 and Q are not af- 
fected by it Ap- 
plying this to a cascade and approximating the passage between 


This is essentially a process of linearization 


any two blades by a curved channel as shown by Fig. 1, it is seen 
that (Q-V)V or the variation along the direction of @ (in this 
case along r) of the velocity vector is 


(a:-V)V = \al\(dV/dr) 


If we further assume that the turning is made in a passage of a 
circular are with an arc angle of e, that the flow in the circular 
passage is that of a free vortex, and that the channel is so narrow 
compared with the radius of turn that the latter is essentially the 
same (ro) for all stream lines, then, 


dV i d (“*) Vo 
dr dr 7 - ri 


Furthermore, since (Q:V)V is a vector along the direction of [, 


V-{((a-V)V] = g?2h/ro 
In Eq. (4) and for this case, (fg), = O and ds = rod@. Also, by 
this approximation, g and Q are both constant. Hence, 


= —2{) (5) 


fre 


which is also obtained in references, 2 and 3. 

Physically, it is interesting to note in connection with Eq. (5) 
that one-half of the secondary vorticity ~ is caused by the turning 
of the streamline by an angle ¢, assuming Q remains exactly con- 
stant (the component of 2 along the new velocity direction is 
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therefore £2). However, we also have the other half of — caused 
by the continuous reorientation of original vorticity, as clearly 
Since reference 1 ignores this second half of 


There is 


indicated by Eq. (1). 
£, the result therefore appears to be less representative. 
also some doubt as to the justification of any refinements based 
on the actual cascade geometry, unless more refined application 
of the exact Eq. (4) is made at the same time. 

Apart from the secondary vorticity discussed above, there is 
the additional secondary vorticity due to the shedding of trailing 
vortices behind the cascade blade, when the circulation along the 
blade span is nonuniform. The total induced flow is therefore 
the sum of two flows induced by the two kinds of vorticity. The 
two kinds of vorticity, however, are at different locations relative 
to the cascade. A simple algebraic addition of the two vortici- 
ties, disregarding their variation from point to point, appears in- 
adequate to prove that induced velocity is a minimum when the 
exit velocity is axial, as is claimed in reference 1. 
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Sweptback Wing Stress Analysis 


L. A. Goran 
McDonnell Aircraft Corporation, St. Louis, Mo. 
November 19, 1951 


theoretical and 


stress- 


_— RECENT PAPERS" ? have reviewed the 
literature on the sweptback 
On this subject, two McDonnell Aircraft 


experimental wing 
analysis problem. 
Corporation reports,®* 4 not discussed in these references, should 
be of general interest. 

In reference 3, test results and a minimum strain energy method 
of analysis are presented for a sweptback box beam having ribs 
parallel to the air stream. 

A summary of the work presented in reference 4 is given below: 

A minimum energy method of analysis is presented for the 
problem of the four-element sweptback box beam. The solution 
consists of a set of self-equilibrating ‘‘correcting’’ loads that are 
added to the loads calculated by ordinary bending theory 

The solution is extended to apply to the case of a multielement 
sweptback wing by assuming that the correcting loads have a 
linear stress distribution across the top and bottom cover plates. 
method of analysis is presented. 


An approximate ‘‘two-bay”’ 


In this method the structure outboard of the root rib is con- 
sidered to be of constant section, having an infinite number of 
closely spaced ribs, and a differential equation is solved to deter- 
mine the stress distribution in this bay. It is shown that good 
agreement exists at the root between the approximate method of 


analysis and the more laborious ‘‘exact’’ method. 
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Slender Body Theory for Supersonic Unsteag 
Flow , 


John W. Miles* 

Fulbright Lecturer, Auckland University College, Auckland New 
Zealand 

November 12, 1951 


oe HERE ARE THE RESULTS of extending Ward's 

slender body theory! to the case of unsteady flow, a question 
treated heuristically in a previous paper.?_ The complete analysis 
is to be presented in a forthcoming N.O.T-S. report, representing 
work supported by the Office of Naval Research 

In the following the body length (/) has been chosen as th, 
characteristic linear dimension and the free-stream velocity a 
time zero as the characteristic velocity (Ul). All subsequent 
quantities are dimensionless, and the characteristic, dimension. 
less parameters of the problem are 6, the larger of the transverse 
MW, the free. 


time rate of 


body dimension or amplitude of transverse motion; 
stream Mach Number; 
change, which for harmonic motion reduces to the well-known 


and k, the measure of 
reduced frequency. As coordinates, the Cartesian set s and : 
= x + 1y were chosen, where s is measured downstream from the 
body nose. 

The restrictions on body shape are identical with those given 
by Ward. 


follow the procedure applied by Lin, Reissner, and Tsien‘ to the 


To investigate the extent of the linearized régime, we 


two-dimensional problem and find the essential restrictions en- 
suring small perturbations to be 


k = O(1): 6< 1 and (Mb)? < 1 la 


k>1: ki<X 1 and (kM,5)?< 1 lb 
It is found that the velocity potential satisfies Laplace's equa 


tion in planes transverse to the flow and is given by the real part 


of the complex potential 


WS, 2, 


t) = ao(s, t) Ins + do(s, t) + >» Guls, t)s 
l 


| 
0 ( Jaren 5 2 
bk | 


subject to the linearized boundary condition 


(;) : : 
+ O 6? In 6 3 
k 


where » denotes the normal to the body contour in a transvers¢ 
plane, (D/Dt) implies 


Dn 
~ Dt 


Q, 


D( j Ol ) of 
l 4 
Dt Os ol 


and u(t) is the increment of flight velocity acquired in the interval 


t. It suffices to invoke this boundary condition at the mean 


The pressure is given by 


[/1 
4\¢2|* +0 ( ) ae Inz6| (5 
kA 


position of the body. 


(’ Po 


where fp and g 


Do 
2 
Dt 


disturbed flow. In each of the error terms, the larger of the fac 


tors 1 or k” is implied 
By an analysis similar to that given by Ward, it is found that 


1 D 


ads. t) S(s.t 6 


isi 2r Dt *’ 


from University of California, Los Angeles; formerly Con 


Naval Ordnance Test Station, Inyokern, Calif 


* On leave 


sultant, 


are the static and dynamic pressures in the un- 
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where S(s, ¢) is the transverse cross section, the results being valid 
for both time and spatial variations thereof. The nonbase drag 


js given by 
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where the line integrals are taken around the body contour at 


station s. For the transverse force and rolling moment per unit 
length, we find 
D? 


1 
= 1 2 2g) 
as\ go ta De? 


1 
0 I(;.) M244 In? | (9) 


rdy)] 4 
1 M.8* In? 
0 Re M,76° In? 6 (10) 


where We 
remark that the contribution of the terms dp In z and hp to the 


_— ; f ix + 
= -2 *(xdx + 
he “DI @*(xadx 4 


a9 


zy denotes the centroid of the cross section at s. 


last integral are null. 
The foregoing results have been applied to the questions of the 
acceleration and transverse motion of winged bodies of revolution 


In the latter case we find 


2 (* + iY wee Sy psy) 2 pom 
= — 27 bh? — a? + (a/b?) Sof 
ds qo Dt | ptf 
l 1 2 sin 4¢ 
N — = plb + (a2/b)}*| & + (1 — ~le+ 
% 4 1 2r 
(5) (1 a } I sin? 2 12 
— cos‘ ¢) - sin + sin* 2 (12 
' , { °° i si , 
¢ = cos”! [2ab f(a? + b?)] (13) 


where a is the body radius and 26 the total wing span (tip to tip), 
evaluated at station s in Equation (11) and at the trailing edge 
in Eq. (12), and p is the angular velocity in roll. These parti 
cular results are subject to the restriction 

b’(s) > 0 (14) 
and, in the case of Eq. (11), to the requirement that the body 
either terminate at the trailing edge or be of constant cross sec- 
tion thereafter. The restriction (14), which ensures no wake ef- 
fects, is more severe than that required in steady flow, where 
b(s) may be negative provided that the body be of constant 
cross section aft of the maximum wing span. 
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Comments on ‘‘The Evaluation of Theoretical 
Critical Compression in Sandwich Plates’’ 


Bernard Budiansky and Charles Libove 
Langley Aeronautical Laboratory, N.A.C.A 
November 12, 1951 


| aggro IN A RECENT PAPER by Goodier and Neou! may 
be misleading. The footnote that wide 
strip must begin to bend with anticlastic curvature 
load for infinitesimal buckling would appear, therefore, not to 
It is indeed true, as shown by 


states “since a 


the critical 


require the factor 1/(1 — v?).” 
Ashwell,? that in pure bending of a wide plate there is a gradual 
transition due to finite deflections from E to E/(1 — v*) in the 
effective Young’s modulus. However, the initial buckling of a 
plate column of finite width supported along its loaded edges 
always occurs with some inhibition of anticlastic curvature; 
Houbolt and Stowell’ have shown that, as the plate width varies 
from zero to infinity, the initial buckling load ranges from the 
column value (based on £) to a value slightly below the full 
plate value [based on E/(1 — v?)]. For a width-height ratio of 
only 6 the initial buckling load of a simply supported aluminum 
plate column is already 99 per cent of the full plate value 

After buckling commences at the Houbolt-Stowell load, it is 
probably true that, for the buckle to grow, the applied load must 
be increased in order to overcome the stiffening that occurs as 
the anticlastic curvature is suppressed still more through the 
finite deflection mechanism discussed by Ashwell. The upper 
limit to this load would be the full plate value based on F./(1 y?). 
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Three-Dimensional Conformal 
Transformations 


Keeve M. Siegel 

Willow Run Research Center, University of Michigan 
Mich. 

November 26, 1951 


Ypsilanti, 


THE SUCCESS in the use of conformal transforma 


oe OF 

tions in two-space for ideal fluid flow, one would be interested 
a priori, in conformal transformations in three-space. The three- 
dimensional method of inversion is valid in electrostatics. How 


ever, it will be shown that this method is not valid in hydro 
dynamics. 

When one attempts to find an exact three-dimensional solution 
to a boundary-value problem, he is almost always limited to the 
contour surfaces of coordinate systems in which the variables in 
the equation of the motion can be separated. When Laplace's 
equation is the equation of the motion and a solution is found for 


, One ustu- 


such a contour surface (in electrostatics, for example 
ally uses Kelvin’s inversion theorem to find solutions for other 


bodies. Thus, from a known electrostatic solution for a cone, a 
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contour surface of a coordinate system for which Laplace’s 
equation is separable, a solution was obtained for an ogive! by 
the method of inversion. However, the ogive is not a contour 
surface of a coordinate system for which Laplace’s equation is 
separable. 

Lemma 1. The inversion transformation (sometimes called 
the reciprocal radii transformation) is the only nontrivial con- 
formal transformation in three-space.? 

Theorem 1. 
transformation which is applicable to boundary-value problems in 
fluid mechanics. 

By Lemma 1, to prove theorem 1, we need only discuss the 
inversion transformation. 


There is no nontrivial three-dimensional conformal 


Let us consider the equation for an 
ideal, incompressible, axially symmetric flow in cylindrical co- 
ordinates (7, x, @). 

Boundary-Value Problem I 


y = stream function 


oO? 1-2 o? 
E*y = = . es 0 
or? r Or Ox? 
Wevaluated at body = O 
Boundary-Value Problem II 
¢@ = velocity potential 
o* 1 oO oO? 
or? ror Ox? 
fa) 
° = () 


O normal / evaluated at body 


Boundary-value problem I cannot be solved by inversion be- 
cause the equation of the motion does not remain invariant under 
Boundary value problem II 
cannot be solved by inversion (even though Laplace’s equation 


the inversion transformation. 


would remain invariant) because it has a Neumann boundary 
condition.* 

In hydrodynamics, for all nontrivial cases, it can be shown 
that, when the equation of the motion is not Laplace’s, the 
equation does not remain invariant under the inversion trans- 
formation. When the three-dimensional Laplace’s equation does 
occur, it occurs with Neumann boundary conditions. 

In two-dimensional theory, the method of inversion becomes 
applicable in both fields because in two-space the stream-function 
equation becomes Laplace’s equation. 

In three-dimensional electrostatics, the method of inversion is 
valid because the Dirichlet boundary conditions, ¢ = 0, can exist 
at the same time as Laplace’s equation. 
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Erratum! 


N. Scholz 

Dr.-Ing., Institute of Fluid Mechanics, Technical University, Braun- 
schweig, Germany 

January 2, 1952 


| ipeaherateornegetg MY REPORT! obtains a misprint that makes 
unintelligible Eq. (10) of this note. This equation should 


read 


APRIL, 1952 


r= Uat X 


=?% 


age" sinh 


n=1 


TC x , 

ir —ir or nf me 

e +e > d, tan h pth 
2t ) j (“. 


2~** means that the summation has to be taken only over the odd 


terms. 
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Comments on ‘‘Horizontal Tail Loads in 
Abrupt Pull-Ups from Level Flight’’' 


R. M. Rosenberg 
Design Specialist, Boeing Airplane Company 
January 9, 1952 


ECKER’S PAPER ON TAIL LOADS! presents an approach to that 
D subject which will yield correct results of interest to the en- 
gineer. Since his results are integrals of the differential equations 
in hand, they will approximate to nearly the same degree the re- 
sponse of an airplane as his differential equations approximate the 
behavior of the system. 

Before discussing the method used by Decker, it would be well 
to clear up an unintentional misunderstanding. Decker implies 
in his “Introduction” that Kelley and Misall? have considered 
the same problem and have developed a method yielding solu- 
tions that are ‘‘not unique.’’ Certainly, Decker did not intend 
to say all that is suggested by this statement. Decker, as well as 
Kelley and Misall, has considered a linear problem of the sort 


Lix™, x("—-0, , 2, x) = F(t) (1) 


If such an equation possesses unique solutions in the accepted 
sense of this term, then it is not possible to construct correct 
methods yielding solutions that are not unique. Decker’s remark 
on lack of uniqueness is meant to say only that there is no unique 
F(t) that satisfies the condition that, at a particular time ti, 
#(ti) = #; (a constant), but no more. Clearly, this remark refers 
to the uniqueness of the equation and not to that of solutions. 

Decker’s mode of attack on the problem is an inverse method. 
It consists in prescribing for Eq. (1) the solution x = x*(t) where 
the * is used to indicate that the solution is picked arbitrarily 
(within suitable restrictions to be defined later). Assuming that 
Eq. (1) satisfies the requirements for uniqueness and existence of 
solutions and that x*(¢) is admissible as a solution, it follows 
that x«*(¢) is the unique solution to Eq. (1) provided 


F(t) Let), gtd... *, 2") 


F(t) is the forcing function in Eq. (1). We also identify it as the 
“object function,”’ because it is that function in Eq. (1) which is 
to be so adjusted that a prescribed solution satisfies the equation. 
On physical grounds, F(t) of Decker’s problem must be continu- 
ous everywhere. Therefore, in general, any x*(t) is admissible if 
it satisfies all prescribed initial and end conditions and if it pos- 
sesses continuous derivatives up to order n. 

Decker’s approach is extremely simple and effective when ap- 
plied to difficult problems because it changes the need for solving 
a differential equation to one of solving an algebraic equation for 
an object function. This algebraic equation can always be 
solved explicitly for any quantity of interest if the original dif- 
ferential equation is linear. Even when that equation is non- 
linear, one can often solve the corresponding algebraic equation 
explicitly for the object function,* but at times this may not be 
possible. 
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It should be pointed out that only a limited amount of informa- 
tion can be gathered by this method of finding equations that are 
satisfied by a prescribed solution. This limitation is brought 
about by the fact that only one particular solution, satisfying one 
particular set of boundary values of one particular equation, is 
found. This precludes, for instance, the possibility of answering 
the following question related to Decker’s problem: What are 
the conditions on the elevator motion F(t) such that a prescribed 
maximum load will not be exceeded at any time? A solution of 
the differential equation in terms of F(t) (say, a convolution in- 
tegral) can provide this answer; such methods as the one used 
by Decker can also be used at times to provide answers of this 
type if they are pooled with other phases of analysis.* 

In the past, the method used by Decker has been applied with 
good success to linear and nonlinear systems, antonomous or 
5 In autonomous systems the object function is usually a 


not.® 
function f(x) of the dependent variable. In these cases, the pre- 
scribed solution takes the form « = #*(x) rather than x = x*(t); 


all other arguments remain unchanged. 

In nonlinear problems, the inverse method is particularly 
powerful because it releases the investigator from the restriction 
that solutions must lie in the neighborhood of those to a related 
linear problem; thus, particular solutions to strongly nonlinear 
equations can be found easily, and solutions neighboring on these 
may still be found by perturbation methods or iterative proce- 
dures. Finally, valuable quantitative information regarding the 
behavior of a nonlinear system can be gathered from the corre- 
spondence between object functions and solutions, obtained when 
prescribed solutions are varied in a systematic fashion over wide 


intervals. 
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Assistant Professor 
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December 18, 1951 


of Aeronautical Engineering, University of 


HE PAPER BY DECKER! is an application of the inverse method 
in airplane dynamics. In this method the airplane motion is 
specified and the required control time histories are determined. 
It appears to have been used first in reference 2. When properly 
used, it is a powerful tool and can be applied to nonlinear, as 
well as linear, problems. It has not in the past received the at- 
tention it merits. Several problems, including the subject of 
the paper, have been, and are being, investigated by this ap- 
proach at the University of Toronto.* 4 

Insofar as tail loads are concerned, the paper is the fourth pub- 
lished solution by the inverse method. The first three are refer- 
ences 5, 6, and 3, the last two of these having appeared within 


the past year and a half 
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In applying the inverse method, it must be remembered that 
there are infinite numbers of ways in which the airplane may per- 
form a given type of maneuver (e.g., attain a given load factor in 
a given time). The success or failure of the method hinges en- 
tirely upon the choice of function that is to describe the maneuver. 

In Decker’s paper, Eq. (13) is the chosen function in this case. 
The function must, above all, represent a physically possible mo- 
tion of the airplane, consistent with the assumed initial conditions 
and the limitations imposed by the machine, the pilot, and the 
control system. It should, in addition, be in accord with the 
manner in which the maneuver is likely, or may be required, to be 
Thus, the solution obtained must be continuous in all dis- 
placements, velocities, and accelerations. (Discontinuous accel- 
erations will involve discontinuous control displacements.) Vio- 
lation of any of these conditions implies infinite forces, which 


done. 


cannot exist. 

In references 5 and 6 this state of affairs has been satisfactorily 
achieved, with the result that the maneuvers are realistic and 
check well with flight records. 

In the present paper, the choice of Eq. (13) is rather unfortu- 
nate in this respect. 

(1) The initial rate of displacement of the elevator is not zero, 
and thus the maneuver is not consistent with an initial state of 
steady unaccelerated flight. It entails an infinite acceleration of 
the control system at time zero; this is, of course, impossible. 

(2) The negative maximum of dé,/dt is too small compared 
with the positive maximum. This results in the maximum 6, 
occurring too late and in underestimation of the positive 6 and 
negative L;, which occur in the initial stages The last effects 
may be important (see references 6 and 7) 

(3) The tail load used by the author as the design value is that 
occurring at 2nazx,, even though it is still increasing. This artificial 
cutting-off of the solution is made necessary by the failure of 
Eq. (13) to represent this type of maneuver adequately. If the 
solution were continued in time with that equation, the maxi- 
mum tail load would exceed that determined by Decker and occur 


at a time later than ”,,a7.. Decker correctly notes that this is not 


in accord with experience. 
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Erratum 


P. P. Bijlaard 
Professor, Cornell University, Ithaca, N.Y. 
December 18, 1951 


N MY PAPER, “‘Analysis of the Elastic and Plastic Stability of 
I Sandwich Plates by the Method of Split Rigidities—II,”’ in 
the December, 1951, issue of the JoURNAL OF THE AERONAUTICAL 
ScIENCEs reference 1 incorrectly referred to Part I in the May, 
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1951, issue of the JouRNAL. This reference, as given in the orig- TABLE 1 1B 
ane ‘ ‘TT » . 1 iila: > 2 y Thy > ‘ . — 6 ° . me fiom 
~— een should read: . eee, B  P., Analysts of the Corrected Disc Frictional Moment Coefficient as Function of lo 
Elastic and Plastic Stability of Sandwich Plates by the Method of Relative Disc and Fluid Motions j - 
Split Rizidities, Preprint No. 259, Institute of the Aeronautical - 5 ” 
Sciences, pp. 1-45, January, 1950. Cu V R., or | 
»R2 
M ~ 
. V5 
+ rs) (1/2 )p(wR)?R: 
0 1.94 | Fre 
0.1 2.03 F 
0.25 2.24 
On the Flow Around a Rotating Disc in a TZ : 7 G.I 
Uniform Stream* 15 420 Air 
2.0 4.82 | 
Arthur N. Tifford and Sheng To Chu 4.0 6.04 Dec! 
: 4 : 6.0 8.32 
Ohio State University C 
January 8, 1952 T 
} 
ii aR | ; 
5 a S hiicht; — . : ™ ABLE 2 vice 
N REFERENCE I, Schlichting and Truckenbrodt summarized a I We 1) Radial Friction Coeffici F / 
‘ ” : vca orrected ) Kadtal Friction Coefficient as F “ ala. om 
“Pohlhausen-type”’ solution of the boundary layer on a rotat- catenins . ee Sane eee OS FORCE as aly sae 
‘ pie ; ia : : tive Disc and Fluid Motions I 
ing disc in a uniform stream. The existence of exact solutions ; 
7 : ~ : : / cri 
of the Navier-Stokes equations in two of the cases covered— c,;(radial)V R., Z 
ee ° * too 
namely, (1) the case of a negligible stream velocity far from the or 
. . . ° y 2 9\.9 met 
disc and (2) the case of zero disc rotation is (as stated) well a Tr mo + w*)r? 
Tr . . = 2 2\y2 tor} 
known. The existence of these solutions led us to suspect that w ("/2)pe(a? + w?)r v ‘l 
° ° 9 cille 
all the cases treated are amenable to exact mathematical analysis. : ' : re tt 
. 0. 05 ie 
And indeed, we have found them so. 0 25 1 14 Th 
. . . . . a 
The details of our analysis will be given in a later paper.f 05 1 38 . 
~ . ° que 
Some of the numerical results, however, are discussed below. 1.0 1.88 ' I 
Table 1 lists the values obtained for the corrected frictional 4 ° 2 17 van 
ae : . : : 2.0 2.33 
moment coefficient as a function of the relative fluid and disc 10 ra ~ pro 
motions. As a partial check on precision, it may be noted that, 6.0 2 57 “‘- 
in the case of negligible free-stream velocity—i.e., a/w = O0— o% 2.61 ae 
the moment coefficient based upon Cochran’s calculations,? — ? ne 
agrees within 0.5 per cent with our value. res 
In Table 2 are given values of the corrected (surface) radial It appears probable, then, that the frictional coefficients of wh 
shearing stress coefficient as a function of the relative fluid Tables 1 and 2 are correct to within at least 1 per cent. On the ae 
and disc motions. Checking again with Cochran’s calculations other hand, the frictional moment coefficient values given by oi 
(for the case a/w = 0), agreement within 1 per cent is found. Schlichting and Truckenbrodt are about 7 per cent lower than pla 
Homann’s analysis* of the stationary disc—i.e., a/w = ©, pro- those of Table 1. (Values of the radial shearing stress coefficient ¥ 
vides an addition check. Making use of Hartree’s more precise were not reported.) These discrepancies are presumably mainly 
evaluation of Homann’s function—available as the case, B = due to the relatively approximate nature of their analysis in 
0.5, in reference 4—we find agreement with present results within As might be anticipated, the same coordinate system that leads eq 
0.4 per cent. to the exact solution of the Navier-Stokes equations, leads to the 
* Work supported by the Flight Research Laboratory, U.S.A.F. The exact solution of the seitithecst equatien. Thus, the analyses . 
numerical results of Tables 1 and 2 were obtained by the Computation Sec- references 5, 6, and 7{ may be unified into a single mathematical 
tion of the Flight Research Laboratory, U.S.A.F. treatment that simultaneously covers all intermediate flow geom- wh 
t It has just been brought to our attention that Batchelor has reported etries 
in reference 8 the existence of such an analysis as British A.R.C. Paper No. 
10,482, 1947, authored by Miss D. M. Hannah. 
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Frequency Response of an Airplane from 
Flight-Test Data 


G. M. Andrew and J. H. Polve 

Air Force Flight Test Center, Edwards Air Force Base, Muroc, 
Calif. 

December 1, 1951 


HE TREND IN MODERN AIRPLANES AND MISSILES toward using 
TF tion stabilized gun platforms, automatic landing de- 
vices, guidance systems, and all kinds of servomechanisms is be- 
coming increasingly more prevalent. 

The dynamic behavior of these complicated systems is de- 
scribed by the transfer functions. As the computations are not 
too accurate, the transfer functions are often determined experi- 
mentally. The transfer function can be determined in the labora- 
tory by putting a component, say the autopilot, on a table os- 
cillating with a known frequency and amplitude and recording 
the magnification factor and phase lag of the forced oscillation. 
The result covering a certain frequency range is called the fre- 
quency response. : 

If the autopilot is built into an airplane, the overall transfer 
function of the system (airplane plus autopilot) is given by the 
product of the airplane transfer function and the autopilot trans- 
fer function. The transfer function of the airplane must be de- 
termined from a flight test. Disturbing the steady-state flight 
by control surface motion (pulse or step) results in a dynamic 
response (angle, angular velocity, acceleration) of the airplane 
which is recorded vs. time. Because the experimental transfer 
function of the autopilot is obtained as a function of frequency, 
and that of the airplane as a time history, the response of the air- 
plane must be now translated in terms of frequency. When that 
is done, the product of both transfer functions can be obtained. 

The Laplace transformation relates the function f(t) of the 
time history to its transform f(s) in terms of frequency by the 


equation 
f(s) = f, e~“'f(t) dt (1) 
where 
s=oa+jw (2) 


Eq. (2) and ¢ = 0 give 


I, e Jot f(t) dt 


ex ® x 
= fi cos wif(t) dt — if sin wf(t) dt (3) 


f(jw) = 


0 


Generally, because a considerable number of time histories are 
recorded in a flight test, some automatic computing machine 
such as an IBM machine is needed to obtain results (say from a 
90 ordinate scheme) in a reasonable time. Where only slide 
rules and desk computers are available, the method of trapezoidal 
approximation, which requires a small number of ordinates, can 
be successfully used as a timesaving method. This method is 
used in reference 1 for the solution of the converse problem of 
getting time history from the frequency response. The deriva- 
tion of the method and the technique of computation are as fol- 
lows: Fig. 1 represents the time history of the airplane response 
The function f(t) is approximated by a few straight lines 
(f,a,b,c,d,¢ 


constant 


Point e corresponds to t = 7, where f(t) assumes a 


value f(7). Each trapezoid must be extended to the 


f(t) axis, like (f,a,b,0), (0,b,c,g), (0,e,d,g). 
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Let the areas above the 
If the area of the 


t axis have the plus sign; areas below, minus. 
trapezoid does not contribute to the area between the curve and 
Thus, the minus sign for 


From the nota- 


the ¢ axis, the sign must be reversed. 
(0,b,c,g) has to be reversed to plus for this reason. 
tion given in Fig. 1, it can be seen that 


i = (ta + ty)/2 
Ai = (ty _ $.) 2 
A = ryt 


In order to illustrate the evaluation of Eq. (3), consider the trape- 
zoid (f,a,b,0) for which the imaginary part becomes 


‘ a-a . : hitA 
Imy,a,b0) = -i ff sin (wt)r, dt — if, 
"a bi weil 
sin (wt )r; dt 


Proceeding similarly with the real part and summing all the trape- 
zoids, it can be demonstrated? that the transformation Eq. (3) 


sin t,w sin A,@ 7 fi. 
> a. I[ | -j> -x 
n=1 tw Aw n=1@ 
sin A,w 
E — ( ) cos 1 (4) 
A,@ 


The tables of the function sin x/x are given in reference 1, p. 357. 
When at time ¢ = 7, the function f(t) has a constant final value 
f(T) (Fig. 1), the end correction must be applied 


+ A; 
2A, 


yields 


f(jw) = 


f(T) sin wT S(T) cos wT 
j 


Af(jw) = — ; (5) 
aw WwW 
The real part of this correction is to be added to the real part (and 
the imaginary added to the imaginary part) of Eq. (4). 
For a satisfactory frequency response plot, it is enough to 
select from six to ten values of w for use in Eqs. (4) and (5). 
Comparing the results of this short method with the laborious 
one using, say, 90 ordinates, it will be found that the approxima- 


tion is surprisingly satisfactory. 
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On the Speed of Sound in Gases 


Irvine |. Glass 


Research Associate, Institute of Aerophysics, University of Toronto, 


Toronto, Canada 
November 4, 1951 


Fagen DISCUSSION of the speed of sound in fluids was conducted 
recently in these pages by Professors Truesdell, Murdochow, 
Morkovin, Millikan, and Munk.! 

In the case of one-dimensional nonlinear unsteady and inviscid 
flow such as through a rarefaction wave, the local speed of sound 
may be defined as the rate of propagation into the gas of a charac- 
teristic line (Mach wave) with respect to the flow velocity at the 
point under consideration.” This definition follows as a conse- 
quence of the solution of the nonlinear differential equations and 
does not depend on the linearized equations used in acoustics. 

This definition was tested experimentally in a wave interaction 
tube (shock tube )* by taking schlieren photographs with a wave 
speed camera of the rarefaction wave that is therein produced. 
In the theory of the shock tube‘ it is assumed that, when a dia- 
phragm, which separates a gas at high pressure from a gas at low 
pressure in a tube of uniform cross section, is suddenly ruptured, 
a shock wave and a contact surface are propagated into the low- 
pressure chamber (expansion chamber) and a centered rarefac- 
tion wave is propagated into the high-pressure chamber (com- 
pression chamber). The waves are assumed to form instantane- 
ously. It is also assumed that the waves are plane and travel at 
uniform velocities and separate regions of constant state. 

The theoretical (x, ¢)-plane of this wave model is shown in Fig. 
1. Along with this model are shown two schlieren photographs 
of the flow in the wave interaction tube, which were taken with 
the wave speed camera. The tube* has a 2.25 in. O.D. and a 1.00- 
in. 1.D. The flow is viewed through a lucite window 24 in. long 
and of a 4/j\-in. tee cross section. The wave speed camera has a 
rotating drum 12 in. in diameter. Seventy-millimeter film is 
used, and it attains a peripheral speed of 215 ft. per sec. The 
longitudinal striations on the photographic records are in main 


due to the lucite windows. 
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Fic. 1. The wave system generated in a shock tube. 
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The experimental results for the shock wave and contact sur. 
face are reported in reference 3. 

The present note deals with the rarefaction wave. 
tered rarefaction wave moving to the left (Fig. 1 


For a cen. 
), the following 
relations may be used: 


x/t=u-—a (1) 
u + [2a/(y — 1)] = 2ao/(y — 1) 9) 
Op 2po0 y-1/x (3 a 
= — ji- +1 (3 
Ox] t (y + 1) ao y + 1 \ac 
where 
a = local speed of sound 
a) = speed of sound in the rest state into which the rarefaction 
wave is advancing 
t = time 
x = distance 
u = flow velocity at a point (x, t) 
y = ratio of the specific heats C,/C, 
po = density of the rest state 


Since the rarefaction wave separates two constant states, the 
density derivative (Op/0x), is discontinuous at the head and at 
the tail of the rarefaction wave. As a result, these two border 
Mach lines should show up in a schlieren photograph of the (x, ¢)- 


plane. 
At the head of the wave mp = 0 and therefore x/t = —ay, the 
speed of sound of the gas in the rest state (the minus sign is con- 


sistent with Fig. 1). At the tail of the wave x/t = 1; — a, and, 
if w, is known, this relation would give a means of measuring the 
speed of sound in gases at extremely low temperatures and pres- 
sures. 

The speed of sound in the rest state is directly determined from 
the trace of the head of the rarefaction wave on the schlieren 
record, when the film speed and the optical magnification factor 
are known. 

Experimental results were obtained for the speeds of sound in 
air, argon, and carbon dioxide by direct measurement of the 
traces of the head of the rarefaction wave from a large number of 
schlieren photographic records, similar to the one that appears 
Distances on the photographic records could be meas 
Timing markers 


in Fig. 1. 
ured to 10-4 cm. with a traveling microscope. 
of 10-microsec. duration also appear on the photographic film 
This makes it possible to de 
The effects of 
A diaphragm 


and give a continuous time scale. 
termine the film speed to five significant figures 
film shrinkage must also be taken into account. 
pressure ratio range of 1.2 to 2,400 was utilized during the experi- 
ments. The sound speed in the rest state is independent of the 
pressure and of the diaphragm pressure ratio 

The following are the values obtained from the head of the 
rarefaction wave for the speed of sound at 0°C. in air, argon, and 
carbon dioxide. 

Air: ao = 1087 + 2 ft. per sec. 
1009 + 1 ft. per sec. 
844 + 1 ft. per sec. 


Argon: ao 
Carbon dioxide: a) = 


The above results are in good agreement with the theoretical 
sound speed obtained from the isentropic relation ay) = W yRTofor 
a perfect gas, or a9 = V yR7T)(K) for areal gas. K isa small cor- 
rection factor depending on the equation of state used for the gas 
The present results are also in good agreement with the values 
for ao obtained from acoustic methods.*~* For further details 
see reference 5 

In the present experiments the tail of the rarefaction wave did 
not appear on the schlieren records. Therefore, it was not pos- 
sible to measure the rate of propagation of this Mach wave. It is 
likely that the discontinuity in the density gradient tends to be 
smoothed out by viscous action where the gas is in motior 
Comments on this point are invited 
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On Interference Factors for Finned Bodies 


John W. Miles* 

Fulbright Lecturer, Auckland University College, Auckland, New 
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HE QUESTION OF WING-BODY INTERFERENCE for a_ slender 
j bad of revolution carrying a low aspect ratio wing has been 
treated by Ward.! We present here the extension of his results 
to the case of multiple fins. 
The conformal transformation 


2/2? + (a?, a yn . r* 2 +. (c?, c" 2 (1) 


maps the profile of the circle |z;| = a and the n, symmetrically 


disposed, radial segments bounded by the circles |2;| = a and 
| = bon the circle |¢| = c, where 
2ch/? = 5/2 + (a2/b)/? (2) 


The correct phases of the various radicals involved are ascer- 
tained by placing all branch cuts inside of the respective circles 
z| = aand ¢| = c¢ and by requiring the transformation to re- 
duce to zs; = ¢ at infinity, in order to preserve the flow there. 

Following the analysis of Ward, the transformation (1) en- 
ables us to determine the flow past a slender pointed body of 
revolution of radius a on which are mounted m isogonal fins of 
maximum total envelope 26, provided that we impose the re- 
striction that either the body radius be constant aft of the maxi- 
mum fin span or that both body and fins terminate at maxi- 
mum fin span. The resulting lift due to an angle of attack a is 
found to be 


24 2(1+6,' + 


fer 


L= 2rqh*a 


t= a/b (4) 


where g is the free-stream dynamic pressure and 6;' the Kronecker 
delta (= lift = 7; = Oif 7 # j). 
The lift due to the body alone is given by setting 6 = a, whence 


Le = 2rqa*a (5) 


* On leave from University of California, Los Angeles; formerly Con 


sultant, U.S. Naval Ordnance Test Station, Inyokern, Calif 
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Further, the lift due to (1/2) wings of span 2(b — a) in the ab- 
sence of interference is given by 


] “« 2 2rm 
Ley = (2rb?q)(1 — &) >” > cos? (6a) 
~ nm J n 


: 1 (1 + bn! + bp? 
= (2rb*ga)(1 — £) 5M - (6b) 


where the quantity in square brackets is the average incidence. 
We then define the interference factor 


«93 —(4/n) ny4 94 4 
L — Lg 2 (1 + ey)" — Qh Mes (7 
_ ‘> 
Lr n (1 — &)? 


nlé) = 
The first factor in Eq. (7)—namely, (2°~ ‘4/™ /n), may be thought 
of as a cascade (mutual fin) interference factor and the second 
as a fin-body interference factor, albeit this division is rather 
arbitrary. 

For very small wings, we have 

lim n,(£) = 4 (8) 
i 
which is a generalization of the result obtained by Ward. Thus, 
as might be expected, the extremely small fins act independently 
of one another and effectively have their aspect ratios doubled 
as a result of reflection in the body. 

A question of particular interest is the effectiveness of increas- 
ing the number of fins. For practical rockets the cases nm = 1 and 
2 are of no interest, since they provide forces that are dependent 
on the angle of roll. Moreover, the cruciform (n = 4) may be 
taken as a standard configuration, and we are led to define the 


aerodynamic effectiveness factor 
Ln — Le 

en(E) = = 21 — ¢*)-? 
L, — Lp 


The limits of interest are 


én(0) = 2'- 4/") (10) 
en(1) = n/4 (11) 


: 4\ /in2-—¢\" | ; 
éen(€) = Al — &)-'* 7 1 - = + O(n~*) |, 
n l— # 


se 21 (12) 


It would appear that, in virtue of representing a ratio of similar 
quantities, €,(£) would have a much wider range of validity than 
the slender body theory per se. 
We remark that the result obtained by allowing m to become 
infinite in Eq. (3), viz., 
lim Ly» = 2rqb*a(2 — é?*) (13) 
n> «x 
does not correspond to a body of revolution of radius 6, since 


the fins have been assumed to have zero cross section. 
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Comments on ‘‘Uniformly Loaded 
Semiinfinite Wedge-Shaped Plates’’ 


Eric Reissner 
Department of Mathematics, Massachusetts Institute of Technology, 


Cambridge,, Mass. 
December 19, 1951 


’ | SHE SOLUTION OF THE PROBLEM considered by Morgan! can be 
obtained a good deal more simply than the author has done 
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by using polar coordinates instead of cartesian coordinates. The 
differential equation 


Viw = P(x, y)/D [with P(x, y) = const. ] 
is reduced by the assumption w = r‘W(@) to an ordinary differen- 
tial equation with constant coefficients for W(@). Mr. Morgan’s 


solution can be written in the form 
w = r'cos‘ 6 f(tan™! @) 


and the function f then follows from a differential equation with 
variable coefficients. It would be difficult to see how Mr. Morgan 
recognized the possibility of solving this equation with variable 
coefficients as simply as it was done were it not for the known 
fact that the assumption w = r”W(@) reduces the expression V 4w 
to an expression LW(6@), where the operator L has constant coef- 


ficients. 
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Use of polar coordinates appears to remove Mr Morgan's re 
striction that wedge angles greater than 180° cannot be consid 
ered. It may also be seen that the same procedure is Possible for 
load functions of the form 

p(x, y) = r™Q(6) 
for which the proper choice for w is 
4W(0) 


w= ym 
w=Tf 


Concerning the problem of finite sector plates, reference should 
be made to a recent paper by Carrier and Shaw in the Proceedings 
of Symposia of Applied Mathematics, Vol. 3, p. 125, McGraw-Hill 
Book Company, Inc., New York, 1950 
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